Machine learning for predictive condensed-phase simulation

We show how machine learning techniques based on Bayesian inference can be used to reach new levels of realism in the computer simulation of molecular materials, focusing here on water. We train our machine-learning algorithm using accurate, correlated quantum chemistry, and predict energies and forces in molecular aggregates ranging from clusters to solid and liquid phases. The widely used electronic-structure methods based on density-functional theory (DFT) give poor accuracy for molecular materials like water, and we show how our techniques can be used to generate systematically improvable corrections to DFT. The resulting corrected DFT scheme gives remarkably accurate predictions for the relative energies of small water clusters and of different ice structures, and greatly improves the description of the structure and dynamics of liquid water.

Computer simulation of a-Si : H/mu c-Si : H diphasic silicon solar cells

Based on our experimental research on diphasic silicon films, the parameters such as absorption coefficient, mobility lifetime product and bandgap were estimated by means of effective-medium theory. And then computer simulation of a-Si: H/mu c-Si: H diphasic thin film solar cells was performed. It was shown that the more crystalline fraction in the diphasic silicon films, the higher short circuit density, the lower open-circuit voltage and the lower efficiency. From the spectral response, we can see that the response in long wave region was improved significantly with increasing crystalline fraction in the silicon films. Taking Lambertian back refraction into account, the diphasic silicon films with 40%-50% crystalline fraction was considered to be the best intrinsic layer for the bottom solar cell in micromorph tandem.

Investigation of periodicity fluctuations in strained (GaNAs)(1)(GaAs)(m) superlattices by the kinematical simulation of x-ray diffraction

Periodicity fluctuations of layer thickness and composition in a superlattice not only decrease the intensity, they also broaden the width of the satellite peaks in the x-ray diffraction pattern. In this letter, we develop a method that is dependent on the width of satellite peaks to assess periodicity fluctuations of a superlattice quickly. A linear relation of the magnitude of fluctuations, peak width and peak order has been derived from x-ray diffraction kinematical theory. By means of this method, periodicity fluctuations in strained (GaNAs)(1)(GaAs)(m) superlattices grown on GaAs substrates by molecular beam epitaxy have been studied. Distinct satellite peaks indicate that the superlattices are of high quality. The N composition of 0.25 and its fluctuation of 20% in a strained GaNxAs1-x monolayer are obtained from simulations of the measured diffraction pattern. The x-ray simulations and in situ observation results of reflection high-energy electron diffraction are in good agreement. (C) 1999 American Institute of Physics. [S0003-6951(99)00828-1].

An analysis and simulation about conditioned launch for very short reach (VSR-4)

VSR4 links use graded index multimode fibers (GIMMFs) as the transmission medium with operation wavelength 850nm. For cost reasons, VCSEL has been selected as the optical source to VSR4. The minimum bandwidth specification for 62.5um GIMMF in VSR4 is only 400 MHz(.)km for over-filled-launch (OFL) condition. The distance of 300 meters is limited over transmission rates of 1.25Gbit/s on the basis of this specification. In order to overcome the OFL bandwidth limit by selective excitation of a limited number of modes, conditioned launch technique is investigated. In this paper, based on a comprehensive dispersion theory of GIMMF, a model is built to simulate the transmission of optical signal in GIMMFs and a comparison between OFL and conditioned launch is analyzed. The result can be the guidelines for the best choice of techniques for various LAN and interconnect systems also.

THE STUDY OF MICROSTRUCTURE IN IRON-DOPED AND SULFUR-DOPED INP CRYSTALS BY MEANS OF HREM AND COMPUTER-SIMULATION

The microstructures in iron- and sulphur-doped InP crystals were studied using both electron microscopy and electron diffraction. A modulated structure has been found in S-doped InP crystal, where the commensurate modulations corresponded to periodicities of 0.68 nm and 0.7 nm in real space and were related to the reflections of the cubic lattice in [111] and [113BAR] directions; they were indexed as q111* = 1/2(a* + b* + c*) and q113BAR* = 1/4(-a* - b* + 3c*), respectively. Single atomic layers of iron precipitate were observed, with preferred orientations along which precipitates are formed. Simulated calculations by means of the dynamical theory of electron diffraction using models for the precipitate structure were in good agreement with our experimental results. The relation between the modulated structure and the precipitates is also discussed.

Comparison of aggregation behaviors between branched and linear block polyethers: MesoDyn simulation study

The comparison of aggregation behaviors between the branched block polyether T1107 (polyether A) and linear polyether (EO)(60)(PO)(40)(EO)(60) (polyether B) in aqueous solution are investigated by the MesoDyn simulation. Polyether A forms micelles at lower concentration and has a smaller aggregation number than B. Both the polyethers show the time-dependent micellar growth behaviors. The spherical micelles appear and then change to rod-like micelles with time evolution in the 10 vol% solution of polyether A. The micellar cluster appears and changes to pseudo-spherical micelles with time evolution in the 20 vol% solution of polyether A. However, the spherical micelles appear and change to micellar cluster with time evolution in the 20 vol% polyether B solution. The shear can induce the micellar transition of both block polyethers. When the shear rate is 1x10(5) s(-1), the shear can induce the sphere-to-rod transition of both polyethers at the concentration of 10 and 20 vol%. When the shear rate is lower than 1x10(5) s(-1), the huge micelles and micellar clusters can be formed in the 10 and 20 vol% polyether A systems under the shear, while the huge micelles are formed and then disaggregated with the time evolution in the 20 vol% polyether B system.

Bump-Surface Multicompartment Micelles from a Linear ABC Triblock Copolymer: A Combination Study by Experiment and Computer Simulation

Novel bump-surface multicompartment micelles formed by a linear amphiphilic ABC triblock copolymer via self-assembly in selective solvent were successfully observed both in simulation and experiment. The results revealed that the block A forms the most inner core, and the blocks B and C form the inner and outer layers, respectively, and the bumps were formed by block A and more likely to be born on curving surfaces. Moreover, the micelle shape could be controlled by changing the solvent selectivity of the blocks A and B. Spherical, cylindrical, and discoidal micelles with bumpy surfaces were obtained both in experiment and simulation.

Monte Carlo simulation of the grafting of maleic anhydride onto polypropylene at higher temperature

The graft of maleic anhydride (MAH) onto isotactic polypropylene (iPP) initiated by dicumyl peroxide (DCP) at 190 degreesC was studied by means of the Monte Carlo method. The ceiling temperature theory, i.e., no possibility for the homopolymerization of MA-H to occur at higher temperatures, was used in this study. The simulation results show that most MAH monomers were grafted onto the radical chain ends arising from beta scission at a lower MAH concentration, whereas the amount of MAH monomers attached to the tertiary carbons was much larger than that grafted onto the radical chain ends at a higher MAH concentration for various DCP concentrations. This conclusion gives a good interpretation for the disagreement on the grafting sites along a PP chain. Moreover, it was found that the grafting degree increased considerably up to a peak value; thereafter, it decreased continuously with increasing MA-H concentration. The peak shifted in the lower MAH concentration direction and became lower and lower with increasing DCP concentration. When the DCP concentration was below 0.1 wt %, the peak was hardly observed. Those results are in good agreement with the experiments.

Simulation of characteristic electrochemical behavior of azobenzene derivative self-assembled monolayers

The characteristic electrochemical mechanics of azobenzene derivative self-assembled monolayers is discussed in present paper. It is presented that the structure inhibition is one of the most important factors in the increase of electrochemical reactive energy. A corresponding mathematical model was established based on Levich and Marcus's theory. Moreover, computational program was written to simulate the decrease of apparent rate constant (k(app)) of electron transfer with increasing surface concentration.

Monte Carlo simulation of phase behavior of polymer blends with special interactions

Full Paper: The phase, behavior of A-B-random copolymer/C-homopolymer, blends with special interaction was studied by a. Monte, Carlo simulation in two dimensions. The interaction between I segment A and segment C was repulsive, whereas it was attractive between segment B and segment C. The simulation results showed that the blend became two large co-continuous phase domains at lower segment-B component compositions, indicating that the blend showed spinodal decomposition. With an increase of the segment-B component, the miscibility between the copolymer,and the polymer was gradually improved up to being miscible. In addition, it was found that segment B tended to move to the surface of the copolymer phase in the case of a lower component of segment B. On the other hand, if was observed that the average, end-to-end distances ((h) over bar) for both copolymer and polymer changed slowly with increasing segment-B component of the copolymer up to 40%, thereafter they increased considerably with increasing segment B component. Moreover, it was found that the (h) over bar of the copolymer was obviously shorter than that of the homopolymer for the segment-B composition, region from 0% to 80%. Finally, a, phase diagram showing I phase and - II phase regions under the condition of constant-temperature is presented.

Molecular dynamics simulation of sub-transition for polyethersulfone

Molecular dynamics (MD) simulations of a polyethersulfone (PES) chain are carried out in the amorphous state by using the Dreiding 2.21 force field at four temperatures. Two types of molecular motion, i.e, rotations of phenylene rings and torsions of large segments containing two oxygen atoms, two sulfur atoms, and five phenylene rings on the backbone, are simulated. The modeling results show that the successive phenylene rings should be in-phase cooperative rotations, whereas the successive large segments should be out-of-phase cooperative torsions. By calculating the diffusion coefficient for the phenylene ring rotations, it is found that this rotation contributes to the beta -transition of PES.

Monte Carlo simulation of phase separation of A/B/A-B ternary mixtures

Monte Carlo simulations were used to model A/B/A-B ternary mixtures with different AB diblock copolymer volume fractions for which both the dispersed and continuous phase volume fractions were kept constant. For concentrations of the diblock copolymer below a critical value, the domain size increment of the dispersed phase decreases linearly with the copolymer concentration. This is in agreement with the predictions of Noolandi and Hong. The dependence of the domain size as a function of the copolymer volume fraction can also be fitted by the equation of Tang and Huang. Our simulations indicate, for the first time, that the micelles form before saturation of the interface occurs. This means that the formation of the micelles is not a result of the saturation of the interface.

Molecular dynamics simulation of polystyrene-block-poly(methyl methacrylate)

Molecular dynamics is applied to the system of polystyrene-block-poly(methyl methacrylate). The simulation shows that for the block copolymer system, a layered structure, which reflects microphase separation, is obtained and this structure is stable. In order to elucidate that the formation of the layered structure is reasonable, some static properties such as the radial distribution function and the dipole moment are analyzed in some detail.

Numerical Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the source region, the seafloor lifts to a designated height by a generation function. The numerical tests show that file linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations eta(max)(0) are carried out by both the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on eta(max)(0) are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences between the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the eta(max)(0) near-linearly varies with the wave amplitudes of the surface waves, and the eta(max)(0) has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, aid these differences are significantly affected by The wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.