Algorithms for the symmetry-adapted energy minimisation of solids using analytical first and second derivatives have been devised and implemented in a new computer program GULP. These new methods are found to lead to an improvement in computational efficiency of up to an order of magnitude over the standard algorithm, which takes no account of symmetry, the largest improvement being obtained from the use of symmetry in the generation of the hessian. Accelerated convergence techniques for the dispersion energy are found to be beneficial in improving the precision at little extra computational cost, particularly when a one centre decomposition is possible or the Ewald sum weighting towards real-space is increased.

Crystallisation is at the heart of various scientific disciplines, but still the understanding of the molecular mechanisms underlying phase separation and the formation of the first solid particles in aqueous solution is rather limited. In this review, classical nucleation theory, as well as established concepts of spinodal decomposition and liquid–liquid demixing, is introduced together with a description of the recently proposed pre-nucleation cluster pathway. The features of pre-nucleation clusters are presented and discussed in relation to recent modifications of the classical and established models for phase separation, together with a review of experimental work and computer simulations on the characteristics of pre-nucleation clusters of calcium phosphate, calcium carbonate, iron(oxy)(hydr)oxide, silica, and also amino acids as an example of small organic molecules. The role of pre-nucleation clusters as solute precursors in the emergence of a new phase is summarized, and the link between the chemical speciation of homogeneous solutions and the process of phase separation via pre-nucleation clusters is highlighted.

Recent developments in and around the SIESTA method of first-principles simulation of condensed matter are described and reviewed, with emphasis on (i) the applicability of the method for large and varied systems, (ii) efficient basis sets for the standards of accuracy of density-functional methods, (iii) new implementations, and (iv) extensions beyond ground-state calculations.

A review of the present status, recent enhancements, and applicability of the Siesta program is presented. Since its debut in the mid-1990s, Siesta’s flexibility, efficiency, and free distribution have given advanced materials simulation capabilities to many groups worldwide. The core methodological scheme of Siesta combines finite-support pseudo-atomic orbitals as basis sets, norm-conserving pseudopotentials, and a real-space grid for the representation of charge density and potentials and the computation of their associated matrix elements. Here, we describe the more recent implementations on top of that core scheme, which include full spin–orbit interaction, non-repeated and multiple-contact ballistic electron transport, density functional theory (DFT)+U and hybrid functionals, time-dependent DFT, novel reduced-scaling solvers, density-functional perturbation theory, efficient van der Waals non-local density functionals, and enhanced molecular-dynamics options. In addition, a substantial effort has been made in enhancing interoperability and interfacing with other codes and utilities, such as wannier90 and the second-principles modeling it can be used for, an AiiDA plugin for workflow automatization, interface to Lua for steering Siesta runs, and various post-processing utilities. Siesta has also been engaged in the Electronic Structure Library effort from its inception, which has allowed the sharing of various low-level libraries, as well as data standards and support for them, particularly the PSeudopotential Markup Language definition and library for transferable pseudopotentials, and the interface to the ELectronic Structure Infrastructure library of solvers. Code sharing is made easier by the new open-source licensing model of the program. This review also presents examples of application of the capabilities of the code, as well as a view of on-going and future developments.

Calcium carbonate is an abundant substance that can be created in several mineral forms by the reaction of dissolved carbon dioxide in water with calcium ions. Through biomineralization, organisms can harness and control this process to form various functional materials that can act as anything from shells through to lenses. The early stages of calcium carbonate formation have recently attracted attention as stable prenucleation clusters have been observed, contrary to classical models. Here we show, using computer simulations combined with the analysis of experimental data, that these mineral clusters are made of an ionic polymer, composed of alternating calcium and carbonate ions, with a dynamic topology consisting of chains, branches and rings. The existence of a disordered, flexible and strongly hydrated precursor provides a basis for explaining the formation of other liquid-like amorphous states of calcium carbonate, in addition to the non-classical behaviour during growth of amorphous calcium carbonate.