Atomistic phase

FIG. 8 Sketch of the four types of degrees of freedom used in the atomistic-continuum model H, Dh, and the scaled coordinates and s of the nodes in the continuum and the atoms in the atomistic phase. The system unit cell matrix H and the scaled coordinates and s entirely determine the finite-element configuration of the continuum, h and s are necessary to determine the atomistic configuration. The nodal points on the inclusion boundary s and the atomistic scaling matrix h are calculated from the variables H and Dh. 

Since and depend only on die valence charge densities, they can be detennined once the valence pseudo- wavefiinctions are known. Because the pseudo-wavefiinctions are nodeless, the resulting pseudopotential is well defined despite the last temi in equation Al.3.78. Once the pseudopotential has been constructed from the atom, it can be transferred to the condensed matter system of interest. For example, the ionic pseudopotential defined by equation Al.3.78 from an atomistic calculation can be transferred to condensed matter phases without any significant loss of accuracy. 

The microstmcture and imperfection content of coatings produced by atomistic deposition processes can be varied over a very wide range to produce stmctures and properties similar to or totally different from bulk processed materials. In the latter case, the deposited materials may have high intrinsic stress, high point-defect concentration, extremely fine grain size, oriented microstmcture, metastable phases, incorporated impurities, and macro-and microporosity. AH of these may affect the physical, chemical, and mechanical properties of the coating. 

In addition to the MD method, a wealth of Monte Carlo methods is used also at the atomistic level [6]. They use essentially the same models, force fields, for polymers. Their main advantage, however, is that by introduction of clever moves one can beat the slow physical dynamics of the systems and can run through phase space much faster than by MD. These methods are still in their infancy, but will certainly become more important. 

Short-time Brownian motion was simulated and compared with experiments [108]. The structural evolution and dynamics [109] and the translational and bond-orientational order [110] were simulated with Brownian dynamics (BD) for dense binary colloidal mixtures. The short-time dynamics was investigated through the velocity autocorrelation function [111] and an algebraic decay of velocity fluctuation in a confined liquid was found [112]. Dissipative particle dynamics [113] is an attempt to bridge the gap between atomistic and mesoscopic simulation. Colloidal adsorption was simulated with BD [114]. The hydrodynamic forces, usually friction forces, are found to be able to enhance the self-diffusion of colloidal particles [115]. A novel MC approach to the dynamics of fluids was proposed in Ref. 116. Spinodal decomposition [117] in binary fluids was simulated. BD simulations for hard spherocylinders in the isotropic [118] and in the nematic phase [119] were done. A two-site Yukawa system [120] was studied with... 

Chemical vapor deposition may be defined as the deposition of a solid on a heated surface from a chemical reaction in the vapor phase. It belongs to the class of vapor-transfer processes which is atomistic in nature, that is the deposition species are atoms or molecules or a combination ofthese. Beside CVD, they include various physical-vapor-deposition processes (PVD) such as evaporation, sputtering, molecular-beam epitaxy, and ion plating. 

This article reviews progress in the field of atomistic simulation of liquid crystal systems. The first part of the article provides an introduction to molecular force fields and the main simulation methods commonly used for liquid crystal systems molecular mechanics, Monte Carlo and molecular dynamics. The usefulness of these three techniques is highlighted and some of the problems associated with the use of these methods for modelling liquid crystals are discussed. The main section of the article reviews some of the recent science that has arisen out of the use of these modelling techniques. The importance of the nematic mean field and its influence on molecular structure is discussed. The preferred ordering of liquid crystal molecules at surfaces is examined, along with the results from simulation studies of bilayers and bulk liquid crystal phases. The article also discusses some of the limitations of current work and points to likely developments over the next few years. 

The question arises as to how useful atomistic models may be in predicting the phase behaviour of real liquid crystal molecules. There is some evidence that atomistic models may be quite promising in this respect. For instance, in constant pressure simulations of CCH5 [25, 26] stable nematic and isotropic phases are seen at the right temperatures, even though the simulations of up to 700 ps are too short to observe spontaneous formation of the nematic phase from the isotropic liquid. However, at the present time one must conclude that atomistic models can only be expected to provide qualitative data about individual systems rather than quantitative predictions of phase transition temperatures. Such predictions must await simulations on larger systems, where the system size dependency has been eliminated, and where constant... 

The rapid rise in computer speed over recent years has led to atom-based simulations of liquid crystals becoming an important new area of research. Molecular mechanics and Monte Carlo studies of isolated liquid crystal molecules are now routine. However, care must be taken to model properly the influence of a nematic mean field if information about molecular structure in a mesophase is required. The current state-of-the-art consists of studies of (in the order of) 100 molecules in the bulk, in contact with a surface, or in a bilayer in contact with a solvent. Current simulation times can extend to around 10 ns and are sufficient to observe the growth of mesophases from an isotropic liquid. The results from a number of studies look very promising, and a wealth of structural and dynamic data now exists for bulk phases, monolayers and bilayers. Continued development of force fields for liquid crystals will be particularly important in the next few years, and particular emphasis must be placed on the development of all-atom force fields that are able to reproduce liquid phase densities for small molecules. Without these it will be difficult to obtain accurate phase transition temperatures. It will also be necessary to extend atomistic models to several thousand molecules to remove major system size effects which are present in all current work. This will be greatly facilitated by modern parallel simulation methods that allow molecular dynamics simulations to be carried out in parallel on multi-processor systems [115]. 

As computer power continues to increase over the next few years, there can be real hope that atomistic simulations will have major uses in the prediction of phases, phase transition temperatures, and key material properties such as diffusion coefficients, elastic constants, viscosities and the details of surface adsorption. 

Two macromolecular computational problems are considered (i) the atomistic modeling of bulk condensed polymer phases and their inherent non-vectorizability, and (ii) the determination of the partition coefficient of polymer chains between bulk solution and cylindrical pores. In connection with the atomistic modeling problem, an algorithm is introduced and discussed (Modified Superbox Algorithm) for the efficient determination of significantly interacting atom pairs in systems with spatially periodic boundaries of the shape of a general parallelepiped (triclinic systems). 

As for any modeling of continuum structures, the properties of the phases must be known for this kind of approach to work. Here, estimates obtained by atomistic methods of other techniques, described in the earlier chapters of this review, may be employed, or empirically known values may be used. It is hoped that the co-development of these continuum techniques and atomistic and coarse-grained approaches will lead to a seamless integration of the different techniques. 

Clearly, the total number of unknowns that need to be determined is m = a + +. .. + z and a solution set for parameters p, P2 pm is determined using the singular value decomposition or any other suitable method. The mean pair-force corresponding to the potential of mean force can be obtained in a systematic manner by averaging a number of sets of solutions for parameters p, P2 Pm obtained along the atomistic MD trajectory in which the phase space is sampled extensively. 

Even when the composition range of a nonstoichiometric phase remains small, complex defect structures can occur. Both atomistic simulations and quantum mechanical calculations suggest that point defects tend to cluster. In many systems isolated point defects have been replaced by aggregates of point defects with a well-defined structure. These materials therefore contain a population of volume defects. 

End-Bridging Monte Carlo A Fast Algorithm for Atomistic Simulation of Condensed Phases of Long Polymer Chains. 

Chemical equilibrium methods provide useful predictions of the EOS of detonation processes and the product molecules formed, but no details of the atomistic mechanisms in the detonation are revealed. We now discuss condensed-phase detonation simulations using atomistic modeling techniques to evaluate reaction mechanisms on the microscopic level. 

Gas-phase results provide insight into the reaction pathways for isolated HE molecules however, the absence of the condensed-phase environment is believed to affect reaction pathways strongly. Some key questions related to condensed-phase decomposition are as follows (1) How do the temperature and pressure affect the reaction pathways (2) Are there temperature or pressure-induced phase-transitions that play a role in the reaction pathways that may occur (3) What happens to the reaction profiles in a shock-induced detonation These questions can be answered with condensed-phase simulations, but such simulations would require large-scale reactive chemical systems consisting of thousands of atoms. Here we present results of condensed-phase atomistic simulations, which are pushing the envelope toward reaching the required simulation goal. 

Atomistic simulations have been performed on condensed-phase HMX, which is a material that is widely used as an ingredient in various explosives and propellants. A molecular solid at standard state, it has four known... 

---