Reply to ‘comment on “extending Hirshfeld-I to bulk and periodic materials”’

UV-Curable Biobased Polyacrylates Based on a Multifunctional Monomer Derived from Furfural

The controlled polymerization of a new biobased monomer, 4-oxocyclopent-2-en-1-yl acrylate (4CPA), was established via reversible addition-fragmentation chain transfer (RAFT) (co)polymerization to yield polymers bearing pendent cyclopentenone units. 4CPA contains two reactive functionalities, namely, a vinyl group and an internal double bond, and is an unsymmetrical monomer. Therefore, competition between the internal double bond and the vinyl group eventually leads to gel formation. With RAFT polymerization, when aiming for a degree of polymerization (DP) of 100, maximum 4CPA conversions of the vinyl group between 19.0 and 45.2% were obtained without gel formation or extensive broadening of the dispersity. When the same conditions were applied in the copolymerization of 4CPA with lauryl acrylate (LA), methyl acrylate (MA), and isobornyl acrylate, 4CPA conversions of the vinyl group between 63 and 95% were reached. The additional functionality of 4CPA in copolymers was demonstrated by model studies with 4-oxocyclopent-2-en-1-yl acetate (1), which readily dimerized under UV light via [2 + 2] photocyclodimerization. First-principles quantum mechanical simulations supported the experimental observations made in NMR. Based on the calculated energetics and chemical shifts, a mixture of head-to-head and head-to-tail dimers of (1) were identified. Using the dimerization mechanism, solvent-cast LA and MA copolymers containing 30 mol % 4CPA were cross-linked under UV light to obtain thin films. The cross-linked films were characterized by dynamic scanning calorimetry, dynamic mechanical analysis, IR, and swelling experiments. This is the first case where 4CPA is described as a monomer for functional biobased polymers that can undergo additional UV curing via photodimerization.

Investigation of structural, electronic and magnetic properties of breathing metal-organic framework MIL-47(Mn): a first principles approach

The structural, electronic and magnetic properties of the MIL-47(Mn) metal-organic framework are investigated using first principles calculations. We find that the large-pore structure is the ground state of this material. We show that upon transition from the large-pore to the narrow-pore structure, the magnetic ground-state configuration changes from antiferromagnetic to ferromagnetic, consistent with the computed values of the intra-chain coupling constant. Furthermore, the antiferromagnetic and ferromagnetic configuration phases have intrinsically different electronic behavior: the former is semiconducting, the latter is a metal or half-metal. The change of electronic properties during breathing posits MIL-47(Mn) as a good candidate for sensing and other applications. Our calculated electronic band structure for MIL-47(Mn) presents a combination of flat dispersionless and strongly dispersive regions in the valence and conduction bands, indicative of quasi-1D electronic behavior. The spin coupling constants are obtained by mapping the total energies onto a spin Hamiltonian. The inter-chain coupling is found to be at least one order of magnitude smaller than the intra-chain coupling for both large and narrow pores. Interestingly, the intra-chain coupling changes sign and becomes five times stronger going from the large pore to the narrow pore structure. As such MIL-47(Mn) could provide unique opportunities for tunable low-dimensional magnetism in transition metal oxide systems.

Vertical ferroelectric switching by in-plane sliding of two-dimensional bilayer WTe2

Based on first-principles calculations, we studied the ferroelectric properties of bilayer 1T'-WTe₂. In this work, we discovered that the polarization stems from uncompensated out-of-plane interlayer charge transfer, which can be switched upon interlayer sliding of an in-plane translation. Our differential charge density results also confirmed that such ferroelectricity in the bilayer WTe₂ is derived from interlayer charge transfer. The ferroelectric polarization directions further control the spin texture of the bilayer WTe₂, which may have important applications in spintronics. Therefore, we propose a spin field effect transistor (spin-FET) design that may effectively improve the spin-polarized injection rate. In addition, the lattice strain has been found to have an important influence on the ferroelectric properties of the bilayer WTe₂. One can effectively increase the polarization with a maximum at 3% tensile strain, whereas a 3% compressive strain can transform the bilayer WTe₂ from the ferroelectric to paraelectric phase.

New scaling relations to compute atom-in-material polarizabilities and dispersion coefficients: part 1. Theory and accuracy

Polarizabilities and London dispersion forces are important to many chemical processes. Force fields for classical atomistic simulations can be constructed using atom-in-material polarizabilities and C n (n = 6, 8, 9, 10…) dispersion coefficients. This article addresses the key question of how to efficiently assign these parameters to constituent atoms in a material so that properties of the whole material are better reproduced. We develop a new set of scaling laws and computational algorithms (called MCLF) to do this in an accurate and computationally efficient manner across diverse material types. We introduce a conduction limit upper bound and m-scaling to describe the different behaviors of surface and buried atoms. We validate MCLF by comparing results to high-level benchmarks for isolated neutral and charged atoms, diverse diatomic molecules, various polyatomic molecules (e.g., polyacenes, fullerenes, and small organic and inorganic molecules), and dense solids (including metallic, covalent, and ionic). We also present results for the HIV reverse transcriptase enzyme complexed with an inhibitor molecule. MCLF provides the non-directionally screened polarizabilities required to construct force fields, the directionally-screened static polarizability tensor components and eigenvalues, and environmentally screened C6 coefficients. Overall, MCLF has improved accuracy compared to the TS-SCS method. For TS-SCS, we compared charge partitioning methods and show DDEC6 partitioning yields more accurate results than Hirshfeld partitioning. MCLF also gives approximations for C8, C9, and C10 dispersion coefficients and quantum Drude oscillator parameters. This method should find widespread applications to parameterize classical force fields and density functional theory (DFT) + dispersion methods.

Ligand Addition Energies and the Stoichiometry of Colloidal Nanocrystals

Experimental nonstoichiometries of colloidal nanocrystals such as CdSe and PbS are accounted for by attributing to each constituent atom and capping ligand a formal charge equal to its most common oxidation state to obtain an overall neutral nanocrystal. In spite of its apparent simplicity, little theoretical support of this approach-called here the oxidation-number sum rule-is present in the current literature. Here, we introduce the ligand addition energy, which we define as the energy gained or expended upon the transfer of one ligand from a reference state to a metal-rich solid surface. For the combination of CdSe, ZnSe and InP with either chalcogen, halogen or hydrochalcogen ligands, we compute successive ligand addition energies using ab initio methods and determine the thermodynamically stable surface composition as that composition where ligand addition turns endothermic. We find that the oxidation-number sum rule is valid in many situations, although exceptions occur for each material studied, most notably when exposed to small oxidative ligands. In the case of InP, however, violations are more severe, extending toward the entire chalcogen ligand family. In addition, we find that electronegativity rather than chemical hardness is a reasonable predictor for ligand addition energies, with the most electronegative ligands yielding the most exothermic addition energies. Finally, we argue that the ligand addition energy will be a most useful quantity for future computational studies on the structure, stability and reactivity of nanocrystal surfaces.

Introducing DDEC6 atomic population analysis: part 1. Charge partitioning theory and methodology

We introduce a new atomic population analysis method that performs exceptionally well across an extremely broad range of periodic and non-periodic material types.

Introducing DDEC6 atomic population analysis: part 2. Computed results for a wide range of periodic and nonperiodic materials

DDEC6 atomic population analysis gives excellent performance for small and large molecules, porous solids, dense solids, solid surfaces, organometallic complexes, nanoclusters, and magnetic materials.

A Flexible Photoactive Titanium Metal-Organic Framework Based on a [Ti(IV)3(μ3-O)(O)2(COO)6] Cluster

The synthesis of titanium-carboxylate metal-organic frameworks (MOFs) is hampered by the high reactivity of the commonly employed alkoxide precursors. Herein, we present an innovative approach to titanium-based MOFs by the use of titanocene dichloride to synthesize COK-69, the first breathing Ti MOF, which is built up from trans-1,4-cyclohexanedicarboxylate linkers and an unprecedented [Ti(IV)3(μ3-O)(O)2(COO)6] cluster. The photoactive properties of COK-69 were investigated in depth by proton-coupled electron-transfer experiments, which revealed that up to one Ti(IV) center per cluster can be photoreduced to Ti(III) while preserving the structural integrity of the framework. The electronic structure of COK-69 was determined by molecular modeling, and a band gap of 3.77 eV was found.

Quasi-1D physics in metal-organic frameworks: MIL-47(V) from first principles

The geometric and electronic structure of the MIL-47(V) metal-organic framework (MOF) is investigated by using ab initio density functional theory (DFT) calculations. Special focus is placed on the relation between the spin configuration and the properties of the MOF. The ground state is found to be antiferromagnetic, with an equilibrium volume of 1554.70 Å(3). The transition pressure of the pressure-induced large-pore-to-narrow-pore phase transition is calculated to be 82 MPa and 124 MPa for systems with ferromagnetic and antiferromagnetic chains, respectively. For a mixed system, the transition pressure is found to be a weighted average of the ferromagnetic and antiferromagnetic transition pressures. Mapping DFT energies onto a simple-spin Hamiltonian shows both the intra- and inter-chain coupling to be antiferromagnetic, with the latter coupling constant being two orders of magnitude smaller than the former, suggesting the MIL-47(V) to present quasi-1D behavior. The electronic structure of the different spin configurations is investigated and it shows that the band gap position varies strongly with the spin configuration. The valence and conduction bands show a clear V d-character. In addition, these bands are flat in directions orthogonal to VO6 chains, while showing dispersion along the the direction of the VO6 chains, similar as for other quasi-1D materials.

Modeling 1D structures on semiconductor surfaces: synergy of theory and experiment

Atomic scale nanowires attract enormous interest in a wide range of fields. On the one hand, due to their quasi-one-dimensional nature, they can act as an experimental testbed for exotic physics: Peierls instability, charge density waves, and Luttinger liquid behavior. On the other hand, due to their small size, they are of interest not only for future device applications in the micro-electronics industry, but also for applications regarding molecular electronics. This versatile nature makes them interesting systems to produce and study, but their size and growth conditions push both experimental production and theoretical modeling to their limits. In this review, modeling of atomic scale nanowires on semiconductor surfaces is discussed, focusing on the interplay between theory and experiment. The current state of modeling efforts on Pt- and Au-induced nanowires on Ge(001) is presented, indicating their similarities and differences. Recently discovered nanowire systems (Ir, Co, Sr) on the Ge(001) surface are also touched upon. The importance of scanning tunneling microscopy as a tool for direct comparison of theoretical and experimental data is shown, as is the power of density functional theory as an atomistic simulation approach. It becomes clear that complementary strengths of theoretical and experimental investigations are required for successful modeling of the atomistic nanowires, due to their complexity.