The curious case of proton migration under pressure in the malonic acid and 4,4'-bipyridine cocrystal.

In the search for new active pharmaceutical ingredients, the precise control of the chemistry of cocrystals becomes essential. One crucial step within this chemistry is proton migration between cocrystal coformers to form a salt, usually anticipated by the empirical ΔpKa rule. Due to the effective role it plays in modifying intermolecular distances and interactions, pressure adds a new dimension to the ΔpKa rule. Still, this variable has been scarcely applied to induce proton-transfer reactions within these systems. In our study, high-pressure X-ray diffraction and Raman spectroscopy experiments, supported by DFT calculations, reveal modifications to the protonation states of the 4,4'-bipyridine (BIPY) and malonic acid (MA) cocrystal (BIPYMA) that allow the conversion of the cocrystal phase into ionic salt polymorphs. On compression, neutral BIPYMA and monoprotonated (BIPYH+MA-) species coexist up to 3.1 GPa, where a phase transition to a structure of P21/c symmetry occurs, induced by a double proton-transfer reaction forming BIPYH22+MA2-. The low-pressure C2/c phase is recovered at 2.4 GPa on decompression, leading to a 0.7 GPa hysteresis pressure range. This is one of a few studies on proton transfer in multicomponent crystals that shows how susceptible the interconversion between differently charged species is to even slight pressure changes, and how the proton transfer can be a triggering factor leading to changes in the crystal symmetry. These new data, coupled with information from previous reports on proton-transfer reactions between coformers, extend the applicability of the ΔpKa rule incorporating the pressure required to induce salt formation.

A mechanism for aragonite to post-aragonite transition in MCO3 (M = Ca, Sr and Ba) carbonates: evidence of a hidden metastable polymorph.

To advance in the understanding of the Earth's carbon cycle, it is necessary to determine thermodynamic boundaries and kinetic barriers associated with the pressure-induced polymorphic sequence of alkaline-earth carbonates. Following a symmetry-based strategy within the martensitic approximation, we propose a two-step mechanism mediated by a hexagonal P63/mmc structure for the aragonite to post-aragonite transformation in the MCO3 (M = Ca, Sr, Ba) crystal family. The calculated transition pressures and activation energies, from ∼7 to 42 GPa and ∼0.3 to 0.6 eV, respectively, are low enough to allow this transformation to occur under mantle conditions. Our analysis reveals that the intermediate hexagonal structure is the early one proposed by Holl et al., Phys. Chem. Miner., 2000, 27, 467-473 for high pressure BaCO3, and later considered as metastable. Phonon calculations inform that this P63/mmc structure is in fact unstable at zero pressure. Remarkably, our molecular dynamics calculations showed that this instability smoothly leads to a dynamically stable P63mc structure, which we confirm is actually the phase observed by Holl et al. This finding allows us to reconcile previous controversial data and contributes to clarifying the role of carbonates in the Earth's interior.

A Chemical-Pressure-Induced Phase Transition Controlled by Lone Electron Pair Activity.

The chemical pressure approach offers a new paradigm for property control in functional materials. In this work, we disclose a correlation between the ß â†’ α pressure-induced phase transition in SnMoO4 and the substitution process of Mo6+ by W6+ in SnMo1-xWxO4 solid solutions (x = 0-1). Special attention is paid to discriminating the role of the lone pair Sn2+ cation from the structural distortive effect along the Mo/W substitution process, which is crucial to disentangle the driven force of the transition phase. Furthermore, the reverse α → ß transition observed at high temperature in SnWO4 is rationalized on the same basis as a negative pressure effect associated with a decreasing of W6+ percentage in the solid solution. This work opens a versatile chemical approach in which the types of interactions along the formation of solid solutions are clearly differentiated and can also be used to tune their properties, providing opportunities for the development of new materials.

High-Pressure and High-Temperature Chemistry of Phosphorus and Nitrogen: Synthesis and Characterization of α- and γ-P3N5.

The direct chemical reactivity between phosphorus and nitrogen was induced under high-pressure and high-temperature conditions (9.1 GPa and 2000-2500 K), generated by a laser-heated diamond anvil cell and studied by synchrotron X-ray diffraction, Raman spectroscopy, and DFT calculations. α-P3N5 and γ-P3N5 were identified as reaction products. The structural parameters and vibrational frequencies of γ-P3N5 were characterized as a function of pressure during room-temperature compression and decompression to ambient conditions, determining the equation of state of the material up to 32.6 GPa and providing insight about the lattice dynamics of the unit cell during compression, which essentially proceeds through the rotation of the PN5 square pyramids and the distortion of the PN4 tetrahedra. Although the identification of α-P3N5 demonstrates for the first time the direct synthesis of this compound from the elements, its detection in the outer regions of the laser-heated area suggests α-P3N5 as an intermediate step in the progressive nitridation of phosphorus toward the formation of γ-P3N5 with increasing coordination number of P by N from 4 to 5. No evidence of a higher-pressure phase transition was observed, excluding the existence of predicted structures containing octahedrally hexacoordinated P atoms in the investigated pressure range.

Pressure and temperature stability boundaries of cubic SiC polymorphs: a first-principles investigation.

A better understanding of the effects of temperature and pressure on the wide gap SiC semiconductor is necessary for both (i) an improvement of the performance of this compound in a variety of technological applications and (ii) a clarification of controversial issues related to the stability of its cubic polymorphs at high pressure and high temperature. Bearing in mind this double demand, we perform first-principles calculations of the phonon band structures, vibrational density of states, and thermal and mode Grüneisen parameters of the zinc blende (B3) and rock-salt (B1) cubic polymorphs of 3C-SiC covering pressures and temperatures up to 120 GPa and 3000 K, respectively. Under a martensitic description of the B3-B1 transformation, we found that the large hysteresis pressure range observed at room temperature (35-100 GPa) disappears at around 1100 K. The calculated Clapeyron slope of this transformation is slightly negative, with average values of -2.9 MPa K-1 in the 0-3000 K interval and -3.7 MPa K-1 at 2000 K. We also study the decomposition reaction of the two cubic polymorphs into their elemental constituents (C and Si), finding a decreasing (increasing) decomposition temperature for the B3 (B1) phase as the pressure is applied. All these results are sustained by good agreement with other recently reported experimental and theoretical thermodynamic data that have also been evaluated under our quasi-harmonic approximation framework.

Structural and Electronic Effect Driven Distortions in Visible Light Absorbing Polar Materials ATa2V2O11 (A = Sr, Pb).

Complex vanadates of tantalum(V), such as ATa2V2O11 (A = Sr, Pb), are rare and underrated materials, which have potential application domains that could be substantially expanded, mitigating the existing controversy on their atomic and electronic organization. Herein, we present a thorough structural examination combining synchrotron powder X-ray diffraction-aided distortion mode analysis with computational methods to study hettotypes of SrTa2V2O11 (STVO) and PbTa2V2O11 (PTVO). Being distinct from the perovskite family due to the presence of [VO4] groups, both compounds are polar dielectric materials with certain similarities to SBT and PBT Aurivillius phases. Applying the model of anions of metallic matrices to the analysis of electron localization functions calculated on top of as-established equilibrium structures helps retrace the effects in the Sr and Pb surroundings on the respective crystal packings of STVO and PTVO.

Comment on "Uncommon structural and bonding properties in Ag16B4O10" by A. Kovalevskiy, C. Yin, J. Nuss, U. Wedig, and M. Jansen, Chem. Sci., 2020, 11, 962.

A thorough systematic study of the Electron Localization Function (ELF) in fcc silver metal, the deficient vacant-type Ag16□4 structure, and the Ag16B4O10 title compound of the Chem. Sci., 2020, 11, 962 edge article leads to a further understanding of the sub-valent characteristics of silver in the silver borate compound. By visualizing the process in three consecutive steps, (fcc)eq-Ag → (fcc)ex-Ag → Ag16□4 → Ag16B4O10, the electron reduction of Ag atoms can be traced to be due to (i) the expansion (ex) of the host metallic array from its equilibrium (eq) geometry and (ii) the vacancy creation and subsequent insertion of guest borate clusters. Our ELF analysis also allows us to identify to what extent metallic features remain in the title compound, providing an alternative explanation of why Ag16B4O10 is not a conductor whereas pure silver is.

Cyclic Single Atom Vertical Manipulation on a Nonmetallic Surface.

Motivated by the quest for experimental procedures capable of controlled manipulation of single atoms on surfaces, we set up a computational strategy that explores the cyclical vertical manipulation of a broad set of single atoms on the GaAs(110) surface. First-principles simulations of atomic force microscope tip-sample interactions were performed considering families of GaAs and Au-terminated tip apexes with varying crystalline termination. We identified a subset of tips capable of both picking up and depositing an adatom (Ga, As, Al, and Au) any number of times via a modify-restore cycle that "resets" the apex of the scanning probe to its original structure at the end of each cycle. Manipulation becomes successful within a certain window of lateral and vertical tip distances that are observed to be different for extracting and depositing each atom. A practical experimental protocol of special utility for potential cyclical manipulation of single atoms on a nonmetallic surface is proposed.

Highs and Lows of Bond Lengths: Is There Any Limit?

Two distinct points on the potential energy curve (PEC) of a pairwise interaction, the zero-energy crossing point and the point where the stretching force constant vanishes, allow us to anticipate the range of possible distances between two atoms in diatomic, molecular moieties and crystalline systems. We show that these bond-stability boundaries are unambiguously defined and correlate with topological descriptors of electron-density-based scalar fields, and can be calculated using generic PECs. Chemical databases and quantum-mechanical calculations are used to analyze a full set of diatomic bonds of atoms from the s-p main block. Emphasis is placed on the effect of substituents in C-C covalent bonds, concluding that distances shorter than 1.14 Šor longer than 2.0 Šare unlikely to be achieved, in agreement with ultra-high-pressure data and transition-state distances, respectively. Presumed exceptions are used to place our model in the correct framework and to formulate a conjecture for chained interactions, which offers an explanation for the multimodal histogram of O-H distances reported for hundreds of chemical systems.

Pressure-Driven Metallization in Hafnium Diselenide.

The quest for new transition metal dichalcogenides (TMDs) with outstanding electronic properties operating under ambient conditions draws us to investigate the 1T-HfSe2 polytype under hydrostatic pressure. Diamond anvil cell (DAC) devices coupled to in situ synchrotron X-ray, Raman, and optical (VIS-NIR) absorption experiments along with density functional theory (DFT)-based calculations prove that (i) bulk 1T-HfSe2 exhibits strong structural and vibrational anisotropies, being the interlayer direction especially sensitive to pressure changes, (ii) the indirect gap of 1T-HfSe2 tends to vanish by a -0.1 eV/GPa pressure rate, slightly faster than MoS2 or WS2, (iii) the onset of the metallic behavior appears at Pmet ∼10 GPa, which is to date the lowest pressure among common TMDs, and finally, (iv) the electronic transition is explained by the bulk modulus B0-Pmet correlation, along with the pressure coefficient of the band gap, in terms of the electronic overlap between chalcogenide p-type and metal d-type orbitals. Overall, our findings identify 1T-HfSe2 as a new efficient TMD material with potential multipurpose technological applications.

Disclosing the behavior under hydrostatic pressure of rhombohedral MgIn2Se4 by means of first-principles calculations.

AM2X4 crystalline materials display important technological electronic, optical and magnetic properties that are sensitive to general stress effects. In this paper, the behavior under hydrostatic pressure of the ambient condition rhombohedral phase of MgIn2Se4 is investigated in detail for the first time. We carried out first-principles calculations within the density functional theory framework aimed at determining the pressure-induced polymorphic sequence of this selenide. To accurately evaluate transition pressures at room temperature, thermal corrections have been included after the computation of phonon dispersion curves in potential candidate phases, namely the initial rhombohedral R3[combining macron]m one, inverse and direct spinels, LiTiO2-type and defective I4[combining macron] structures. Only the transition from the R3[combining macron]m to the inverse spinel phase was found to fulfill the thermodynamic and mechanical stability criteria. The direct spinel could appear as metastable if kinetic effects hinder the above transition. Additionally, electronic band structures and chemical bonding properties were analyzed from the outcome of our quantum-mechanical solutions reporting band gap values and ionicity and noncovalent interaction indexes. It is shown that the investigated compound keeps behaving as a semiconductor, loses its van der Waals interactions, and becomes more covalent as hydrostatic pressure is applied.

Generalized Stress-Redox Equivalence: A Chemical Link between Pressure and Electronegativity in Inorganic Crystals.

The crystal structure of many inorganic compounds can be understood as a metallic matrix playing the role of a host lattice in which the nonmetallic atomic constituents are located, the Anions in Metallic Matrices (AMM) model stated. The power and utility of this model lie in its capacity to anticipate the actual positions of the guest atoms in inorganic crystals using only the information known from the metal lattice structure. As a pertinent test-bed for the AMM model, we choose a set of common metallic phases along with other nonconventional or more complex structures (face-centered cubic (fcc) and simple cubic Ca, CsCl-type BaSn, hP4-K, and fcc-Na) and perform density functional theory electronic structure calculations. Our topological analysis of the chemical pressure (CP) scalar field, easily derived from these standard first-principles electronic computations, reveals that CP minima appear just at the precise positions of the nonmetallic elements in typical inorganic crystals presenting the above metallic subarrays: CaF2, rock-salt and CsCl-type phases of CaX (X = O, S, Se, Te), BaSnO3, K2S, and NaX (X = F, Cl, Br, I). A theoretical basis for this correlation is provided by exploring the equivalence between hydrostatic pressure and the oxidation (or reduction) effect induced by the nonmetallic element on the metal structure. Indeed, our CP analysis leads us to propose a generalized stress-redox equivalence that is able to account for the two main observed phenomena in solid inorganic compounds upon crystal formation: (i) the expansion or contraction experienced by the metal structure after hosting the nonmetallic element while its topology is maintained and (ii) the increasing or decreasing of the effective charge associated with the anions in inorganic compounds with respect to the charge already present in the interstices of the metal network. We demonstrate that a rational explanation of this rich behavior is provided by means of Pearson-Parr's electronegativity equalization principle.

Computational Modeling of Tensile Stress Effects on the Structure and Stability of Prototypical Covalent and Layered Materials.

Understanding the stability limit of crystalline materials under variable tensile stress conditions is of capital interest for technological applications. In this study, we present results from first-principles density functional theory calculations that quantitatively account for the response of selected covalent and layered materials to general stress conditions. In particular, we have evaluated the ideal strength along the main crystallographic directions of 3C and 2H polytypes of SiC, hexagonal ABA stacking of graphite and 2H-MoS 2 . Transverse superimposed stress on the tensile stress was taken into account in order to evaluate how the critical strength is affected by these multi-load conditions. In general, increasing transverse stress from negative to positive values leads to the expected decreasing of the critical strength. Few exceptions found in the compressive stress region correlate with the trends in the density of bonds along the directions with the unexpected behavior. In addition, we propose a modified spinodal equation of state able to accurately describe the calculated stress-strain curves. This analytical function is of general use and can also be applied to experimental data anticipating critical strengths and strain values, and for providing information on the energy stored in tensile stress processes.

Controlling the preferential motion of chiral molecular walkers on a surface.

Molecular walkers standing on two or more "feet" on an anisotropic periodic potential of a crystal surface may perform a one-dimensional Brownian motion at the surface-vacuum interface along a particular direction in which their mobility is the largest. In thermal equilibrium the molecules move with equal probabilities both ways along this direction, as expected from the detailed balance principle, well-known in chemical reactivity and in the theory of molecular motors. For molecules that possess an asymmetric potential energy surface (PES), we propose a generic method based on the application of a time-periodic external stimulus that would enable the molecules to move preferentially in a single direction thereby acting as Brownian ratchets. To illustrate this method, we consider a prototypical synthetic chiral molecular walker, 1,3-bis(imidazol-1-ylmethyl)-5(1-phenylethyl)benzene, diffusing on the anisotropic Cu(110) surface along the Cu rows. As unveiled by our kinetic Monte Carlo simulations based on the rates calculated using ab initio density functional theory, this molecule moves to the nearest equivalent lattice site via the so-called inchworm mechanism in which it steps first with the rear foot and then with the front foot. As a result, the molecule diffuses via a two-step mechanism, and due to its inherent asymmetry, the corresponding PES is also spatially asymmetric. Taking advantage of this fact, we show how the external stimulus can be tuned to separate molecules of different chirality, orientation and conformation. The consequences of these findings for molecular machines and the separation of enantiomers are also discussed.

Correction: Controlling the preferential motion of chiral molecular walkers on a surface.

[This corrects the article DOI: 10.1039/C9SC01135H.].

Chemical Pressure Maps of Molecules and Materials: Merging the Visual and Physical in Bonding Analysis.

The characterization of bonding interactions in molecules and materials is one of the major applications of quantum mechanical calculations. Numerous schemes have been devised to identify and visualize chemical bonds, including the electron localization function, quantum theory of atoms in molecules, and natural bond orbital analysis, whereas the energetics of bond formation are generally analyzed in qualitative terms through various forms of energy partitioning schemes. In this Article, we illustrate how the chemical pressure (CP) approach recently developed for analyzing atomic size effects in solid state compounds provides a basis for merging these two approaches, in which bonds are revealed through the forces of attraction and repulsion acting between the atoms. Using a series of model systems that include simple molecules (H2, CO2, and S8), extended structures (graphene and diamond), and systems exhibiting intermolecular interactions (ice and graphite), as well as simple representatives of metallic and ionic bonding (Na and NaH, respectively), we show how CP maps can differentiate a range of bonding phenomena. The approach also allows for the partitioning of the potential and kinetic contributions to the interatomic interactions, yielding schemes that capture the physical model for the chemical bond offered by Ruedenberg and co-workers.

The role of isomerization in the kinetics of self-assembly: p-terphenyl-m-dicarbonitrile on the Ag(111) surface.

Using a toolkit of theoretical techniques comprising ab initio density functional theory calculations, the nudged elastic band method and kinetic Monte Carlo (KMC) modeling, we investigate in great detail how para-terphenyl-meta-dicarbonitrile (pTmDC) molecules diffuse and isomerize to self-assemble on the Ag(111) surface. We show that molecules "walk" on the surface via a pivoting mechanism moving each of its two "legs" one at a time. We then identify a peculiar "under-side" isomerization mechanism capable of changing the molecules chirality, and demonstrate that it is fundamental in understanding the growth of hydrogen bonding assembles of ribbons, linkers, clusters and brickwall islands on the Ag(111) surface, as observed in recent scanning tunneling microscopy experiments (ChemPhysChem, 2010, 11, 1446). The discovered underlying atomistic mechanism of self-assembly may be behind the growth of other hydrogen bonding structures of chiral molecules on metal surfaces.

Anions in metallic matrices model: application to the aluminium crystal chemistry.

We introduce and discuss an interpretative model of the structure and bonding of inorganic crystals containing metallic elements. The central idea is the conception of the crystal structure of such an inorganic compound as a metallic matrix whose geometric and electronic structures govern the formation and localization of the anions in the lattice. This is the reason for labelling the model anions in metallic matrices (AMM). Taking the AlX3 crystal family (X = F, Cl, OH) as a suitable test-bed class of compounds, we illustrate how this approach gives a direct interpretation of the crystalline structures and explains the variable coordination that Al exhibits in crystalline materials. An exhaustive analysis of the topology of the electron density allows us to provide a quantum-mechanical assessment of the main hypotheses of the AMM model and to uncover, using microscopic arguments, the behavior of anions as chemical pressure agents.