Source author record

Igor A. Abrikosov

Igor A. Abrikosov appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2024arXiv

Strong electron-phonon coupling and phonon-induced superconductivity in tetragonal C$_3$N$_4$ with hole doping

C$_3$N$_4$ is a recently discovered phase of carbon nitrides with the tetragonal crystal structure [D.Laniel $\textit{et al.}$, Adv. Mater. 2023, 2308030] that is stable at ambient conditions. C$_3$N$_4$ is a semiconductor exhibiting flat-band anomalies in the valence band, suggesting the emergence of many-body instabilities upon hole doping. Here, using state-of-the-art first-principles calculations we show that hole-doped C$_3$N$_4$ reveals strong electron-phonon coupling, leading to the formation of a gapped superconducting state. The phase transition temperatures turn out to be strongly dependent on the hole concentration. We propose that holes could be injected into C$_3$N$_4$ via boron doping which induces, according to our results, a rigid shift of the Fermi energy without significant modification of the electronic structure. Based on the electron-phonon coupling and Coulomb pseudopotential calculated from first principles, we conclude that the boron concentration of 6 atoms per nm$^3$ would be required to reach the critical temperature of $\sim$36 K at ambient pressure.

preprint2022arXiv

Electronic and optical properties of crystalline nitrogen versus black phosphorus: A comparative first-principles study

Crystalline black nitrogen (BN) is an allotrope of nitrogen with the black phosphorus (BP) structure recently synthesized at high pressure by two independent research groups [Ji et al., Sci. Adv. 6, eaba9206 (2020); Laniel et al., Phys. Rev. Lett. 124, 216001 (2020)]. Here, we present a systematic study of the electronic and optical properties of BN focusing on its comparison with BP. To this end, we use the state-of-the-art quasiparticle self-consistent $GW$ approach with vertex corrections in both the electronic and optical channels. Despite many similarities, the properties of BN are found to be considerably different. Unlike BP, BN exhibits a larger optical gap (2.5 vs 0.26 eV), making BN transparent in the visible spectral region with a highly anisotropic optical response. This difference can be primarily attributed to a considerably reduced dielectric screening in BN, leading to enhancement of the effective Coulomb interaction. Despite relatively strong Coulomb interaction, exciton formation is largely suppressed in both materials. Our analysis of the elastic properties shows exceptionally high stiffness of BN, comparable to that of diamond.

preprint2022arXiv

Isotope purification induced reduction of spin relaxation and spin coherence times in semiconductors

Paramagnetic defects and nuclear spins are often the major sources of decoherence and spin relaxation in solid-state qubits realized by optically addressable point defect spins in semiconductors. It is commonly accepted that a high degree of depletion of nuclear spins can enhance the coherence time by reducing magnetic noise. Here we show that the isotope purification beyond a certain optimal level becomes contra-productive, when both electron and nuclear spins are present in the vicinity of the qubits. Using state-of-the-art numerical tools and considering the silicon vacancy qubit in various spin environments, we demonstrate that the coupling to spin-1/2 point defects in the lattice can be significantly enhanced by isotope purification. The enhanced coupling shortens the spin relaxation time that in turn may limit the the coherence time of spin qubits. Our results can be straightforwardly generalized to triplet point defect qubits, such as the NV center in diamond and the divacancy in SiC.

preprint2022arXiv

Revealing Phosphorus Nitrides up to the Megabar Regime: Synthesis of alpha-P3N5, delta-P3N5 and PN2

Non-metal nitrides are an exciting field of chemistry, featuring a significant number of compounds that can possess outstanding material properties. This characteristic relies on maximizing the number of strong covalent bonds, with crosslinked XN6 octahedra frameworks being particularly intriguing. In this study, the phosphorus-nitrogen system was studied up to 137 GPa in laser-heated diamond anvil cells and three previously unobserved phases were synthesized and characterized by single-crystal X-ray diffraction, Raman spectroscopy measurements and density functional theory calculations. Delta-P3N5 and PN2 were found to form at 72 and 134 GPa, respectively, and both feature dense 3D networks of the so far elusive PN6 units. The two are ultra-incompressible, having a bulk modulus of K0 = 322 GPa for delta-P3N5 and K0 = 339 GPa for PN2. Upon decompression below 7 GPa, delta-P3N5 undergoes a transformation into a novel alpha'-P3N5 solid, stable at ambient conditions, that has a unique structure type based on PN4 tetrahedra. The formation of alpha'-P3N5 underlines that a phase space otherwise inaccessible can be explored through high-pressure formed phases.

preprint2020arXiv

Photoluminescence at the ground state level anticrossing of the nitrogen-vacancy center in diamond

The nitrogen-vacancy center (NV center) in diamond at magnetic fields corresponding to the ground state level anticrossing (GSLAC) region gives rise to rich photoluminescence (PL) signals due to the vanishing energy gap between the electron spin states, which enables to have an effect on the NV center's luminescence for a broad variety of environmental couplings. In this article we report on the GSLAC photoluminescence signature of NV ensembles in different spin environments at various external fields. We investigate the effects of transverse electric and magnetic fields, P1 centers, NV centers, and the $^{13}$C nuclear spins, each of which gives rise to a unique PL signature at the GSLAC. The comprehensive analysis of the couplings and related optical signal at the GSLAC provides a solid ground for advancing various microwave-free applications at the GSLAC, including but not limited to magnetometry, spectroscopy, dynamic nuclear polarization (DNP), and nuclear magnetic resonance (NMR) detection. We demonstrate that not only the most abundant $^{14}$NV center but the $^{15}$NV can also be utilized in such applications and that nuclear spins coupled to P1 centers can be polarized directly by the NV center at the GSLAC, through a giant effective nuclear $g$-factor arising from the NV center-P1 center-nuclear spin coupling. We report on new alternative for measuring defect concentration in the vicinity of NV centers and on the optical signatures of interacting, mutually aligned NV centers.

preprint2019arXiv

Anharmonicity and Ultra-Low Thermal Conductivity in Lead-Free Halide Double Perovskites

The lead-free halide double perovskite class of materials offers a promising venue for resolving issues related to toxicity of Pb and long-term stability of the lead-containing halide perovskites. We present a first-principles study of the lattice vibrations in Cs$_2$AgBiBr$_6$ , the prototypical compound in this class, and show that the lattice dynamics of Cs$_2$AgBiBr$_6$ is highly anharmonic, largely in regards to tilting of AgBr$_6$ and BiBr$_6$ octahedra. Using an energy and temperature dependent phonon spectral function, we then show how the experimentally observed cubic-to-tetragonal phase transformation is caused by the collapse of a soft phonon branch. We finally reveal that the softness and anharmonicity of Cs$_2$AgBiBr$_6$ yield an ultra-low thermal conductivity, unexpected of high symmetry cubic structures.

preprint2019arXiv

High-pressure synthesis of ultraincompressible hard rhenium nitride pernitride Re$_{2}$(N$_{2}$)N$_{2}$ stable at ambient conditions

Here we report the synthesis of metallic, ultraincompressible (bulk modulus $K_{0}$ = 428(10) GPa) and very hard (nanoindentation hardness 36.7(8) GPa) rhenium (V) nitride pernitride Re$_{2}$(N$_{2}$)N$_{2}$. While the empirical chemical formula of the compound, ReN$_{2}$, is the same as for other known transition metals pernitrides, e.g. IrN$_{2}$, PtN$_{2}$, PdN$_{2}$ and OsN$_{2}$, its crystal chemistry is unique. The known pernitrides of transition metals consist of a metal in the oxidation state +IV and pernitride anions N$_{2}^{4-}$. ReN$_{2}$ contains both pernitride N$_{2}^{4-}$ and discrete N$^{3-}$ anions, which explains its exceptional properties. Moreover, in the original experimental synthesis of Re$_{2}$(N$_{2}$)N$_{2}$ performed in a laser-heated diamond anvil cell via a direct reaction between rhenium and nitrogen at pressures from 40 to 90 GPa we observed that the material was recoverable at ambient conditions. Consequently, we developed a route to scale up its synthesis through a reaction between rhenium and ammonium azide, NH$_{4}$N$_{3}$, in a large-volume press at 33 GPa. Our work resulted not only in a discovery of a novel material with unusual crystal chemistry and a set of properties attractive for potential applications, but also demonstrated a feasibility of surmounting conceptions common in material sciences.

preprint2016arXiv

Efficient and accurate determination of lattice-vacancy diffusion coefficients via non equilibrium ab initio molecular dynamics

We revisit the color-diffusion algorithm [P. C. Aeberhard et al., Phys. Rev. Lett. 108, 095901 (2012)] in nonequilibrium ab initio molecular dynamics (NE-AIMD), and propose a simple efficient approach for the estimation of monovacancy jump rates in crystalline solids at temperatures well below melting. Color-diffusion applied to monovacancy migration entails that one lattice atom (colored-atom) is accelerated toward the neighboring defect-site by an external constant force F. Considering bcc molybdenum between 1000 and 2800 K as a model system, NE-AIMD results show that the colored-atom jump rate k_{NE} increases exponentially with the force intensity F, up to F values far beyond the linear-fitting regime employed previously. Using a simple model, we derive an analytical expression which reproduces the observed k_{NE}(F) dependence on F. Equilibrium rates extrapolated by NE-AIMD results are in excellent agreement with those of unconstrained dynamics. The gain in computational efficiency achieved with our approach increases rapidly with decreasing temperatures, and reaches a factor of four orders of magnitude at the lowest temperature considered in the present study.

preprint2016arXiv

Finite temperature elastic constants of paramagnetic materials within the disordered local moment picture from ab initio molecular dynamics calculations

We present a theoretical scheme to calculate the elastic constants of magnetic materials in the high-temperature paramagnetic state. Our approach is based on a combination of disordered local moments picture and ab initio molecular dynamics (DLM-MD). Moreover, we investigate a possibility to enhance the efficiency of the simulations of elastic properties using recently introduced method: symmetry imposed force constant temperature dependent effective potential (SIFC-TDEP). We have chosen cubic paramagnetic CrN as a model system. This is done due to its technological importance and its demonstrated strong coupling between magnetic and lattice degrees of freedom. We have studied the temperature dependent single-crystal and polycrystalline elastic constants of paramagentic CrN up to 1200 K. The obtained results at T= 300 K agree well with the experimental values of polycrystalline elastic constants as well as Poisson ratio at room temperature. We observe that the Young's modulus is strongly dependent on temperature, decreasing by ~14% from T=300 K to 1200 K. In addition we have studied the elastic anisotropy of CrN as a function of temperature and we observe that CrN becomes substantially more isotropic as the temperature increases. We demonstrate that the use of Birch law may lead to substantial errors for calculations of temperature induced changes of elastic moduli. The proposed methodology can be used for accurate predictions of mechanical properties of magnetic materials at temperatures above their magnetic order-disorder phase transition.

preprint2016arXiv

High fidelity bi-directional nuclear qubit initialization in SiC

Dynamic nuclear polarization (DNP) is an attractive method for initializing nuclear spins that are strongly coupled to optically active electron spins because it functions at room temperature and does not require strong magnetic fields. In this Letter, we demonstrate that DNP, with near-unity polarization efficiency, can be generally realized in weakly coupled hybrid registers, and furthermore that the nuclear spin polarization can be completely reversed with only sub-Gauss magnetic field variations. This mechanism offers new avenues for DNP-based sensors and radio-frequency free control of nuclear qubits.

preprint2016arXiv

Structural and electronic properties of Li intercalated graphene on SiC(0001)

We investigate the structural and electronic properties of Li-intercalated monolayer graphene on SiC(0001) using combined angle-resolved photoemission spectroscopy and first-principles density functional theory. Li intercalates at room temperature both at the interface between the buffer layer and SiC and between the two carbon layers. The graphene is strongly $n$-doped due to charge transfer from the Li atoms and two $π$-bands are visible at the $\bar{K}$-point. After heating the sample to 300$^\circ$C, these $π$-bands become sharp and have a distinctly different dispersion to that of Bernal-stacked bilayer graphene. We suggest that the Li atoms intercalate between the two carbon layers with an ordered structure, similar to that of bulk LiC$_6$. An AA-stacking of these two layers becomes energetically favourable. The $π$-bands around the $\bar{K}$-point closely resemble the calculated band structure of a C$_6$LiC$_6$ system, where the intercalated Li atoms impose a super-potential on the graphene electronic structure that opens pseudo-gaps at the Dirac points of the two $π$-cones.

preprint2015arXiv

Anharmonicity changes the solid solubility of an alloy at high temperatures

We have developed a method to accurately and efficiently determine the vibrational free energy as a function of temperature and volume for substitutional alloys from first principles. Taking Ti$_{1-x}$Al$_x$N alloy as a model system, we calculate the isostructural phase diagram by finding the global minimum of the free energy, corresponding to the true equilibrium state of the system. We demonstrate that the anharmonic contribution and temperature dependence of the mixing enthalpy have a decisive impact on the calculated phase diagram of a Ti$_{1-x}$Al$_x$N alloy, lowering the maximum temperature for the miscibility gap from 6560 K to 2860 K. Our local chemical composition measurements on thermally aged Ti$_{0.5}$Al$_{0.5}$N alloys agree with the calculated phase diagram.

preprint2015arXiv

Dirac points with giant spin-orbit splitting in the electronic structure of two-dimensional transition-metal carbides

Two-dimensional (2D) materials, especially their most prominent member, graphene, have greatly influenced many scientific areas. Moreover, they have become a base for investigating the relativistic properties of condensed matter within the emerging field of Dirac physics. This has ignited an intense search for new materials where charge carriers behave as massless or massive Dirac fermions. Here, we theoretically show the existence of Dirac electrons in a series of 2D transition-metal carbides, known as MXenes. They possess twelve conical crossings in the 1st Brillouin zone with giant spin-orbit splitting. Our findings indicate that the 2D band structure of MXenes is protected against external perturbations and preserved even in multilayer phases. These results, together with the broad possibilities to engineer the properties of these materials phases, make Dirac MXenes a potential candidate for studying and developing novel Dirac-physics-based technologies.

preprint2015arXiv

Pressure induced crossing of the core-levels in 5d metals

Pressure induced interaction between core electrons, the core level crossing (CLC) transition has been observed in hcp Os at P~400 GPa [L. Dubrovinksy, et al., Nature 525, 226-229 (2015)]. In this work, we carry out a systematic study of the influence of pressure on the electronic structure in all metals of the 5d series (Hf,Ta,W,Re,Os,Ir,Pt,Au) using first-principles electronic structure calculations. We have found that CLC is a general effect for this series of metals. While in Pt it occurs at ~1500 GPa, at a pressure substantially higher than in Os, in Ir it occurs already at 80 GPa. Moreover, we predict that in Re the CLC transition may appear at ambient pressure. We analyze the shifts of the CLC transition pressure across the series within the Thomas-Fermi model, and show that the effect has many common features to the atomic collapse in the rare-earth elements.

preprint2015arXiv

The morphology transition mechanism from icosahedral to decahedral phase during growth of nanoclusters

The morphology transition from the thermodynamically favorable to the unfavorable phase during growth of free-standing copper nanoclusters is studied by molecular dynamics simulations. We give a detailed description of the kinetics and thermodynamics of the process. A universal mechanism of a solid-solid transition, from icosahedral to decahedral morphology in nanoclusters, is proposed. We show that a formation of distorted NC during the growth process with islands of incoming atoms localized in certain parts of the grown particle may shift the energy balance between Ih and Dh phases in favour of the latter leading to the morphology transition deep within the thermodynamic stability field of the former. The role of diffusion in the morphology transition is revealed. In particular, it is shown that fast diffusion should suppress the morphology transition and favour homogeneous growth of the nanoclusters.

preprint2015arXiv

Theoretical model of the dynamic spin polarization of nuclei coupled to paramagnetic point defects in diamond and silicon carbide

Dynamic nuclear spin polarization (DNP) mediated by paramagnetic point defects in semiconductors is a key resource for both initializing nuclear quantum memories and producing nuclear hyperpolarization. DNP is therefore an important process in the field of quantum-information processing, sensitivity-enhanced nuclear magnetic resonance, and nuclear-spin-based spintronics. DNP based on optical pumping of point defects has been demonstrated by using the electron spin of nitrogen-vacancy (NV) center in diamond, and more recently, by using divacancy and related defect spins in hexagonal silicon carbide (SiC). Here, we describe a general model for these optical DNP processes that allows the effects of many microscopic processes to be integrated. Applying this theory, we gain a deeper insight into dynamic nuclear spin polarization and the physics of diamond and SiC defects. Our results are in good agreement with experimental observations and provide a detailed and unified understanding. In particular, our findings show that the defects' electron spin coherence times and excited state lifetimes are crucial factors in the entire DNP process.

preprint2014arXiv

Molecular dynamics simulation of growth of Cu nanoclusters from Cu-ions in a plasma

A recently developed method of nanoclusters growth in a pulsed plasma is studied by means of molecular dynamics. A model that allows one to consider high-energy charged particles in classical molecular dynamics is suggested, and applied for studies of single impact events in nanoclusters growth. In particular, we provide a comparative analysis of the well-studied inert gas aggregation method and the growth from ions in a plasma. The importance to consider of the angular distribution of incoming ions in the simulations of the nanocluster growth is underlined. A detailed study of the energy transfer from the incoming ions to a nanocluster, as well as the diffusion of incoming ions on the cluster surface is carried out. Our results are important for understanding and control of the nanocluster growth process.

preprint2014arXiv

The role of N defects in paramagnetic CrN at finite temperatures from first-principles

Simulations of defects in paramagnetic materials at high temperature constitute a formidable challenge to solid state theory due to the interaction of magnetic disorder, vibrations, and structural relaxations. CrN is a material where these effects are particularly large due to a strong magneto-lattice coupling and a tendency for deviations from the nominal 1:1 stoichiometry. In this work we present a first-principles study of nitrogen vacancies and nitrogen interstitals in CrN at elevated temperature. We report on formation energetics, the geometry of interstital nitrogen dimers, and impact on the electronic structure caused by the defects. We find a vacancy formation energy of 2.28 eV with a small effect of temperature, a formation energy for N interstitial in the form of a $\big<111\big>$ oriented split-bond of 3.77 eV with an increase to 3.97 at 1000 K. Vacancies are found to add three electrons while split bond interstitial adds one electron to the conduction band. The band gap of defect-free CrN is smeared out due to vibrations, though it is difficult to draw conclusion about the exact temperature at which the band gap closes from our calculations. However, it is clear that at 900 K there is a non-zero density of electronic states at the Fermi level. At 300 K our results indicate a boarder case were the band-gap is about to close.

preprint2014arXiv

Theoretical unification of hybrid-DFT and DFT+U methods for the treatment of localized orbitals

We formulate the on-site occupation dependent exchange correlation energy and effective potential of hybrid functionals for localized states and connect them to the on-site correction term of the DFT+U method. Our derivation provides a theoretical justification for adding a DFT+U-like onsite potential in hybrid DFT calculations to resolve issues caused by overscreening of localized states. The resulting scheme, hybrid- DFT+Vw, is tested for chromium impurity in wurtzite AlN and vanadium impurity in 4H-SiC, which are paradigm examples of systems with different degree of localization between host and impurity orbitals.

preprint2014arXiv

Vibrational free energy and phase stability of paramagnetic and antiferromagnetic CrN from ab-initio molecular dynamics

We present a theoretical first-principles method to calculate the free energy of a magnetic system in its high-temperature paramagnetic phase, including vibrational, electronic, and magnetic contributions. The method for calculating free energies is based on ab-initio molecular dynamics and combines a treatment of disordered magnetism using disordered local moments molecular dynamics (DLM-MD) with the temperature dependent effective potential (TDEP) method to obtain the vibrational contribution to the free energy. We illustrate the applicability of the method by obtaining the anharmonic free energy for the paramagnetic cubic and the antiferromagnetic orthorhombic phases of chromium nitride. The influence of lattice dynamics on the transition between the two phases is demonstrated by constructing the temperature-pressure phase diagram.

preprint2013arXiv

Electrically and mechanically tunable electron spins in silicon carbide color centers

The electron spins of semiconductor defects can have complex interactions with their host, particularly in polar materials like SiC where electrical and mechanical variables are intertwined. By combining pulsed spin resonance with ab-initio simulations, we show that spin-spin interactions within SiC neutral divacancies give rise to spin states with an enhanced Stark effect, sub-10**-6 strain sensitivity, and highly spin-dependent photoluminescence with intensity contrasts of 15-36%. These results establish SiC color centers as compelling systems for sensing nanoscale fields.

preprint2013arXiv

Temperature dependent effective potential method for accurate free energy calculations of solids

We have developed a thorough and accurate method of determining anharmonic free energies, the temperature dependent effective potential technique (TDEP). It is based on \emph{ab initio} molecular dynamics followed by a mapping onto a model Hamiltonian that describes the lattice dynamics. The formalism and the numerical aspects of the technique are described in details. A number of practical examples are given, and results are presented, which confirm the usefulness of TDEP within \emph{ab initio} and classical molecular dynamics frameworks. In particular, we examine from first-principles the behavior of force constants upon the dynamical stabilization of body centered phase of Zr, and show that they become more localized. We also calculate phase diagram for $^4$He modeled with the Aziz \emph{et al.} potential and obtain results which are in favorable agreement both with respect to experiment and established techniques.

preprint2013arXiv

Temperature dependent effective third order interatomic force constants from first principles

The temperature dependent effective potential (TDEP) method is generalized beyond pair interactions. The second and third order force constants are determined consistently from ab initio molecular dynamics simulations at finite temperature. The reliability of the approach is demonstrated by calculations of the Mode Grüneisen parameters for Si. We show that the extension of TDEP to higher order allows for an efficient calculation of the phonon life time, in Si as well as in $ε$-FeSi, a system that exhibits anomalous softening with temperature.

preprint2011arXiv

Electronic structure investigation of the cubic inverse perovskite Sc3AlN

The electronic structure and chemical bonding of the recently discovered inverse perovskite Sc3AlN, in comparison to ScN and Sc metal have been investigated by bulk-sensitive soft x-ray emission spectroscopy. The measured Sc L, N K, Al L1, and Al L2,3 emission spectra are compared with calculated spectra using first principle density-functional theory including dipole transition matrix elements. The main Sc 3d - N 2p and Sc 3d - Al 3p chemical bond regions are identified at -4 eV and -1.4 eV below the Fermi level, respectively. A strongly modified spectral shape of 3s states in the Al L2,3 emission from Sc3AlN in comparison to pure Al metal is found, which reflects the Sc 3d - Al 3p hybridization observed in the Al L1 emission. The differences between the electronic structure of Sc3AlN, ScN, and Sc metal are discussed in relation to the change of the conductivity and elastic properties.

preprint2011arXiv

Equation of state of paramagnetic CrN from ab initio molecular dynamics

Equation of state for chromium nitride has been debated in the literature in connection with a proposed collapse of its bulk modulus following the pressure induced transition from the paramagnetic cubic phase to the antiferromagnetic orthorhombic phase [F. Rivadulla et al., Nat Mater 8, 974 (2009); B. Alling et al., Nat Mater 9, 283 (2010)]. Experimentally the measurements are complicated due to the low transition pressure, while theoretically the simulation of magnetic disorder represent a major challenge. Here a first-principles method is suggested for the calculation of thermodynamic properties of magnetic materials in their high temperature paramagnetic phase. It is based on ab-initio molecular dynamics and simultaneous redistributions of the disordered but finite local magnetic moments. We apply this disordered local moments molecular dynamics method to the case of CrN and simulate its equation of state. In particular the debated bulk modulus is calculated in the paramagnetic cubic phase and is shown to be very similar to that of the antiferromagnetic orthorhombic CrN phase for all considered temperatures.

preprint2010arXiv

Origin of the anomalous piezoelectric response in wurtzite Sc$_x$Al$_{1-x}$N alloys

The origin of the anomalous, 400% increase of the piezoelectric coefficient in Sc$_x$Al$_{1-x}$N alloys is revealed. Quantum mechanical calculations show that the effect is intrinsic. It comes from a strong change in the response of the internal atomic coordinates to strain and pronounced softening of C$_{33}$ elastic constant. The underlying mechanism is the flattening of the energy landscape due to a competition between the parent wurtzite and the so far experimentally unknown hexagonal phases of the alloy. Our observation provides a route for the design of materials with high piezoelectric response.

preprint2010arXiv

Significant elastic anisotropy in Ti$_{1-x}$Al$_x$N alloys

Strong compositional-dependent elastic properties have been observed theoretically and experimentally in Ti$_{1-x}$Al$_x$ N alloys. The elastic constant, C$_{11}$, changes by more than 50% depending on the Al-content. Increasing the Al-content weakens the average bond strength in the local octahedral arrangements resulting in a more compliant material. On the other hand, it enhances the directional (covalent) nature of the nearest neighbor bonds that results in greater elastic anisotropy and higher sound velocities. The strong dependence of the elastic properties on the Al-content offers new insight into the detailed understanding of the spinodal decomposition and age hardening in Ti$_{1-x}$Al$_x$N alloys.

preprint2010arXiv

Wavefunction extended Lagrangian Born-Oppenheimer molecular dynamics

Extended Lagrangian Born-Oppenheimer molecular dynamics [Niklasson, Phys. Rev. Lett. 100 123004 (2008)] has been generalized to the propagation of the electronic wavefunctions. The technique allows highly efficient first principles molecular dynamics simulations using plane wave pseudopotential electronic structure methods that are stable and energy conserving also under incomplete and approximate self-consistency convergence. An implementation of the method within the planewave basis set is presented and the accuracy and efficiency is demonstrated both for semi-conductor and metallic materials.