Pressure-induced reconstructive phase transitions, polarization with metallicity, and enhanced hardness in antiperovskite MgCNi3.

In general, hydrostatic pressure can suppress electrical polarization, instead of creating and/or enhancing polarization like strain engineering. Here, a combination of first-principles calculations and CALYPSO crystal structures prediction is used to point out that hydrostatic pressure applied on antiperovskite MgCNi3 can stabilize polarization with metallicity, and thus a polar metal can exist under high pressure. Strikingly, the metallic polar phase of MgCNi3 exhibits an original linear-cubic coupling between polar and nonpolar modes, resulting in an asymmetrical double-well when the polarization is switched. Moreover, another novel phase of MgCNi3 under high pressure possesses an enhanced hardness stemming from a robust s-s electrons interaction of an unexpected C-C bond, rather than typical sp3 orbital hybridization. These discoveries open new routes to design superhard materials and polar metals.

Structural Stability and Evolution of Medium-Sized Tantalum-Doped Boron Clusters: A Half-Sandwich-Structured TaB12- Cluster.

Transition-metal (TM)-doped boron clusters have received considerable attention in recent years, in part, because of their remarkable size-dependent structural and electronic properties. However, the structures of medium-sized boron clusters doped with TM atoms are still not well-known because of the much increased complexity of the potential surface as well as the rapid increase in the number of low-energy isomers, which are the challenges in cluster structural searches. Here, by means of an unbiased structure search, we systematically investigated the structural evolution of medium-sized tantalum-doped boron clusters, TaBn0/- (n = 10-20). The results revealed that TaBn0/- (n = 10-15) clusters adopt half-sandwich molecular geometries, with the notable exception of TaB10-, while for n = 16-18 and 19-20, the lowest-energy clusters are characterized by drum-type geometries and tubular molecules with two B atoms on the top, respectively. Good agreement between the calculated and experimental photoelectron spectra strongly support the validity of our global minimum structures. Molecular orbital and adaptive natural density partitioning analyses indicate that the enhanced stability of half-sandwich TaB12- is due to the strong interaction of the Ta atom (5d orbitals) with surrounding B atoms (2p orbitals) and σ B-B bonds in the B12 moiety.

Evolution of the Structural and Electronic Properties of Medium-Sized Sodium Clusters: A Honeycomb-Like Na20 Cluster.

Sodium is one of the best examples of a free-electron-like metal and of a certain technological interest. However, an unambiguous determination of the structural evolution of sodium clusters is challenging. Here, we performed an unbiased structure search among neutral and anionic sodium clusters in the medium size range of 10-25 atoms, using the Crystal structure AnaLYsis by Particle Swarm Optimization (CALYPSO) method. Geometries are determined by CALYPSO structure searches, followed by reoptimization of a large number of candidate structures. For most cluster sizes the simulated photoelectron spectra of the lowest-energy structures are in excellent agreement with the experimental data, indicating that the current ground-state structures are the true minima. The equilibrium geometries show that, for both neutral and anionic species, the structural evolution from bilayer structures to layered outsides with interior atoms occurs at n = 16. A novel unprecedented honeycomb-like structure of Na20 cluster with C3 symmetry is uncovered, which is more stable than the prior suggested structure based on pentagonal structural motifs.

Prediction of hypervalent molecules: investigation on MnC (M = Li, Na, K, Rb and Cs; n = 1-8) clusters.

New hypervalent molecules have emerged from a systematic exploration of the structure and bonding of MnC (M = Li, Na, K, Rb and Cs; n = 1-8) clusters via an unbiased CALYPSO structure investigation combined with density functional theory. The global minimum structures are obtained at the B3LYP/6-311+G* and CCSD(T)/6-311+G* levels of theory. The observed growth behavior clearly indicates that the ground state of MnC (M = Li, Na, K, Rb and Cs; n = 1-8) is transformed from a planar to a three-dimensional (3D) structure at n = 4. A maximum of six alkali atoms can be bound atomically to a carbon atom. The determination of the averaged binding energies Eb(n), fragmentation energies ΔE(n) and HOMO-LUMO energy gaps unambiguously supports the stability of M6C. This indicated conclusively that 6 is a magic Li-coordination number for C. The nature of bonding is further investigated by an insightful analysis of the highest occupied molecular orbital (HOMO) and the topology of chemical bonds for the most stable clusters. In the final step, electron localization functions (ELF) and density of states (DOS) are determined in order to consolidate the acquired information on the studied electronic structures.

Large polarization and dielectric response in epitaxial SrZrO3 films.

First-principles calculations are performed to investigate the ferroelectric and dielectric properties of (001) epitaxial SrZrO3 thin films under misfit strain. A rich phase diagram is predicted. By condensing the polar instability, the ferroelectric Pmc21 and Ima2 phases can coexist under tensile strain (about 3.7%-5.2%/5.7%). Combining in-plane ferroelectric (FExy) and out-of-plane in-phase antiferrodistortive (IAFDz) modes, another new Pmc21 state (P > 56 µC cm(-2)) occurs with increase in the tensile strain. The paraelectric I4/mcm and ferroelectric P4mm phases emerge around -3.2%/-3.7% and -6.4%/-7.4% compressive strain, respectively. The former exhibits an intense out-of-plane dielectric response, while the latter possesses a rather large polarization (∼ 110 µC cm(-2)). The large polarization and dielectric response are discussed in relationship to strain-driven structural distortion.

Structural and relative stabilities, electronic properties, and hardness of iron tetraborides from first prinicples.

First-principles calculations were carried out to investigate the structure, phase stability, electronic property, and roles of metallicity in the hardness for recently synthesized FeB4 with various different structures. Our calculation indicates that the orthorhombic phase with Pnnm symmetry is the most energetically stable one. The other four new dynamically stable phases belong to space groups monoclinic C2/m, orthorhombic Pmmn, trigonal R3̅m, and hexagonal P63/mmc. Their mechanical and thermodynamic stabilities are verified by calculating elastic constants, formation enthalpies, and phonon dispersions. We found that all phases are stabilized further under pressure. Above the pressure of about 50 GPa, the formation enthalpy of Pmmn is almost equal to that of P63/mmc phase. The analysis on density of states not only demonstrates that formation of strong covalent bonding in these compounds contributes greatly to their stabilities but also that they all exhibit metallic behavior which does not relate to the approach used. By considering metallic contributions, the estimated Vickers hardness values based on the semiempirical model show that the OsB4-structured FeB4, with a hardness of 48.1 GPa, well exceeding the limitation of superhardness (40 GPa), is more hard than the most stable phase. The others are predicted to be potential hard materials. Moreover, the atomic configuration and strong B-B covalent bonds are found to play important roles in the hardness of materials.

Probing the structural, bonding, and magnetic properties of cobalt coordination complexes: co-benzene, co-pyridine, and co-pyrimidine.

Neutral and anionic Co1,2(benzene)1,2, Co1,2(pyridine)1,2, and Co1,2(pyrimidine)1,2 complexes have been investigated within the framework of an all-electron gradient-corrected density functional theory. The ground-state structures for each size clusters were identified based on the geometry optimization. Meanwhile, their electron affinities and vertical detachment energies were predicted and compared with the experimental values. By analyzing the pattern of highest occupied molecular orbitals (HOMOs), we found that the bond formation of these Co-organic complexes mainly arises from the 3d/4s electrons of the cobalt atoms and the π-cloud of the organic molecules. More importantly, we presented an approach to map and analyze the Co-organic interactions from another perspective. The scatter plots of the reduced density gradient (RDG) versus ρ allow us to identify the different types of interactions, and the maps of the gradient isosurfaces show a rich visualization of chemical bond and steric effects. Their magnetic properties were studied by determining the spin magnetic moments and visualizing the spin density distributions. Finally, the natural population analysis (NPA) charge was calculated to achieve a deep insight into the distribution of electron density and the reliable charge-transfer information.

Structures, electrophilic properties, and hydrogen bonds of cytidine, uridine, and their radical anions: Microhydration effects.

Structures, electrophilic properties, and hydrogen bonds of the neutral and anionic monohydrated nucleoside, (cytidine)H2O, and (uridine)H2O have been systematically investigated using density functional theory. Various water-binding sites were predicted by explicitly considering the optimized monohydrated structures. Meanwhile, predictions of electron affinities and vertical detachment energies were also carried out to investigate their electrophilic properties. By examining the singly occupied molecular orbital and natural population analysis, we found the excess negative charge is localized on the cytidine and uridine moiety in anionic monohydrates. This may be the reason why the strength of hydrogen bonding undergoes an obvious change upon the extra electron attachment. Based on the electron density (ρ) and reduced density gradient (RDG), we present an approach to map and analyze the weak interaction (especially hydrogen bond) in monohydrated cytidine and uridine. The scatter plots of RDG versus ρ allow us to identify the different type interactions. Meanwhile, the maps of the gradient isosurfaces show a rich visualization of hydrogen bond, van der Waals interaction, and steric effect.

Formation and properties of iron-based magnetic superhalogens: a theoretical study.

In order to explore new magnetic superhalogens, we have systematically investigated the structures, electrophilic properties, stabilities, magnetic properties, and fragmentation channels of neutral and anionic Fe(m)F(n) (m = 1, 2; n = 1-7) clusters using density functional theory. Our results show that a maximum of six F atoms can be bound atomically to one Fe atom, and the Fe-Fe bonding is not preferred in Fe2F(n)(0/-) clusters. The computed electron affinities (EAs) indicate that FeF(n) with n ≥ 3 are superhalogens, while Fe2F(n) can be classified as superhalogens for n ≥ 5. To further understand their superhalogen characteristic, the natural population analysis charge distribution and the HOMOs of anionic clusters were also analyzed. When the extra negative charge and the content of HOMO are mainly located on F atoms, the clusters could be classified as superhalogens with EAs substantially surpass that of Cl. By calculating the binding energies per atom and the HOMO-LUMO gaps, FeF3, FeF4(-), Fe2F4, Fe2F5(-), and Fe2F7(-) clusters were found to have higher stabilities, corresponding to the Fe atoms that are attained at their favorite +2 and +3 oxidation states. Furthermore, we also predicted the most preferred fragmentation channel and product for all the ground state clusters. Even more striking is the fact that both neutral and anionic Fe(m)F(n) (m = 1, 2; n = 1-7) clusters carry large magnetic moments which mainly come from 3d orbital of iron atom.

Phase stability, mechanical properties, hardness, and possible reactive routing of chromium triboride from first-principle investigations.

The first-principles calculations are employed to provide a fundamental understanding of the structural features and relative stability, mechanical and electronic properties, and possible reactive route for chromium triboride. The predicted new phase of CrB3 belongs to the rhombohedral phase with R-3m symmetry and it transforms into a hexagonal phase with P-6m2 symmetry at 64 GPa. The mechanical and thermodynamic stabilities of CrB3 are verified by the calculated elastic constants and formation enthalpies. Also, the full phonon dispersion calculations confirm the dynamic stability of predicted CrB3. Considering the role of metallic contributions, the calculated hardness values from our semiempirical method for rhombohedral and hexagonal phases are 23.8 GPa and 22.1 GPa, respectively. In addition, the large shear moduli, Young's moduli, low Poisson's ratios, and small B∕G ratios indicate that they are potential hard materials. Relative enthalpy calculations with respect to possible constituents are also investigated to assess the prospects for phase formation and an attempt at high-pressure synthesis is suggested to obtain chromium triboride.

Probing the structural and electronic properties of small vanadium monoxide clusters.

The structural evolution and bonding of a series of early transition-metal oxide clusters, V(n)O(q) (n = 3-9, q = 0,-1), have been investigated with the aid of previous photoelectron spectroscopy (PES) and theoretical calculations. For each vanadium monoxide cluster, many low-lying isomers are generated using the Saunders "Kick" global minimum stochastic search method. Theoretical electron detachment energies (both vertical and adiabatic) were compared with the experimental measurements to verify the ground states of the vanadium monoxide clusters obtained from the DFT calculations. The results demonstrate that the combination of photoelectron spectroscopy experiments and DFT calculation is not only powerful for obtaining the electronic and atomic structures of size-selected clusters, but also valuable in resolving structurally and energetically close isomers. The second difference energies and adsorption energies as a function of the cluster size exhibit a pronounced even-odd alternation phenomenon. The adsorption energies of one O atom on the anionic (6.64 → 8.16 eV) and neutral (6.41 → 8.13 eV) host vanadium clusters are shown to be surprisingly high, suggesting strong capabilities to activate O by structural defects in vanadium oxides.

Geometries, stabilities, and electronic properties of Pt-group-doped gold clusters, their relationship to cluster size, and comparison with pure gold clusters.

A systematic study of bimetallic Au(n)M(2) (n = 1-6, M = Ni, Pd, and Pt) clusters is performed by using density functional theory at the B3LYP level. The geometric structures, relative stabilities, HOMO-LUMO gaps, natural charges and electronic magnetic moments of these clusters are investigated, and compared with pure gold clusters. The results indicate that the properties of Au(n)M(2) clusters for n = 1-3 diverge more from pure gold clusters, while those for n = 4-6 show good agreement with Au(n) clusters. The dissociation energies, the second-order difference of energies, and HOMO-LUMO energy gaps, exhibiting an odd-even alternation, indicate that the Au(4)M(2) clusters are the most stable structures for Au(n)M(2) (n = 1-6, M = Ni, Pd, and Pt) clusters. Moreover, we predict that the average atomic binding energies of these clusters should tend to a limit in the range 1.56-2.00 eV.

Equilibrium geometries, stabilities, and electronic properties of the bimetallic M2-doped Au(n) (M = Ag, Cu; n = 1-10) clusters: comparison with pure gold clusters.

The density functional method with relativistic effective core potential has been employed to investigate systematically the geometrical structures, relative stabilities, growth-pattern behaviors, and electronic properties of small bimetallic M(2)Au(n) (M = Ag, Cu; n = 1-10) and pure gold Au(n) (n ≤ 12) clusters. The optimized geometries reveal that M(2) substituted Au(n+2) clusters and one Au atom capped M(2)Au(n-1) structures are dominant growth patterns of the stable alloyed M(2)Au(n) clusters. The calculated averaged atomic binding energies, fragmentation energies, and the second-order difference of energies as a function of the cluster size exhibit a pronounced even-odd alternation phenomenon. The analytic results exhibit that the planar structure Ag(2)Au(4) and Cu(2)Au(2) isomers are the most stable geometries of Ag(2)Au(n) and Cu(2)Au(n) clusters, respectively. In addition, the HOMO-LUMO gaps, charge transfers, chemical hardnesses and polarizabilities have been analyzed and compared further.

Density functional study of structural and electronic properties of bimetallic copper-gold clusters: comparison with pure and doped gold clusters.

The geometrical structures, relative stabilities, and electronic properties of small bare gold clusters Au(n)(lambda) and bimetallic complexes of bare metal clusters with one copper atom Au(n-1)Cu(lambda) (charge lambda = 0, +1, -1; 2 < or = n < or = 9) have been systematically investigated by means of first-principles density functional calculations at the B3LYP level. The results show that the most stable isomers have a planar structure and resemble pure gold clusters in shape, and no three-dimensional isomers were obtained for neutral and anionic doped gold clusters. However, the geometries of Au(n-1)Cu(+) are found to undergo a structural change from two dimensional to three dimensional when the cluster contains 7 atoms. The calculated dissociation energy and second difference energy as a function of the cluster size exhibit a pronounced even-odd alternation phenomenon. Ionization potentials and electron detachment energies (both vertical and adiabatic) of Au(n)(lambda) and Au(n-1)Cu(lambda) clusters are discussed and compared with available experimental results. A good agreement between experimental and theoretical results suggests good prediction of the lowest energy structures for all clusters calculated in the present study.

Geometries, stabilities, and electronic properties of small anion Mg-doped gold clusters: a density functional theory study.

First-principle density functional theory is used for studying the anion gold clusters doped with magnesium atom. By performing geometry optimizations, the equilibrium geometries, relative stabilities, and electronic and magnetic properties of [Au(n)Mg]⁻ (n = 1-8) clusters have been investigated systematically in comparison with pure gold clusters. The results show that doping with a single Mg atom dramatically affects the geometries of the ground-state Au(n+1)⁻ clusters for n = 2-7. Here, the relative stabilities are investigated in terms of the calculated fragmentation energies, second-order difference of energies, and highest occupied−lowest unoccupied molecular orbital energy gaps, manifesting that the ground-state [Au(n)Mg]⁻ and Au(n+1)⁻ clusters with odd-number gold atoms have a higher relative stability. In particular, it should be noted that the [Au3Mg]⁻ cluster has the most enhanced chemical stability. The natural population analysis reveals that the charges in [Au(n)Mg]⁻ (n = 2-8) clusters transfer from the Mg atom to the Au frames. In addition, the total magnetic moments of [Au(n)Mg]⁻ clusters exhibit an odd-even oscillation as a function of cluster size, and the magnetic effects mainly come from the Au atoms.

Strain-induced structural phase transition, electric polarization and unusual electric properties in photovoltaic materials CsMI3 (M = Pb, Sn).

The structural phase transition, ferroelectric polarization, and electric properties have been investigated for photovoltaic films CsMI3 (M = Pb, Sn) epitaxially grown along (001) direction based on the density functional theory. The calculated results indicate that the phase diagrams of two epitaxial CsPbI3 and CsSnI3 films are almost identical, except critical transition strains varying slightly. The epitaxial tensile strains induce two ferroelectric phases Pmc21, and Pmn21, while the compressive strains drive two paraelectric phases P212121, P21212. The larger compressive strain enhances the ferroelectric instability in these two films, eventually rendering them another ferroelectric state Pc. Whether CsPbI3 or CsSnI3, the total polarization of Pmn21 phase comes from the main contribution of B-position cations (Pb or Sn), whereas, for Pmc21 phase, the main contributor is the I ion. Moreover, the epitaxial strain effects on antiferrodistortive vector, polarization and band gap of CsMI3 (M = Pb, Sn) are further discussed. Unusual electronic properties under epitaxial strains are also revealed and interpreted.

Structural, electronic, and magnetic properties of gold cluster anions doped with zinc: Au(n)Zn- (2 < or = n < or = 10).

The geometrical structures, relative stabilities, electronic, and magnetic properties of zinc-doped gold cluster anions Au(n)Zn(-) (2 < or = n < or = 10) have been systematically investigated by means of first-principles density functional calculations at the B3LYP level. The results show that the most stable isomers have a planar structure and resemble pure gold cluster anions in shape, and no 3D isomers were obtained. Calculated dissociation energy, the second difference energy, and the highest occupied-lowest unoccupied molecular orbital gaps as a function of the cluster size exhibit a pronounced even-odd alternation phenomenon. The first vertical detachment energy of the anion clusters were calculated and compared with available experimental results. A good agreement between experimental and theoretical results suggests good prediction of the lowest-energy structures for all clusters calculated in the present study. Furthermore, both the local and total magnetic moments display a pronounced odd-even oscillation with the number of gold atoms.

Theoretical study on local defect structure of (FeO4)(5-) clusters in YGG and LGG crystals.

The optical spectrum and the local defect structure of tetrahedral (FeO4)(5-) clusters are investigated in yttrium gallium garnet (YGG) and lutetium gallium garnet (LGG) crystals by means of the 252x252 complete energy matrices for d(5) configuration ions in tetragonal ligand field. The results show that the local defect structures around tetrahedral Fe3+ centers display an expansion effect. Simultaneously, the (FeO4(5-) clusters in the two different crystals have very similar local structures, which are close to those in YIG garnet crystal. Finally the relative curves of the zero-field splitting energies DeltaE1 and DeltaE2 in the ground state (6)A(1) varying with the structure parameter theta are plotted.

Hydrostatic pressure induced structural phase transition and mechanical properties of fluoroperovskite.

We perform the first-principles calculations combined with the particle swarm optimization algorithm to investigate the high-pressure phase diagrams of Na[Formula: see text]F3 ([Formula: see text] = Mn, Ni, Zn). Two reconstructive phase transitions are predicted from Pv-[Formula: see text] to pPv-[Formula: see text] at about 9 GPa and pPv-[Formula: see text] to ppPv-[Formula: see text] at around 26 GPa for NaZnF3. That is not the case for NaMnF3-a direct transition (reconstructive transition in nature but with the same Pnma space group) from Pv-[Formula: see text] to ppPv-[Formula: see text] phase around 12 GPa. Strikingly, our simulated results manifest that a disproportionation phase of NaZnF3 post-perovskite is uncovered along the way, which provides a successful explanation for the observed results in experiment. Additionally, the mechanical and thermal properties, especially the dynamical property, of the four NaZnF3 phases have also been studied. Here, we reveal the obvious softening of [Formula: see text]-wave velocity and bulk sound speed in pPv-[Formula: see text]-to-ppPv-[Formula: see text] transition, which may result in the discontinuity of seismic waves propagation through the Earth's interior.

Characterization of the high-pressure structural transition and thermodynamic properties in sodium chloride: a computational investigation on the basis of the density functional theory.

Using first-principles calculations, the elastic constants, the thermodynamic properties, and the structural phase transition between the B1 (rocksalt) and the B2 (cesium chloride) phases of NaCl are investigated by means of the pseudopotential plane-waves method. The calculations are performed within the generalized gradient approximation to density functional theory with the Perdew-Burke-Ernzerhof exchange-correlation functional. On the basis of the third-order Birch-Murnaghan equation of states, the transition pressure Pt between the B1 phase and the B2 phase of NaCl is determined. The calculated values are generally speaking in good agreement with experiments and with similar theoretical calculations. From the theoretical calculations, the shear modulus, Young's modulus, rigidity modulus, and Poisson's ratio of NaCl are derived. According to the quasi-harmonic Debye model, we estimated the Debye temperature of NaCl from the average sound velocity. Moreover, the pressure derivatives of elastic constants, partial differentialC11/partial differentialP, partial differentialC12/partial differentialP, partial differentialC44/partial differentialP, partial differentialS11/partial differential P, partial differentialS12/partial differentialP, and partial differentialS44/partial differentialP, for NaCl crystal are investigated for the first time. This is a quantitative theoretical prediction of the elastic and thermodynamic properties of NaCl, and it still awaits experimental confirmation.