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: Dielectric properties of nano-confined water are important in several areas of science, i.e., it is relevant in the dielectric double layer that exists in practically all heterogeneous fluid-based systems. Molecular dynamics simulations are used to predict the in-plane dielectric properties of confined water in planar channels of width ranging from sub-nanometer to bulk. Because of suppressed rotational degrees of freedom near the confining walls, the dipole of the water molecules tends to be aligned parallel to the walls, which results in a strongly enhanced in-plane dielectric constant (epsilon (parallel to)) reaching values of about 120 for channels with height 8 angstrom < h < 10 angstrom. With the increase in the width of the channel, we predict that epsilon (parallel to) decreases nonlinearly and reaches the bulk value for h > 70 angstrom. A stratified continuum model is proposed that reproduces the h > 10 angstrom dependence of epsilon (parallel to). For sub-nanometer height channels, abnormal behavior of epsilon (parallel to) is found with two orders of magnitude reduction of epsilon (parallel to) around h similar to 7.5 angstrom, which is attributed to the formation of a particular ice phase that exhibits long-time (similar to mu s) stable ferroelectricity. This is of particular importance for the understanding of the influence of confined water on the functioning of biological systems.

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: Using first principles calculations, we study the structural and phononic properties of the recently synthesized Janus type single layers of molybdenum dichalcogenides. The Janus MoSSe single layer possesses 2H crystal structure with two different chalcogenide sides that lead to out-of-plane anisotropy. By virtue of the asymmetric structure of the ultra-thin Janus type crystal, we induced the out-of-plane anisotropy to show the distinctive vertical pressure effect on the vibrational properties of the Janus material. It is proposed that for the corresponding Raman active optical mode of the Janus structure, the phase modulation and the magnitude ratio of the strained atom and its first neighbor atom adjust the distinctive change in the eigen-frequencies and Raman activity. Moreover, a strong variation in the Raman activity of the Janus structure is obtained under bivertical and univertical strains. Not only eigen-frequency shifts but also Raman activities of the optical modes of the Janus structure exhibit distinguishable features. This study reveals that the vertical anisotropic feature of the Janus structure under Raman measurement allows us to distinguish which side of the Janus crystal interacts with the externals (substrate, functional adlayers, or dopants). Published by AIP Publishing.

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: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: The electronic properties of graphene nano-flakes (GNFs) with different edge passivation are investigated by using density functional theory. Passivation with F and H atoms is considered: C-Nc X-Nx (X = F or H). We studied GNFs with 10 < N-c < 56 and limit ourselves to the lowest energy configurations. We found that: (i) the energy difference Delta between the highest occupied molecular orbital and the lowest unoccupied molecular orbital decreases with N-c, (ii) topological defects (pentagon and heptagon) break the symmetry of the GNFs and enhance the electric polarization, (iii) the mutual interaction of bilayer GNFs can be understood by dipole-dipole interaction which were found sensitive to the relative orientation of the GNFs, (iv) the permanent dipoles depend on the edge terminated atom, while the energy gap is independent of it, and (v) the presence of heptagon and pentagon defects in the GNFs results in the largest difference between the energy of the spin-up and spin-down electrons which is larger for the H-passivated GNFs as compared to F-passivated GNFs. Our study shows clearly the effect of geometry, size, termination, and bilayer on the electronic properties of small GNFs. This study reveals important features of graphene nano-flakes which can be detected using Raman spectroscopy. (C) 2014 AIP Publishing LLC.

Abstract: Photocatalytic soot oxidation is studied on P25 TiO2 as an important model reaction for self-cleaning processes by means of electron paramagnetic resonance (EPR) and Fourier transform infrared (FTIR) spectroscopy. Contacting of carbon black with P25 leads on the one hand to a reduction of the local dioxygen concentration in the powder. On the other hand, the weakly adsorbed radicals on the carbon particles are likely to act as alternative traps for the photogenerated conduction-band electrons. We find furthermore that the presence of dioxygen and oxygen-related radicals is vital for the photocatalytic soot degradation. The complete oxidation of soot to CO2 is evidenced by in situ FTIR spectroscopy, no intermediate CO is detected during the photocatalytic process.

Abstract: The present work is devoted to the preparation of Ag/TiO2 nanosystems by an original synthetic strategy, based on the radio-frequency (RF) sputtering of silver particles on titania-based xerogels prepared by the sol-gel (SG) route. This approach takes advantage of the synergy between the microporous xerogel structure and the infiltration power characterizing RF-sputtering, whose combination enables the obtainment of a tailored dispersion of Ag-containing particles into the titania matrix. In addition, the systems chemico-physical features can be tuned further through proper ex situ thermal treatments in air at 400 and 600 °C. The synthesized composites are extensively characterized by the joint use of complementary techniques, that is, X-ray photoelectron and X-ray excited Auger electron spectroscopies (XPS, XE-AES), secondary ion mass spectrometry (SIMS), glancing incidence X-ray diffraction (GIXRD), field emission scanning electron microscopy (FE-SEM), transmission electron microscopy (TEM), electron diffraction (ED), high-angle annular dark field scanning TEM (HAADF-STEM), energy-filtered TEM (EF-TEM) and optical absorption spectroscopy. Finally, the photocatalytic performances of selected samples in the decomposition of the azo-dye Plasmocorinth B are preliminarily investigated. The obtained results highlight the possibility of tailoring the system characteristics over a broad range, directly influencing their eventual functional properties.

Abstract: We study the structure and phase behavior of fullerene-cubane C60·C8H8 by Monte Carlo simulation. Using a simple potential model capturing the icosahedral and cubic symmetries of its molecular constituents, we reproduce the experimentally observed phase transition from a cubic to an orthorhombic crystal lattice and the accompanying rotational freezing of the C60 molecules. We elaborate a scheme to identify the low-temperature orientations of individual molecules and to detect a pattern of orientational ordering similar to the arrangement of C60 molecules in solid C60. Our configuration of orientations supports a doubled periodicity along one of the crystal axes.

Abstract: Single-chain elasticity of polyethylene at theta point up to 90% of stretching with respect to its contour length is computed by Monte Carlo simulation of an atomistic model in continuous space. The elasticity law together with the free-energy and the internal energy variations with stretching are found to be very well represented by the wormlike chain model up to 65% of the chain elongation, provided the persistence length is treated as a temperature-dependent parameter. Beyond this value of elongation simple ideal chain models are not able to describe the Monte Carlo data in a thermodynamic consistent way. This study reinforces the use of the wormlike chain model to interpret experimental data on the elasticity of synthetic polymers in the finite extensibility regime, provided the chain is not yet in its fully stretched regime. Specific solvent effects on the elasticity law and the partition between energetic and entropic contributions to single chain elasticity are investigated. (C) 2002 American Institute of Physics.

Abstract: A model of the intermolecular potential in solid C-60, which is based on Born-Mayer repulsions, van der Waals attractions, and electrostatic multipoles, is presented. The potential is expanded in terms of multipolar rotator functions. The orientation-orientation interaction and the crystal field are calculated. The orientational phase transition to the Pa3 phase is studied with the methods of statistical mechanics. The discontinuity of the order parameter at the transition and the temperature evolution of the order parameter are calculated. The lattice contraction at the phase transition is evaluated. The influence of the lattice contraction on the crystal field and on the orientational order is studied.