Accurate first-principles structures and energies of diversely bonded systems from an efficient density functional

Regularized SCAN functional

We propose modifications to the functional form of the Strongly Constrained and Appropriately Normed (SCAN) density functional to eliminate numerical instabilities. This is necessary to allow reliable, automatic generation of pseudopotentials (including projector augmented-wave potentials). The regularized SCAN is designed to match the original form very closely, and we show that its performance remains comparable.

Room temperature electrofreezing of water yields a missing dense ice phase in the phase diagram

Water can freeze into diverse ice polymorphs depending on the external conditions such as temperature (T) and pressure (P). Herein, molecular dynamics simulations show evidence of a high-density orthorhombic phase, termed ice χ, forming spontaneously from liquid water at room temperature under high-pressure and high external electric field. Using free-energy computations based on the Einstein molecule approach, we show that ice χ is an additional phase introduced to the state-of-the-art T–P phase diagram. The χ phase is the most stable structure in the high-pressure/low-temperature region, located between ice II and ice VI, and next to ice V exhibiting two triple points at 6.06 kbar/131.23 K and 9.45 kbar/144.24 K, respectively. A possible explanation for the missing ice phase in the T–P phase diagram is that ice χ is a rare polarized ferroelectric phase, whose nucleation/growth occurs only under very high electric fields.

Water can crystallize in different ice polymorphs according to temperature and pressure conditions. Here the authors predict by molecular dynamics simulations a new ice phase spontaneously forming at room temperature under high pressure and high electric field.

Atomistic Insight into Glide-Driven Phase Transformations in Layered Oxides for Sodium-Ion Batteries: A Case Study on Na xVO2

Among high-capacity, low-cost cathode contenders for Na-ion batteries, layered transition-metal oxides are particularly promising materials. Yet there is a strong need to improve their long-term stability and capacity retention due to unwanted phase transitions occurring during sodium insertion and extraction cycles. Here, using density functional theory calculations and thoroughly examining P ↔ O transitions in Na _xVO₂ as a case study, we provide atomic-level understanding into the glide-driven processes in these compounds. In particular, we rationalize the commonly observed effect of high Na contents on preventing gliding processes at low states of charge. And we identify the beneficial impact of Na ions at face-shared prismatic sites in P-type phases to specifically hinder P → O transitions. We put forward that the gained new insights can help to focus and guide the current research efforts to design layered transition-metal oxides with enhanced long-term stability for practical applications.

Pancakes under Pressure: A Case Study on Isostructural Dithia- and Diselenadiazolyl Radical Dimers

The isostructural dimers of the 1,4-phenylene-bridged bis-1,2,3,5-dithia- and bis-1,2,3,5-diselenadiazolyl diradicals 1,4-S/Se are small band gap semiconductors. The response of their molecular and solid state electronic structures to pressure has been explored over the range 0-10 GPa. The crystal structures, which consist of cofacially aligned (pancake) π-dimers packed into herringbone arrays, experience a continuous, near-isotropic compression. While the intramolecular covalent E-E (E = S/Se) bonds remain relatively unchanged with pressurization, the intradimer E···E separations are significantly shortened. Molecular and band electronic structure calculations using density functional theory methods indicate that compression of the π-dimers leads to a widening of the gap Δ E between the highest occupied and lowest unoccupied molecular orbitals of the dimer, an effect that offsets the expected decrease in the valence-to-conduction band gap E_g occasioned by pressure-induced spreading of the valence and conduction bands. Consistent with the predicted consequences of this competition between intra- and interdimer interactions, variable temperature high pressure conductivity measurements reveal at best an order-of-magnitude increase in conductivity with pressure for the two compounds over the pressure range 0-10 GPa. While a small reduction in the thermal activation energy E_act with increasing pressure is observed, extrapolation of the rate of decrease suggests a projected onset of metallization ( E_act ≈ 0) in excess of 20 GPa.

In-depth structural analysis of glycans in the gas phase

Although there have been substantial improvements in glycan analysis over the past decade, the lack of both high-resolution and high-throughput methods hampers progress in glycomics. This perspective article highlights the current developments of liquid chromatography, mass spectrometry, ion-mobility spectrometry and cryogenic IR spectroscopy for glycan analysis and gives a critical insight to their individual strengths and limitations. Moreover, we discuss a novel concept in which ion mobility-mass spectrometry and cryogenic IR spectroscopy is combined in a single instrument such that datasets consisting of m/z, collision cross sections and IR fingerprints can be obtained. This multidimensional data will then be compared to a comprehensive reference library of intact glycans and their fragments to accurately identify unknown glycans on a high-throughput scale with minimal sample requirements. Due to the complementarity of the obtained information, this novel approach is highly diagnostic and also suitable for the identification of larger glycans; however, the workflow and instrumentation is straightforward enough to be implemented into a user-friendly setup.

Pressure induced excellent thermoelectric behavior in skutterudites CoSb3 and IrSb3

Utilizing the first-principle calculations combined with Boltzmann transport equation (BTE) and semiclassical analysis, we have systematically investigated the electronic structure, lattice thermal conductivity κL, Seebeck coefficient S, and the dimensionless figure of merit zT as a function of hydrostatic pressure P in crystalline skutterudites CoSb3 and IrSb3. Interestingly, as the pressure increases, the band gap and κL show an approximate parabolic trend, which results in extraordinarily high S and excellent thermoelectric properties, and zT even exceeds 1.4(1.09) in IrSb3(CoSb3) at 54(58) GPa. This anomalous behavior arises from the electron distribution and intrinsic scattering processes. Further analyses indicate that (i) nonbonding electron pairs of Sb atoms are gradually transferred to the region between Co(Ir) and Sb atoms as the pressure increases, which leads to the formation of a partial metallic bond and thus the band gap first expands and then shrinks; (ii) the change of the strength of the anharmonic phonon scattering process results in the variation of κL. As a result, these behaviors cause excellent thermoelectric properties. Our results provide insight into the thermal transport properties of skutterudites, meanwhile, forecast potential high pressure applications for thermoelectric materials.

Silicon Liquid Structure and Crystal Nucleation from Ab Initio Deep Metadynamics

Studying the crystallization process of silicon is a challenging task since empirical potentials are not able to reproduce well the properties of both a semiconducting solid and metallic liquid. On the other hand, nucleation is a rare event that occurs in much longer timescales than those achievable by ab initio molecular dynamics. To address this problem, we train a deep neural network potential based on a set of data generated by metadynamics simulations using a classical potential. We show how this is an effective way to collect all the relevant data for the process of interest. In order to efficiently drive the crystallization process, we introduce a new collective variable based on the Debye structure factor. We are able to encode the long-range order information in a local variable which is better suited to describe the nucleation dynamics. The reference energies are then calculated using the strongly constrained and appropriately normed (SCAN) exchange-correlation functional, which is able to get a better description of the bonding complexity of the Si phase diagram. Finally, we recover the free energy surface with a density functional theory accuracy, and we compute the thermodynamics properties near the melting point, obtaining a good agreement with experimental data. In addition, we study the early stages of the crystallization process, unveiling features of the nucleation mechanism.

Shape-controllable synthesis of GdVO4 photocatalysts and their tunable properties in photocatalytic hydrogen generation

Novel visible light responsive materials for water splitting are essential for the efficient conversion of solar energy into hydrogen bond energy. Among other semiconductors, gadolinium orthovanadate has appropriate conduction and valence band edges positioned to split water molecules and a narrow band gap that allows the use of visible light for hydrogen generation. Thus, we present here that hydrogen evolution under visible light (λ > 420 nm) could be accomplished using hierarchical 3D GdVO₄ particles, obtained by a simple, one pot hydrothermal synthesis. We found that applying various reaction components, such as EDTA-Na₂ and EDTA, and adjusting the pH of the solution allow one to tune the shape of GdVO₄ (such as short nanowires, long nanowires, short nanorods, long nanorods, nanoparticles and spheres - all having a tetragonal crystal structure) as well as optical and photocatalytic properties. The highest ability to photocatalytically split methanol solution into hydrogen under UV-Vis irradiation was detected for the long nanowire sample (42 μmol h^-1), having almost 11 times higher efficiency in comparison with the weakest sample - short nanowires. In addition, GdVO₄ spheres generated H₂ more than 2 times (5.75 μmol h^-1) in comparison with the short nanorod sample (2.5 μmol h^-1) under visible light excitation. Photostable in three-hour work cycles, long nanowires and spheres were even able to generate hydrogen from pure water, reaching values of 17 and 3 μmol under UV-Vis and Vis light, respectively.

Structure, Polarization, and Sum Frequency Generation Spectrum of Interfacial Water on Anatase TiO2

The photocatalytic activity of TiO₂ for water splitting has been known for decades, yet the adsorption structure and hydrogen bonding of water at the interface with TiO₂ have remained controversial. We investigate the prototypical aqueous interface with anatase TiO₂ (101) using ab initio molecular dynamics (AIMD) with the strongly constrained and appropriately normed (SCAN) density functional, recently shown to provide an excellent description of the properties of bulk liquid water. We find that water forms a stable bilayer of intact molecules with ice-like dynamics and enhanced dipole moment and polarizability on the anatase surface. The orientational order and H-bond environment of interfacial water are reflected in the computed sum frequency generation (SFG) spectrum, which agrees well with recent measurements in the OH stretching frequency range (3000-3600 cm^-1). Additional AIMD simulations for a model interface with 66% of dissociated water in the contact layer show that surface hydroxyls disrupt the order in the bilayer and lead to a much faster orientational dynamics of interfacial water. Nonetheless, the computed SFG spectrum for the hydroxylated surface also agrees with experiment, suggesting that SFG measurements in a wider frequency range would be necessary to unambiguously identify the character of interfacial water on anatase.

Structure stabilization effect of configuration entropy in cubic WN

The microscopic structure of cubic WN has been studied combining scanning transmission electron microscopy and first-principles calculations. Because of the contribution of configurational entropy, NaCl-type WN with disordered vacancies becomes more stable at high temperatures than NbO-type WN. Moreover, electron beam irradiation can induce an order-disorder transition in cubic WN. It is suggested that the ordered NbO-type WN can be obtained after annealing below the transition temperature. The results shed light on the stability of materials synthesized at high pressures and high temperatures.

Tuning the balance between dispersion and entropy to design temperature-responsive flexible metal-organic frameworks

Temperature-responsive flexibility in metal-organic frameworks (MOFs) appeals to the imagination. The ability to transform upon thermal stimuli while retaining a given crystalline topology is desired for specialized sensors and actuators. However, rational design of such shape-memory nanopores is hampered by a lack of knowledge on the nanoscopic interactions governing the observed behavior. Using the prototypical MIL-53(Al) as a starting point, we show that the phase transformation between a narrow-pore and large-pore phase is determined by a delicate balance between dispersion stabilization at low temperatures and entropic effects at higher ones. We present an accurate theoretical framework that allows designing breathing thermo-responsive MOFs, based on many-electron data for the dispersion interactions and density-functional theory entropy contributions. Within an isoreticular series of materials, MIL-53(Al), MIL-53(Al)-FA, DUT-4, DUT-5 and MIL-53(Ga), only MIL-53(Al) and MIL-53(Ga) are proven to switch phases within a realistic temperature range.

Rational design of metal organic frameworks (MOFs) with shape-memory nanopores is a formidable challenge. Here the authors use an accurate theoretical approach to design thermo-responsive MOFs based on a balance of van der Waals and entropy contributions.

Applicability of the Strongly Constrained and Appropriately Normed Density Functional to Transition-Metal Magnetism

We find that the recently developed self-consistent and appropriately normed meta-generalized gradient approximation, which has been found to provide highly accurate results for many materials, is, however, not able to describe the stability and properties of phases of Fe important for steel. This is due to an overestimated tendency toward magnetism and exaggeration of magnetic energies, which we also find in other transition metals.

Melting Si: Beyond Density Functional Theory

The melting point of silicon in the cubic diamond phase is calculated using the random phase approximation (RPA). The RPA includes exact exchange as well as an approximate treatment of local as well as nonlocal many body correlation effects of the electrons. We predict a melting temperature of about 1735 and 1640 K without and with core polarization effects, respectively. Both values are within 3% of the experimental melting temperature of 1687 K. In comparison, the commonly used gradient approximation to density functional theory predicts a melting point that is 200 K too low, and hybrid functionals overestimate the melting point by 150 K. We correlate the predicted melting point with the energy difference between cubic diamond and the beta-tin phase of silicon, establishing that this energy difference is an important benchmark for the development of approximate functionals. The current results demonstrate that the RPA can be used to predict accurate finite temperature properties and underlines the excellent predictive properties of the RPA for condensed matter.

Parity-breaking in single-element phases: ferroelectric-like elemental polar metals

Polar metals based on binary and ternary compounds have been demonstrated in literature. Here, we propose a design principle for ferroelectric-like elemental polar metals and relate it to real materials. The design principle is that, to be an elemental polar metal, atoms should occupy at least two inequivalent Wyckoff positions in a crystal with a polar space group, where inversion symmetry is spontaneously broken. According to this rule, we propose the first class of potential ferroelectric-like elemental polar metals in a distorted α-La-like structure with a polar space group P6₃ mc in which two inequivalent Wyckoff positions 2a (0, 0, z) and 2b (1/3, 2/3, z) are occupied by group-V elements (phosphorus, arsenic, antimony, and bismuth). Analyses based on first-principles calculations indicate that the dynamically stable polar phase results from a lone pair driven polar distortion of the nonploar phase in P6₃/mmc symmetry where two inequivalent Wyckoff positions 2a (0, 0, 0) and 2c (1/3, 2/3, 1/4) are occupied. This ferroelectric-like transition involves a transition from a metallic state to a semimetallic state. These predicted polar phases are metastable with respect to their corresponding ground phases. Moreover, ionic bonding characters are found due to the inequivalence in Wyckoff positions between group-V atoms. Our work opens a route to single-element parity-breaking phases.

A Systematic Protocol for Benchmarking Guest-Host Interactions by First-Principles Computations: Capturing CO2 in Clathrate Hydrates

Clathrate hydrates of CO₂ have been proposed as potential molecular materials in tackling important environmental problems related to greenhouse gases capture and storage. Despite the increasing interest in such hydrates and their technological applications, a molecular-level understanding of their formation and properties is still far from complete. Modeling interactions is a challenging and computationally demanding task, essential to reliably determine molecular properties. First-principles calculations for the CO₂ guest in all sI, sII, and sH clathrate cages were performed, and the nature of the guest-host interactions, dominated by both hydrogen-bond and van der Waals forces, was systematically investigated. Different families of density functionals, as well as pairwise CO₂ @H₂ O model potentials versus wavefunction-based quantum approaches were studied for CO₂ clathrate-like systems. Benchmark energies for new distance-dependent datasets, consisting of potential energy curves sampling representative configurations of the systems at the repulsive, near-equilibrium, and asymptotic/long-range regions of the full-dimensional surface, were generated, and a general protocol was proposed to assess the accuracy of such conventional and modern approaches at minimum and non-minimum orientations. Our results show that dispersion interactions are important in the guest-host stabilization energies of such clathrate cages, and the encapsulation of the CO₂ into guest-free clathrate cages is always energetically favorable. In addition, the orientation of CO₂ inside each cage was explored, and the ability of current promising approaches to accurately describe non-covalent CO₂ @H₂ O guest-host interactions in sI, sII, and sH clathrates was discussed, providing information for their applicability to future multiscale computer simulations.

Simple Modifications of the SCAN Meta-Generalized Gradient Approximation Functional

We analyzed various possibilities to improve upon the SCAN meta-generalized gradient approximation density functional obeying all known properties of the exact functional that can be satisfied at this level of approximation. We examined the necessity of locally satisfying a strongly tightened lower bound for the exchange energy density in single-orbital regions, the nature of the error cancellation between the exchange and correlation parts in two-electron regions, and the effect of the fourth-order term in the gradient expansion of the correlation energy density. We have concluded that the functional can be modified to separately reproduce the exchange and correlation energies of the helium atom by locally releasing the strongly tightened lower bound for the exchange energy density in single-orbital regions, but this leads to an unbalanced improvement in the single-orbital electron densities. Therefore, we decided to keep the F_X ≤ 1.174 exact condition for any single-orbital density, where F_X is the exchange enhancement factor. However, we observed a general improvement in the single-orbital electron densities by revising the correlation functional form to follow the second-order gradient expansion in a wider range. Our new revSCAN functional provides more-accurate atomization energies for the systems with multireference character, compared to the SCAN functional. The nonlocal VV10 dispersion-corrected revSCAN functional yields more-accurate noncovalent interaction energies than the VV10-corrected SCAN functional. Furthermore, its global hybrid version with 25% of exact exchange, called revSCAN0, generally performs better than the similar SCAN0 for reaction barrier heights. Here, we also analyzed the possibility of the construction of a local hybrid from the SCAN exchange and a specific locally bounded nonconventional exact exchange energy density. We predict compatibility problems since this nonconventional exact exchange energy density does not really obey the strongly tightened lower bound for the exchange energy density in single-orbital regions.

Perspective: Chemical Information Encoded in Electron Density

In this perspective, we review the chemical information encoded in electron density and other ingredients used in semilocal functionals. This information is usually looked at from the functional point of view: the exchange density or the enhancement factor are discussed in terms of the reduced density gradient. However, what parts of a molecule do these 3D functions represent? We look at these quantities in real space, aiming to understand the electronic structure information they encode and provide an insight from the quantum chemical topology (QCT). Generalized gradient approximations (GGAs) provide information about the presence of chemical interactions, whereas meta-GGAs can differentiate between the different bonding types. By merging these two techniques, we show new insight into the failures of semilocal functionals owing to three main errors: fractional charges, fractional spins, and non-covalent interactions. We build on simple models. We also analyze the delocalization error in hydrogen chains, showing the ability of QCT to reveal the delocalization error introduced by semilocal functionals. Then, we show how the analysis of localization can help understand the fractional spin error in alkali atoms, and how it can be used to correct it. Finally, we show that the poor description of GGAs of isodesmic reactions in alkanes is due to 1,3-interactions.

Nuclear quantum effect with pure anharmonicity and the anomalous thermal expansion of silicon

Despite the widespread use of silicon in modern technology, its peculiar thermal expansion is not well understood. Adapting harmonic phonons to the specific volume at temperature, the quasiharmonic approximation, has become accepted for simulating the thermal expansion, but has given ambiguous interpretations for microscopic mechanisms. To test atomistic mechanisms, we performed inelastic neutron scattering experiments from 100 K to 1,500 K on a single crystal of silicon to measure the changes in phonon frequencies. Our state-of-the-art ab initio calculations, which fully account for phonon anharmonicity and nuclear quantum effects, reproduced the measured shifts of individual phonons with temperature, whereas quasiharmonic shifts were mostly of the wrong sign. Surprisingly, the accepted quasiharmonic model was found to predict the thermal expansion owing to a large cancellation of contributions from individual phonons.

Expanding frontiers in materials chemistry and physics with multiple anions

During the last century, inorganic oxide compounds laid foundations for materials synthesis, characterization, and technology translation by adding new functions into devices previously dominated by main-group element semiconductor compounds. Today, compounds with multiple anions beyond the single-oxide ion, such as oxyhalides and oxyhydrides, offer a new materials platform from which superior functionality may arise. Here we review the recent progress, status, and future prospects and challenges facing the development and deployment of mixed-anion compounds, focusing mainly on oxide-derived materials. We devote attention to the crucial roles that multiple anions play during synthesis, characterization, and in the physical properties of these materials. We discuss the opportunities enabled by recent advances in synthetic approaches for design of both local and overall structure, state-of-the-art characterization techniques to distinguish unique structural and chemical states, and chemical/physical properties emerging from the synergy of multiple anions for catalysis, energy conversion, and electronic materials.

Origin of the size-dependence of the equilibrium van der Waals binding between nanostructures

Nanostructures can be bound together at equilibrium by the van der Waals (vdW) effect, a small but ubiquitous many-body attraction that presents challenges to density functional theory. How does the binding energy depend upon the size or number of atoms in one of a pair of identical nanostructures? To answer this question, we treat each nanostructure as a whole object, not as a collection of atoms. Our calculations start from an accurate static dipole polarizability for each considered nanostructure, and an accurate equilibrium center-to-center distance for the pair (the latter from experiment or from the vdW-DF-cx functional). We consider the competition in each term -C_2k/d^2k (k = 3, 4, 5) of the long-range vdW series for the interaction energy, between the size dependence of the vdW coefficient C_2k and that of the 2kth power of the center-to-center distance d. The damping of these vdW terms can be negligible, but in any case, it does not affect the size dependence for a given term in the absence of non-vdW binding. To our surprise, the vdW energy can be size-independent for quasi-spherical nanoclusters bound to one another by vdW interaction, even with strong nonadditivity of the vdW coefficient, as demonstrated for fullerenes. We also show that, for low-dimensional systems, the vdW interaction yields the strongest size-dependence, in stark contrast to that of fullerenes. We illustrate this with parallel planar polycyclic aromatic hydrocarbons. The size dependences of other morphologies or bonding types lie between, as shown by sodium clusters.

Assessing Intermolecular Interactions in Guest-Free Clathrate Hydrate Systems

Recently, empty hydrate structures sI, sII, sH, and others have been proposed as low-density ice structures by both experimental observations and computer simulations. Some of them have been synthesized in the laboratory, which motivates further investigations on the stability of such guest-free clathrate structures. Using semiempirical and ab initio-based water models, as well as dispersion-corrected density functional theory approaches, we predict their stability, including cooperative many-body effects, in comparison with reference data from converged wave function-based DF-MP2 electronic structure calculations. We show that large basis sets and counterpoise corrections are required to improve convergence in the interaction/binding energies for such systems. Therefore, extrapolation schemes based on triple/quadruple and quadruple/quintuple ζ quality basis sets are used to reach high accuracy. Eleven different water structures corresponding to dodecahedron, edge sharing, face sharing, and fused cubes, as a part of the WATER27 database, as well as cavities from the sI, sII, and sH clathrate hydrates formed by 20, 24, 28, and 36 water molecules, are employed, and new benchmark energies are reported. Using these benchmark sets of interaction energies, we assess the performance of both analytical models and direct DFT calculations for such clathrate-like systems. In particular, seven popular water models (TIP4P/ice, TIP4P/2005, q-TIP4P/F, TTM2-F, TTM3-F, TTM4-F, and MB-pol) available in the literature, and nine density functional approximations (3 meta-GGAs, 3 hybrids, and 3 range separated functionals) are used to investigate their accuracy. By including dispersion corrections, our results show that errors in the interaction energies are reduced for most of the DFT functionals. Despite the difficulties faced by current water models and DFT functionals to accurately describe the interactions in such water systems, we found some general trends that could serve to extend their applicability to larger systems.

Fast and accurate quantum Monte Carlo for molecular crystals

Computational approaches based on the fundamental laws of quantum mechanics are now integral to almost all materials design initiatives in academia and industry. If computational materials science is genuinely going to deliver on its promises, then an electronic structure method with consistently high accuracy is urgently needed. We show that, thanks to recent algorithmic advances and the strategy developed in our manuscript, quantum Monte Carlo yields extremely accurate predictions for the lattice energies of materials at a surprisingly modest computational cost. It is thus no longer a technique that requires a world-leading computational facility to obtain meaningful results. While we focus on molecular crystals, the significance of our findings extends to all classes of materials.

Computer simulation plays a central role in modern-day materials science. The utility of a given computational approach depends largely on the balance it provides between accuracy and computational cost. Molecular crystals are a class of materials of great technological importance which are challenging for even the most sophisticated ab initio electronic structure theories to accurately describe. This is partly because they are held together by a balance of weak intermolecular forces but also because the primitive cells of molecular crystals are often substantially larger than those of atomic solids. Here, we demonstrate that diffusion quantum Monte Carlo (DMC) delivers subchemical accuracy for a diverse set of molecular crystals at a surprisingly moderate computational cost. As such, we anticipate that DMC can play an important role in understanding and predicting the properties of a large number of molecular crystals, including those built from relatively large molecules which are far beyond reach of other high-accuracy methods.

Adsorption and decomposition of dimethyl methylphosphonate on size-selected (MoO3)3 clusters

The adsorption and decomposition of dimethyl methylphosphonate (DMMP), a chemical warfare agent (CWA) simulant, on size-selected molybdenum oxide trimer clusters, i.e. (MoO₃)₃, was studied both experimentally and theoretically.

Dramatic enhancement of spin–spin coupling and quenching of magnetic dimensionality in compressed silver difluoride

The high-pressure HP2 form of AgF₂ features Ag₂F₇³⁻ units that are theoretically predicted to host extremely strong antiferromagnetic interactions, surpassing those seen in copper(ii) oxides.

Low-temperature activation of methane on doped single atoms: descriptor and prediction

We predict that Pt and several other single atoms on rutile TiO₂(110) can chemisorb and activate methane at low temperatures.

Clathrate ice sL: a new crystalline phase of ice with ultralow density predicted by first-principles phase diagram computations

A new crystalline ice phase (clathrate ice sL) with ultralow density under negative pressure predicted by the first-principles phase diagram computations.

Chemically driven surface effects in polar intermetallic topological insulators A3Bi

Surface electronic spectra, surface and bulk properties as well as the underlying chemical bonding characteristics in topological insulators with complex bonding patterns are considered for the example of cubic, polar intermetallics KNa₂Bi, K₃Bi and Rb₃Bi (with the general formula A₃Bi, A – alkali metal).

Computational investigation of the α2β1 integrin–collagen triple helix complex interaction

In this paper, quantum biochemistry methods have been used to describe important protein–protein interactions for the complex integrin–collagen.

Electrochemical phase diagrams of Ni from ab initio simulations: role of exchange interactions on accuracy

The stabilities of Ni metal and its derived compounds, including oxides, hydroxides, and oxyhydroxides under electrochemical conditions, can be readily predicted from the Ni Pourbaix diagram, where the formation free energies of the involved species are utilized to construct the phase stability map with respect to electrode potential and pH. We calculate and analyze the crystal structures, electronic structures, and thermodynamic energies of Ni metal and its compounds using different exchange-correlation functionals to density-functional-theory (DFT), including the semilocal LDA and GGA density functionals, the nonlocal metaGGA, and the hybrid density functionals. Next, we simulate the corresponding Ni Pourbaix diagrams to compare systematically the performance of the functional to each other and to experimental observations. We show that the structures and energies obtained from experimental databases may not be sufficiently accurate to describe direct electrochemical observations, and we explain how the electronic exchange within the density functionals plays a key role in determining the accuracy of the DFT calculated electronic, thermodynamic, and electrochemical properties. We find that only the hybrid density functional produces reliable results owing to the fractional contribution of exact Fock exchange included therein. Last, based on our accurate Ni Pourbaix diagram, we construct band-gap and magnetic electrochemical maps which can facilitate more experimental measurements and property assessments under variable potential and pH in the future.

Properties of real metallic surfaces: Effects of density functional semilocality and van der Waals nonlocality

We have computed the surface energies, work functions, and interlayer surface relaxations of clean (111), (100), and (110) surfaces of Al, Cu, Ru, Rh, Pd, Ag, Pt, and Au. We interpret the surface energy from liquid metal measurements as the mean of the solid-state surface energies over these three lowest-index crystal faces. We compare experimental (and random phase approximation) reference values to those of a family of nonempirical semilocal density functionals, from the basic local density approximation (LDA) to our most advanced general purpose meta-generalized gradient approximation, strongly constrained and appropriately normed (SCAN). The closest agreement is achieved by the simplest density functional LDA, and by the most sophisticated one, SCAN+rVV10 (Vydrov-Van Voorhis 2010). The long-range van der Waals interaction, incorporated through rVV10, increases the surface energies by about 10%, and increases the work functions by about 3%. LDA works for metal surfaces through two known error cancellations. The Perdew-Burke-Ernzerhof generalized gradient approximation tends to underestimate both surface energies (by about 24%) and work functions (by about 4%), yielding the least-accurate results. The amount by which a functional underestimates these surface properties correlates with the extent to which it neglects van der Waals attraction at intermediate and long range. Qualitative arguments are given for the signs of the van der Waals contributions to the surface energy and work function. A standard expression for the work function in Kohn-Sham (KS) theory is shown to be valid in generalized KS theory. Interlayer relaxations from different functionals are in reasonable agreement with one another, and usually with experiment.