Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEfficient computation of free energy of crystal phases due to external potentials by error-biased Bennett acceptance ratio method
148
0
0.0
(
0
)
Authors
Pankaj A. Apte
Ask ChatGPT about the research
No Arabic abstract
Free energy of crystal phases is commonly evaluated by thermodynamic integration (TDI) along a reversible path that involves an external potential. A persistent problem in this method is that a significant hysteresis is observed due to differences in the center of mass position of the crystal phase in the presence and absence of the external potential. To alleviate this hysteresis, a constraint on the translational degrees of freedom of the crystal phase is imposed along the path and subsequently a correction term is added to the free energy to account for such a constraint. In this work, we propose a new methodology termed as error-biased Bennett Acceptance ratio (EBAR) method that effectively solves this problem without the need to impose any constraint. This method is simple to implement as it does not require any modification to the path or to the simulation code. We show the applicability of this method in the computation of crystal-melt interfacial energy by cleaving wall method [J. Chem. Phys., 118, 7651 (2003)] and bulk crystal-melt free energy difference by constrained fluid $lambda$-integration method [J. Chem. Phys., 120, 2122 (2004)] for a model potential of silicon.
rate research
Read More
Direct correlation between temporal structural fluctuations and electron wind force is demonstrated, for the first time, by STM imaging and analysis of atomically-resolved motion on a thin film surface under large applied current (10e5 Amp/sqare cm). The magnitude of the momentum transfer between current carriers and atoms in the fluctuating structure is at least five to fifteen times (plus or minus one sigma range) larger than for freely diffusing adatoms. The corresponding changes in surface resistivity will contribute significant fluctuation signature to nanoscale electronic properties.
Very recently we developed an efficient method to calculate the free energy of 2D materials on substrates and achieved high calculation precision for graphene or $gamma$-graphyne on copper substrates. In the present work, the method was further confirmed to be accurate by molecular dynamic simulations of silicene on Ag substrate using empirical potential and was applied to predict the optimum conditions based on emph{ab initio} calculations for silicene growth on Ag (110) and Ag (111) surface, which are in good agreement with previous experimental observations.
A method based on the Gibbs adsorption isotherm is developed to calculate the decrease in interfacial free energy resulting from solute segregation at an internal interface, built on measured concentration profiles. Utilizing atom-probe tomography (APT), we first measure a concentration profile of the relative interfacial excess of solute atoms across an interface. To accomplish this we utilize a new method based on J. W. Cahns formalism for the calculation of the Gibbs interfacial excess. We also introduce a method to calculate the decrease in interfacial free energy that is caused by the segregating solute atoms. This method yields a discrete profile of the decrease in interfacial free energies, which takes into account the measured concentration profile and calculated Gibbsian excess profile. We demonstrate that this method can be used for both homo- and hetero-phase interfaces and takes into account the actual distribution of solute atoms across an interface as determined by APT. It is applied to the case of the semiconducting system PbTe-PbS 12 mol.%-Na 1 mol.%, where Na segregation at the PbS/PbTe interface is anticipated to reduce the interfacial free energy of the {100} facets. We also consider the case of the nickel-based Alloy 600, where B and Si segregation are suspected to impede inter-granular stress corrosion cracking (IGSCC) at homo- (GB) and hetero-phase metal carbide (M7C3) interfaces. The concentration profiles associated with internal interfaces are measured by APT using an ultraviolet (wavelength = 355 nm) laser to dissect nanotips on an atom-by-atom and atomic plane-by-plane basis.
By suitably adapting a recent approach [A. Laio and M. Parrinello, PNAS, 99, 12562 (2002)] we develop a powerful molecular dynamics method for the study of pressure-induced structural transformations. We use the edges of the simulation cell as collective variables. In the space of these variables we define a metadynamics that drives the system away from the local minimum towards a new crystal structure. In contrast to the Parrinello-Rahman method our approach shows no hysteresis and crystal structure transformations can occur at the equilibrium pressure. We illustrate the power of the method by studying the pressure-induced diamond to simple hexagonal phase transition in a model of silicon.
We introduce an approach to exploit the existence of multiple levels of description of a physical system to radically accelerate the determination of thermodynamic quantities. We first give a proof of principle of the method using two empirical interatomic potential functions. We then apply the technique to feed information from an interatomic potential into otherwise inaccessible quantum mechanical tight-binding calculations of the reconstruction of partial dislocations in silicon at finite temperature. With this approach, comprehensive ab initio studies at finite temperature will now be possible.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
