Simulation population kinetics

The motion of particles of the film and substrate were calculated by standard molecular dynamics techniques. In the simulations discussed here, our purpose is to calculate equilibrium or metastable configurations of the system at zero Kelvin. For this purpose, we have applied random and dissipative forces to the particles. Finite random forces provide the thermal motion which allows the system to explore different configurations, and the dissipation serves to stabilize the system at a fixed temperature. The potential energy minima are populated by reducing the random forces to zero, thus permitting the dissipation to absorb the kinetic energy. 

Figure 2. Kinetic schemes for Na channel gating (right) and the graphed time-course for single channels (solid lines, the higher position is open ) and for the population of many channels (broken line, the fraction open increases upwardly). Numbers at the arrows of the kinetic scheme are the rate constants, in 10 sec" The period of simulation is 5 msec. Computerized model courtesy of Dr. Daniel Chemoff.

The determination of the laser-generated populations rij t) is infinitely more delicate. Computer simulations can certainly be applied to study population relaxation times of different electronic states. However, such simulations are no longer completely classical. Semiclassical simulations have been invented for that purpose, and the methods such as surface hopping were proposed. Unfortunately, they have not yet been employed in the present context. Laser spectroscopic data are used instead the decay of the excited state populations is written n (t) = exp(—t/r ), where Xj is the experimentally determined population relaxation time. The laws of chemical kinetics may also be used when necessary. Proceeding in this way, the rapidly varying component of AS q, t) can be determined. 

The resin system selected to initiate these studies is a step-growth anhydride cured epoxy. The approach to the kinetic analysis is that which is prevalent in the chemical engineering literature on reactor design and analysis. Numerical simulations of oligomeric population density distributions approximate experimental data during the early stages of the cure. Future research will... 

Self-consistent analysis of the 1-D fluid and kinetic models presented here can be developed through numerical simulations, where (17.30) and (17.32) or (17.36) and (17.37) can be used to define the initial conditions for the electric field at time 7 = 0, when the ion acceleration process begins. Results from these studies that use the fluid and the kinetic descriptions can be found in the literature, where either a single electron population [88,96-99]... 

---