Search

Simulations of water near extended hydrophobic spherical solutes have revealed the presence of a region of depleted density and accompanying enhanced density fluctuations. The physical origin of both phenomena has remained somewhat obscure. We investigate these effects employing a mesoscopic binding potential analysis, classical density functional theory (DFT) calculations for a simple Lennard-Jones solvent, and Grand Canonical Monte Carlo (GCMC) simulations of a monatomic water (mw) model. We argue that the density depletion and enhanced fluctuations are near-critical phenomena. Specifically, we show that they can be viewed as remnants of the critical drying surface phase transition that occurs at bulk liquid–vapor coexistence in the macroscopic planar limit, i.e., as the solute radius Rs → ∞. Focusing on the radial density profile ρ(r) and a sensitive spatial measure of fluctuations, the local compressibility profile χ(r), our binding potential analysis provides explicit predictions for the manner in which the key features of ρ(r) and χ(r) scale with Rs, the strength of solute–water attraction ɛsf, and the deviation from liquid–vapor coexistence of the chemical potential, δμ. These scaling predictions are confirmed by our DFT calculations and GCMC simulations. As such, our theory provides a firm basis for understanding the physics of hydrophobic solvation. [ABSTRACT FROM AUTHOR]

Titanium dioxide has been extensively studied in the rutile or anatase phase, while its high-pressure phases are less well-understood, despite that many are thought to have interesting optical, mechanical, and electrochemical properties. First-principles methods, such as density functional theory (DFT), are often used to compute the enthalpies of TiO2 phases at 0 K, but they are expensive and, thus, impractical for long time scale and large system-size simulations at finite temperatures. On the other hand, cheap empirical potentials fail to capture the relative stabilities of various polymorphs. To model the thermodynamic behaviors of ambient and high-pressure phases of TiO2, we design an empirical model as a baseline and then train a machine learning potential based on the difference between the DFT data and the empirical model. This so-called Δ-learning potential contains long-range electrostatic interactions and predicts the 0 K enthalpies of stable TiO2 phases that are in good agreement with DFT. We construct a pressure–temperature phase diagram of TiO2 in the range 0 < P < 70 GPa and 100 < T < 1500 K. We then simulate dynamic phase transition processes by compressing anatase at different temperatures. At 300 K, we predominantly observe an anatase-to-baddeleyite transformation at about 20 GPa via a martensitic two-step mechanism with a highly ordered and collective atomic motion. At 2000 K, anatase can transform into cotunnite around 45–55 GPa in a thermally activated and probabilistic manner, accompanied by diffusive movement of oxygen atoms. The pressures computed for these transitions show good agreement with experiments. Our results shed light on how to synthesize and stabilize high-pressure TiO2 phases, and our method is generally applicable to other functional materials with multiple polymorphs. [ABSTRACT FROM AUTHOR]

If, in a hard sphere fluid, a single (test) particle is fixed, the other particles display a density profile that possesses long-ranged oscillations. Surprisingly, one can show via classical density functional theory that it takes a simple, purely repulsive (external) potential with a finite range in addition to the fixed hard sphere that forces these oscillations to vanish completely. This can give rise to interesting phenomena; however, it gained little attention in the past. In this work, we use the potential in question as an inter-component interaction in a binary hard-sphere mixture, where it is shown that the effective interaction induced by one component resembles qualitatively the well-known Asakura–Oosawa–Vrij potential and can lead to a liquid–gas phase transition in the other component. [ABSTRACT FROM AUTHOR]

Coarse-grained models of lyotropic solutions of semiflexible polymers are studied by both molecular dynamics simulations and density functional theory calculations, using an implicit solvent bead-spring model with a bond-angle potential. We systematically vary the monomer density, persistence length, and contour length over a wide range and explore the full range from the isotropic-nematic transition to the nematic-smectic transition. In the nematic regime, we span the entire regime from rigid-rod like polymers to thin wormlike chains, confined in effective straight tubes caused by the collective nematic effective ordering field. We show that the distribution of bond angles relative to the director is well described by a Gaussian, irrespective of whether the chains are rod-like or rather flexible. However, the related concept of "deflection length" is shown to make sense only in the latter case for rather dilute solutions since otherwise the deflection length is of the order of about two bond lengths only. When the solution is semi-dilute, a substantial renormalization of the persistence length occurs, while this effect is absent in the isotropic phase even at rather high monomer densities. The effective radii of the "tubes" confining the chains in the related description of orientational ordering are significantly larger than the distances between neighboring chains, providing evidence for a pronounced collective character of orientational fluctuations. Hairpins can be identified close to the isotropic-nematic transition, and their probability of occurrence agrees qualitatively with the Vroege-Odijk theory. The corresponding theoretical predictions for the elastic constants, however, are not in good agreement with the simulations. We attribute the shortcomings of the theories to their neglect of the coupling between local density and orientational fluctuations. Finally, we detected for this model a transition to a smectic phase for reduced monomer densities near 0.7. [ABSTRACT FROM AUTHOR]

We report here a new force field for water based solely on quantum mechanics (QM) calculations with no empirical data. The QM was at a high level, coupled cluster single double triple, for all orientations and distances for water dimer plus X3LYP density functional theory (DFT) on 19 larger water clusters. In addition, we included charge and polarization based on the polarizable charge equilibration method and nonbond interactions from DFT-D3 calculations on the H2 and O2 crystal. This model, denoted as RexPoN, provides quite excellent agreement with experimental (expr) data for the solid and liquid phase of water: Tmelt= 273.3 K (expr = 273.15 K) and properties at 298 K: ΔHvap = 10.36 kcal/mol (expr = 10.52), density = 0.9965 gr/cm3 (expr = 0.9965), entropy = 68.4 (J/mol)/K (expr = 69.9), dielectric constant = 76.1 (expr = 78.4), and ln Ds (self-diffusion coef) = −10.08 (expr = −11.24). Such an accurate force field for water will, we believe, be useful for full solvent calculations of electrocatalysis, where we can restrict QM water to just the first one or two layers involving reactions, using RexPoN to provide the polarization for a more distant solvent. Also, RexPoN may provide a better description of the solvent for proteins, DNA, polymers, and inorganic systems for applications to biomolecular, pharma, electrocatalysis (fuel cells and water splitting), and batteries where interaction with explicit water molecules plays a significant role. [ABSTRACT FROM AUTHOR]

In this study, we perform a systematic search to find the possible lowest energy structure of silicon nanoclusters Sin (n = 8-80) by means of an evolutionary algorithm. The fitness function for this search is the total energy of density functional tight binding (DFTB). To be on firm ground, we take several low energy structures of DFTB and perform further geometrical optimization by density functional theory (DFT). Then we choose structures with the lowest DFT total energy and compare them with the reported lowest energy structures in the literature. In our search, we found several lowest energy structures that were previously unreported. We further observe a geometrical transition at n = 27 from elongated to globular structures. In addition, the optical gap of the lowest energy structures is investigated by time-dependent DFTB (TD-DFTB) and time-dependent DFT (TD-DFT). The results show the same trend in TD-DFTB and TD-DFT for the optical gap. We also find a sudden drop in the optical gap at n = 27, precisely where the geometrical transition occurs. [ABSTRACT FROM AUTHOR]

We propose a new formulation of the correlation energy functional derived from the transcorrelated method in use in density functional theory (TC-DFT). An effective Hamiltonian, HTC, is introduced by a similarity transformation of a many-body Hamiltonian, H, with respect to a complex function F: HTC=1/F HF. It is proved that an expectation value of HTC for a normalized single Slater determinant, Dn, corresponds to the total energy: E[n]=⟨Ψn|H|Ψn⟩/⟨Ψn|Ψn⟩=⟨Dn|HTC|Dn⟩ under the two assumptions: (1) The electron density n(r) associated with a trial wave function Ψn = DnF is v-representable and (2) Ψn and Dn give rise to the same electron density n(r). This formulation, therefore, provides an alternative expression of the total energy that is useful for the development of novel correlation energy functionals. By substituting a specific function for F, we successfully derived a model correlation energy functional, which resembles the functional form of the screened exchange method. The proposed functional, named the extended screened exchange (ESX) functional, is described within two-body integrals and is parametrized for a numerically exact correlation energy of the homogeneous electron gas. The ESX functional does not contain any ingredients of (semi-)local functionals and thus is totally free from self-interactions. The computational cost for solving the self-consistent-field equation is comparable to that of the Hartree-Fock method. We apply the ESX functional to electronic structure calculations for a solid silicon, H- ion, and small atoms. The results demonstrate that the TC-DFT formulation is promising for the systematic improvement of the correlation energy functional. [ABSTRACT FROM AUTHOR]

Enhancing Raman signatures of molecules by self-assembled metal nanoparticles, nanolithography patterning, or by designing plasmonic nanostructures is widely used for detection of low abundance biological systems. Self-assembled peptide nanostructures provide a natural template for tethering Au and Ag nanoparticles due to its fractal surface. Here, we show the use of L,L-diphenylalanine micro-nanostructures (FF-MNSs) for the organization of Ag and Au nanoparticles (Nps) and its potential as surface-enhanced Raman scattering (SERS)-active substrates. The FF-MNSs undergo an irreversible phase transition from hexagonally packed (hex) micro-nanotubes to an orthorhombic (ort) structure at ~150 °C. The metal Nps form chains on hex FF-MNSs as inferred from transmission electron microscopy images and a uniform non-aggregated distribution in the ort phase. The high luminescence from the ort FF-MNS phase precludes SERS measurements with AgNps. The calculated Raman spectra using density-functional theory shows a higher intensity from rhodamine 6G (R6G) molecule in the presence of an Ag atom bound to ort FF compared with hex FF. The SERS spectra obtained from R6G bound to FF-MNSs with AuNps clearly show a higher enhancement for the ort phase compared with hex FF, corroborating our theoretical calculations. Our results indicate that FFMNSs both in the hex and ort phases can be used as substrates for the SERS analysis with different metal nanoparticles, opening up a novel class of optically active bio-based substrates. [ABSTRACT FROM AUTHOR]

We report a theoretical and simulation study of the drying and wetting phase transitions of a truncated Lennard-Jones fluid at a flat structureless wall. Binding potential calculations predict that the nature of these transitions depends on whether the wall-fluid attraction has a long ranged (LR) power law decay or is instead truncated, rendering it short ranged (SR). Using grand canonical Monte Carlo simulation and classical density functional theory, we examine both cases in detail. We find that for the LR case wetting is first order, while drying is continuous (critical) and occurs exactly at zero attractive wall strength, i.e., in the limit of a hardwall. In the SR case, drying is also critical but the order of the wetting transition depends on the truncation range of the wall-fluid potential. We characterize the approach to critical drying and wetting in terms of the density and local compressibility profiles and via the finitesize scaling properties of the probability distribution of the overall density. For the LR case, where the drying point is known exactly, this analysis allows us to estimate the exponent ν║, which controls the parallel correlation length, i.e., the extent of vapor bubbles at the wall. Surprisingly, the value we obtain is over twice that predicted by mean field and renormalization group calculations, despite the fact that our three dimensional system is at the upper critical dimension where mean field theory for critical exponents is expected to hold. Possible reasons for this discrepancy are discussed in the light of fresh insights into the nature of near critical finite-size effects. [ABSTRACT FROM AUTHOR]

Classical density functional theory is applied to investigate the validity of a phenomenological force-balance description of the stability of the Cassie state of liquids on substrates with nanoscale corrugation. A bulk free-energy functional of third order in local density is combined with a square-gradient term, describing the liquid-vapor interface. The bulk free energy is parameterized to reproduce the liquid density and the compressibility of water. The square-gradient term is adjusted to model the width of the water-vapor interface. The substrate is modeled by an external potential, based upon the Lennard-Jones interactions. The three-dimensional calculation focuses on substrates patterned with nanostripes and square-shaped nanopillars. Using both the force-balance relation and density-functional theory, we locate the Cassie-to-Wenzel transition as a function of the corrugation parameters. We demonstrate that the force-balance relation gives a qualitatively reasonable description of the transition even on the nanoscale. The force balance utilizes an effective contact angle between the fluid and the vertical wall of the corrugation to parameterize the impalement pressure. This effective angle is found to have values smaller than the Young contact angle. This observation corresponds to an impalement pressure that is smaller than the value predicted by macroscopic theory. Therefore, this effective angle embodies effects specific to nanoscopically corrugated surfaces, including the finite range of the liquid-solid potential (which has both repulsive and attractive parts), line tension, and the finite interface thickness. Consistently with this picture, both patterns (stripes and pillars) yield the same effective contact angles for large periods of corrugation. [ABSTRACT FROM AUTHOR]

Recent technical advances in dealing with finite-size errors make quantum Monte Carlo methods quite appealing for treating extended systems in electronic structure calculations, especially when commonly used density functional theory (DFT) methods might not be satisfactory. We present a theoretical study of martensitic phase transition energetics of a two-dimensional phosphorene by employing diffusion Monte Carlo (DMC) approach. The DMC calculation supports DFT prediction of having a rather diffusive barrier that is characterized by having two transition states, in addition to confirming that the so-called black and blue phases of phosphorene are essentially degenerate. At the same time, the DFT calculations do not provide the quantitative accuracy in describing the energy changes for the martensitic phase transition even when hybrid exchange-correlation functional is employed. We also discuss how mechanical strain influences the stabilities of the two phases of phosphorene. [ABSTRACT FROM AUTHOR]

A free energy functional for a crystal which contains both the symmetry-conserved and symmetrybroken parts of the direct pair correlation function has been used to investigate the fluid-solid transition in systems interacting via purely repulsive Weeks-Chandler-Anderson Lennard-Jones potential and the full Lennard-Jones potential. The results found for freezing parameters for the fluid-face centred cubic crystal transition are in very good agreement with simulation results. It is shown that although the contribution made by the symmetry broken part to the grand thermodynamic potential at the freezing point is small compared to that of the symmetry conserving part, its role is crucial in stabilizing the crystalline structure and on values of the freezing parameters. [ABSTRACT FROM AUTHOR]

In this article, we develop systematically second random phase approximations (RPA) and Tamm-Dancoff approximations (TDA) of particle-hole and particle-particle channels for calculating molecular excitation energies. The second particle-hole RPA/TDA can capture double excitations missed by the particle-hole RPA/TDA and time-dependent density-functional theory (TDDFT), while the second particle-particle RPA/TDA recovers non-highest-occupied-molecular-orbital excitations missed by the particle-particle RPA/TDA. With proper orbital restrictions, these restricted second RPAs and TDAs have a formal scaling of only O(N4). The restricted versions of second RPAs and TDAs are tested with various small molecules to show some positive results. Data suggest that the restricted second particle-hole TDA (r2ph-TDA) has the best overall performance with a correlation coefficient similar to TDDFT, but with a larger negative bias. The negative bias of the r2ph-TDA may be induced by the unaccounted ground state correlation energy to be investigated further. Overall, the r2ph-TDA is recommended to study systems with both single and some low-lying double excitations with a moderate accuracy. Some expressions on excited state property evaluations, such as {S²} are also developed and tested [ABSTRACT FROM AUTHOR]

A deuterated sample of CO2 structure I (sI) clathrate hydrate (CO2 • 8.3 D2O) has been formed and neutron diffraction experiments up to 1.0 GPa at 240 K were performed. The sI CO2 hydrate transformed at 0.7 GPa into the high pressure phase that had been observed previously by Hirai et al. [J. Phys. Chem. 133, 124511 (2010)] and Bollengier et al. [Geochim. Cosmochim. Acta 119, 322 (2013)], but which had not been structurally identified. The current neutron diffraction data were successfully fitted to a filled ice structure with CO2 molecules filling the water channels. This CO2+water system has also been investigated using classical molecular dynamics and density func-tional ab initio methods to provide additional characterization of the high pressure structure. Both models indicate the water network adapts a MH-III "like" filled ice structure with considerable dis-order of the orientations of the CO2 molecule. Furthermore, the disorder appears to be a direct result of the level of proton disorder in the water network. In contrast to the conclusions of Bollengier et al., our neutron diffraction data show that the filled ice phase can be recovered to ambient pres-sure (0.1 MPa) at 96 K, and recrystallization to sI hydrate occurs upon subsequent heating to 150 K, possibly by first forming low density amorphous ice. Unlike other clathrate hydrate systems, which transform from the sI or sII structure to the hexagonal structure (sH) then to the filled ice structure, CO2 hydrate transforms directly from the sI form to the filled ice structure. [ABSTRACT FROM AUTHOR]

Soft porous crystals present a challenge to molecular dynamics simulations with flexible size and shape of the simulation cell (i.e., in the NPT ensemble), since their framework responds very sensitively to small external stimuli. Hence, all interactions have to be described very accurately in order to obtain correct equilibrium structures. Here, we report a methodological study on the nanoporous metal-organic framework MIL-53(Ga), which undergoes a large-amplitude transition between a narrow- and a large-pore phase upon a change in temperature. Since this system has not been investigated by density functional theory (DFT)-based NPT simulations so far, we carefully check the convergence of the stress tensor with respect to computational parameters. Furthermore, we demonstrate the importance of dispersion interactions and test two different ways of incorporating them into the DFT framework. As a result, we propose two computational schemes which describe accurately the narrow- and the large-pore phase of the material, respectively. These schemes can be used in future work on the delicate interplay between adsorption in the nanopores and structural flexibility of the host material. [ABSTRACT FROM AUTHOR]