Simulation with Empirical Models

The observed empirical models should now be employed to simulate and predict kinetic behaviors obtained with administration protocols other than that used for observation. Moreover, we must develop pharmacokinetics in a multicompartment system by including the presence of a fractal organ. We have argued that the liver, where most of the enzymatic processes of drug elimination take place, has a fractal structure. Hence, we expect transport processes as well as chemical reactions taking place in the liver to carry a signature of its fractality. 

Little has so far been done to predict the effect of different modes of administration, according to inhomogeneous conditions, on the observed c(t) when this contains a power function. In fact, the availability of the drug in the process was simply expressed by an initial condition c (to) = cq. Later on, exponential, power-law, or gamma profiles were observed according to the inherent heterogeneity of the process. 

---

Rene Fournier is studying atomic clusters238 and transition metal complexes.239 He is using a combination of density functional methods, tight-binding models, and molecular simulations with empirical interaction potentials, as part of a research program designed to study materials by computations on simple model systems. 

Nonequilibrium sorption due to mass-transfer limitations (including slow external or internal diffusion) and sorption to two different sorbents have been incorporated into a single ADE to evaluate the conditions under which mass-transfer processes may be important [206]. Simulations with this model, using mass-transfer parameters estimated from empirical correlations, reveal nonequilibrium conditions (i.e., mass-transfer limitations) when groundwater velocities increase (such as those that might occur in a funnel-and-gate system). 

Molecular modeling and computer simulation with empirical potential energy function (force field) are now routinely carried out to help understand and predict structures and dynamics of proteins and other macromolecules of biological relevance in water and membrane environments. After over 40 years of development, popular force fields such as AMBER, CHARMM, OPLS and GROMOS have been widely employed in biomolecular simulations. These force fields are used dominantly in highly optimized molecular dynamics... 

Due to the development of advanced numerical methods in the last decades, quantum approaches are now able to accurately describe the chemical bonds formed between two reactants. Nevertheless, when a surface is involved, the actual systems met in practice, for example a dense polymeric layer adsorbed on a rough surface, cannot yet be simulated, because this would require too large a memory size or too long a computation time. Quantum calculations, thus, cannot compete with empirical models in the prediction of adhesion strengths. However, they may allow one to check their validity in model cases, for example small molecules adsorbed on a substrate, or large molecules adsorbed on a cluster of a few atoms which simulates the substrate. This has been done in a number of cases but, to the author s knowledge, mostly for adsorption processes on metallic surfaces. Numerical results for the adsorption of molecules on oxide surfaces may be found in the literature (Henrich and Cox, 1994), but there exists no systematic discussion in the framework of acid-base interactions. 

There are many large molecules whose mteractions we have little hope of detemiining in detail. In these cases we turn to models based on simple mathematical representations of the interaction potential with empirically detemiined parameters. Even for smaller molecules where a detailed interaction potential has been obtained by an ab initio calculation or by a numerical inversion of experimental data, it is usefid to fit the calculated points to a functional fomi which then serves as a computationally inexpensive interpolation and extrapolation tool for use in fiirtlier work such as molecular simulation studies or predictive scattering computations. There are a very large number of such models in use, and only a small sample is considered here. The most frequently used simple spherical models are described in section Al.5.5.1 and some of the more common elaborate models are discussed in section A 1.5.5.2. section Al.5.5.3 and section Al.5.5.4. 

Means of zonal wind and temperature, simulated with 3-D ozone, after 70 days of model simulations, that corresponds to the end of February, are plotted in Figures 6.a) and 6.b). The values and structure of zonal wind and temperature are well comparable with the empirical model CIRA-86 [1], A comparison of corresponded Figures 5.a),b) and 6a),b) exhibit characteristic features of annual variability of zonal wind and... 

One of the concerns regarding the use of COSMO-RS thermodynamics directly in simulations is the considerably larger computation time that is required for the evaluation of the activity coefficients compared to simpler empiricEd models with... 

As a proof of the feasibility of such direct COSMO-RS process simulation, Taylor et al. [100] have linked the COSMOtherm program into their simulation program CHEMSEP [101] for distillation separation processes. For a number of typical separation problems they report very satisfying results, which are comparable with simulations based on empirical models. The simulation times were only a factor of 2 greater than those using empirical models. The quality of the simulations was considered as comparable to empirical models, although those were based on fitted experimental data. 

Bostrom—Kunz—Ninham (BKN) Model for Ion-Solvation Forces. A new model for the short-ranged ion-solvation forces has been proposed recently.12 It is based on the observation that the air/water interface is not sharp, the water density increasing gradually, over a distance of a few angstroms, from zero (in air) to the density of bulk water. The water density profile was obtained by fitting the results of the molecular dynamics simulations with the empirical expression12... 

Future directions in the development of polarizable models and simulation algorithms are sure to include the combination of classical or semiempir-ical polarizable models with fully quantum mechanical simulations, and with empirical reactive potentials. The increasingly frequent application of Car-Parrinello ab initio simulations methods " may also influence the development of potential models by providing additional data for the validation of models, perhaps most importantly in terms of the importance of various interactions (e.g., polarizability, charge transfer, partially covalent hydrogen bonds, lone-pair-type interactions). It is also likely that we will see continued work toward better coupling of charge-transfer models (i.e., EE and semiem-pirical models) with purely local models of polarization (polarizable dipole and shell models). 

Molecular modeling applied to polymers is, in principle, an extension of the concepts applied to small molecules. Readers familiar with energy minimization with empirical force fields (MM2, etc.), Monte Carlo, or molecular dynamics techniques for simulations already know a significant part of what is required to model polymer systems. What we present here is a description of the various techniques available for the simulation of polymers. The discussion is mostly limited to homopolymers, although we briefly mention some exciting topics outside this area. 

---