No Arabic abstract
Equilibrium polyethylene crystal structure, cohesive energy, and elastic constants are calculated by density-functional theory applied with a recently proposed density functional (vdW-DF) for general geometries [Phys. Rev. Lett. 92, 246401 (2004)] and with a pseudopotential-planewave scheme. The vdW-DF with its account for the long-ranged van der Waals interactions gives not only a stabilized crystal structure but also values of the calculated lattice parameters and elastic constants in quite good agreement with experimental data, giving promise for successful application to a wider range of polymers.
We report a study of the anisotropic exchange interactions in bulk CrO_2 calculated from first principles within density functional theory. We determine the exchange coupling energies, using both the experimental lattice parameters and those obtained within DFT, within a modified Heisenberg model Hamiltonian in two ways. We employ a supercell method in which certain spins within a cell are rotated and the energy dependence is calculated and a spin-spiral method that modifies the periodic boundary conditions of the problem to allow for an overall rotation of the spins between unit cells. Using the results from each of these methods, we calculate the spin-wave stiffness constant D from the exchange energies using the magnon dispersion relation. We employ a Monte Carlo method to determine the DFT-predicted Curie temperature from these calculated energies and compare with accepted values. Finally, we offer an evaluation of the accuracy of the DFT-based methods and suggest implications of the competing ferro- and antiferromagnetic interactions.
We present the implementation of GAtor, a massively parallel, first principles genetic algorithm (GA) for molecular crystal structure prediction. GAtor is written in Python and currently interfaces with the FHI-aims code to perform local optimizations and energy evaluations using dispersion-inclusive density functional theory (DFT). GAtor offers a variety of fitness evaluation, selection, crossover, and mutation schemes. Breeding operators designed specifically for molecular crystals provide a balance between exploration and exploitation. Evolutionary niching is implemented in GAtor by using machine learning to cluster the dynamically updated population by structural similarity and then employing a cluster-based fitness function. Evolutionary niching promotes uniform sampling of the potential energy surface by evolving several sub-populations, which helps overcome initial pool biases and selection biases (genetic drift). The various settings offered by GAtor increase the likelihood of locating numerous low-energy minima, including those located in disconnected, hard to reach regions of the potential energy landscape. The best structures generated are re-relaxed and re-ranked using a hierarchy of increasingly accurate DFT functionals and dispersion methods. GAtor is applied to a chemically diverse set of four past blind test targets, characterized by different types of intermolecular interactions. The experimentally observed structures and other low-energy structures are found for all four targets. In particular, for Target II, 5-cyano-3-hydroxythiophene, the top ranked putative crystal structure is a $Z^prime$=2 structure with P$bar{1}$ symmetry and a scaffold packing motif, which has not been reported previously.
Finite-temperature Kohn--Sham density-functional theory (KS-DFT) is a widely-used method in warm dense matter (WDM) simulations and diagnostics. Unfortunately, full KS-DFT-molecular dynamics models scale unfavourably with temperature and there remains uncertainty regarding the performance of existing approximate exchange-correlation (XC) functionals under WDM conditions. Of particular concern is the expected explicit dependence of the XC functional on temperature, which is absent from most approximations. Average-atom (AA) models, which significantly reduce the computational cost of KS-DFT calculations, have therefore become an integral part of WDM modelling. In this paper, we present a derivation of a first-principles AA model from the fully-interacting many-body Hamiltonian, carefully analysing the assumptions made and terms neglected in this reduction. We explore the impact of different choices within this model -- such as boundary conditions and XC functionals -- on common properties in WDM, for example equation-of-state data. Furthermore, drawing upon insights from ground-state KS-DFT, we speculate on likely sources of error in KS-AA models and possible strategies for mitigating against such errors.
In this paper the relationship between the density functional theory of freezing and phase field modeling is examined. More specifically a connection is made between the correlation functions that enter density functional theory and the free energy functionals used in phase field crystal modeling and standard models of binary alloys (i.e., regular solution model). To demonstrate the properties of the phase field crystal formalism a simple model of binary alloy crystallization is derived and shown to simultaneously model solidification, phase segregation, grain growth, elastic and plastic deformations in anisotropic systems with multiple crystal orientations on diffusive time scales.
Chemisorbed molecules at a fuel cell electrode are a very sensitive probe of the surrounding electrochemical environment, and one that can be accurately monitored with different spectroscopic techniques. We develop a comprehensive electrochemical model to study molecular chemisorption at either constant charge or fixed applied voltage, and calculate from first principles the voltage dependence of vibrational frequencies -- the vibrational Stark effect -- for CO adsorbed on close-packed platinum electrodes. The predicted vibrational Stark slopes are found to be in very good agreement with experimental electrochemical spectroscopy data, thereby resolving previous controversies in the quantitative interpretation of in-situ experiments and elucidating the relation between canonical and grand-canonicaldescriptions of vibrational surface phenomena.