Abstract: In order to increase the accuracy of the linearized augmented plane wave (LAPW) method, we present a new approach where the plane wave basis function is augmented by two different atomic radial components constructed at two different linearization energies corresponding to two different electron bands (or energy windows). We demonstrate that this case can be reduced to the standard treatment within the LAPW paradigm where the usual basis set is enriched by the basis functions of the tight binding type, which go to zero with zero derivative at the sphere boundary. We show that the task is closely related with the problem of extended core states which is currently solved by applying the LAPW method with local orbitals (LAPW+LO). In comparison with LAPW+LO, the number of supplemented basis functions in our approach is doubled, which opens up a new channel for the extension of the LAPW and LAPW+LO basis sets. The appearance of new supplemented basis functions absent in the LAPW+LO treatment is closely related with the existence of the ul-component in the canonical LAPW method. We discuss properties of additional tight binding basis functions and apply the extended basis set for computation of electron energy bands of lanthanum (face and body centered structures) and hexagonal close packed lattice of cadmium. We demonstrate that the new treatment gives lower total energies in comparison with both canonical LAPW and LAPW+LO, with the energy difference more pronounced for intermediate and poor LAPW basis sets.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: The equilibrium bond length, harmonic frequency, first and second order anharmonicity constants, rotational and centrifugal distortion constants, as well as the rotation-vibrational and centrifugal coupling constants for the ground X(2) Sigma(+) and first excited A(2) Pi states of the SiN radical have been calculated at the complete active space SCF (CASSCF), multireference CI (MRCI) and coupled cluster (CCSD(T)) levels using Dunning's correlation-consistent basis sets. The excitation energy of the A(2) Pi State has also been computed at these theoretical levels. Dipole moments of SiN in the X(2) Sigma(+) and A(2) Pi states are given. Our study shows that core correlation must be considered in order to obtain satisfactory accuracy for the spectroscopic constants.

Abstract: Multi-wall carbon nanotubes have been produced by the catalytic decomposition of acetylene. Go-Mo, Co-V and Co-Fe mixtures supported either on zeolite or corundum alumina were used as catalysts. When Fe or V is added to Co, the carbon deposit increases. The nanotubes were characterized by both low and high resolution TEM. From histograms representing the outer diameter distributions, it is clear that the outer diameter of the nanotubes can be controlled by choosing the appropriate catalyst. (C) 2000 Elsevier Science B.V. All rights reserved.

Abstract: A molecular dynamics (MD) study of MgxAlyOz thin films grown by magnetron sputtering is presented using an ionic model and comparing two potential sets with formal and partial charges. The applicability of the model and the reliability of the potential sets for the simulation of thin film growth are discussed. The formal charge potential set was found to reproduce the thin film structure in close agreement with the structure of the experimentally grown thin films. Graphical abstract A molecular dynamics study of growth of MgxAlyOz thin films is presented using an ionic model and comparing two potential sets with formal and partial charges. The simulation results with the formal charge potential set showed a transition in the film from a crystalline to an amorphous structure, when the Mg metal content decreases below 50% in very close agreement with the structure of the experimentally deposited films.

Abstract: Despite significant progress in single walled carbon nanotube (SWCNT) production by plasma enhanced chemical vapor deposition (PECVD), the growth mechanism in this method is not clearly understood. We employ reactive molecular dynamics simulations to investigate how plasma-based deposition allows growth at low temperature. We first investigate the SWCNT growth mechanism at low and high temperatures under conditions similar to thermal CVD and PECVD. We then show how ion bombardment during the nucleation stage increases the carbon solubility in the catalyst at low temperature. Finally, we demonstrate how moderate energy ions sputter amorphous carbon allowing for SWCNT growth at 500 K. (C) 2013 Elsevier B. V. All rights reserved.

Abstract: Pristine and oxygen plasma functionalised carbon nanotubes (CNTs) were studied after the evaporation of Pt and Pd atoms. High resolution transmission electron microscopy shows the formation of metal nanoparticles at the CNT surface. Oxygen functional groups grafted by the plasma functionalization act as nucleation sites for metal nanoparticles. Analysis of the C1s core level spectra reveals that there is no covalent bonding between the Pt or Pd atoms and the CNT surface. Unlike other transition metals such as titanium and copper, neither Pd nor Pt show strong oxygen interaction or surface oxygen scavenging behaviour.