Stabilization of Ab Initio Molecular Dynamics Simulations at Large Time Steps

The Verlet method is still widely used to integrate the equations of motion in ab initio molecular dynamics simulations. We show that the stability limit of the Verlet method may be significantly increased by setting an upper limit on the kinetic energy of each atom with only a small loss in accuracy. The validity of this approach is demonstrated for molten lithium fluoride.Comment: 9 pages, 3 figure

Ideal, Defective, and Gold--Promoted Rutile TiO2(110) Surfaces: Structures, Energies, Dynamics, and Thermodynamics from PBE+U

Extensive first principles calculations are carried out to investigate gold-promoted TiO2(110) surfaces in terms of structure optimizations, electronic structure analyses, ab initio thermodynamics calculations of surface phase diagrams, and ab initio molecular dynamics simulations. All computations rely on density functional theory in the generalized gradient approximation (PBE) and account for on-site Coulomb interactions via inclusion of a Hubbard correction, PBE+U, where U is computed from linear response theory. This approach is validated by investigating the interaction between TiO2(110) surfaces and typical probe species (H, H2O, CO). Relaxed structures and binding energies are compared to both data from the literature and plain PBE results. The main focus of the study is on the properties of gold-promoted titania surfaces and their interactions with CO. Both PBE+U and PBE optimized structures of Au adatoms adsorbed on stoichiometric and reduced TiO2 surfaces are computed, along with their electronic structure. The charge rearrangement induced by the adsorbates at the metal/oxide contact are also analyzed and discussed. By performing PBE+U ab initio molecular dynamics simulations, it is demonstrated that the diffusion of Au adatoms on the stoichiometric surface is highly anisotropic. The metal atoms migrate either along the top of the bridging oxygen rows, or around the area between these rows, from one bridging position to the next along the [001] direction. Approximate ab initio thermodynamics predicts that under O-rich conditions, structures obtained by substituting a Ti5c atom with an Au atom are thermodynamically stable over a wide range of temperatures and pressures.Comment: 20 pages, 12 figures, accepted for publication in Phys. Rev.

Quantum Fluctuations Driven Orientational Disordering: A Finite-Size Scaling Study

The orientational ordering transition is investigated in the quantum generalization of the anisotropic-planar-rotor model in the low temperature regime. The phase diagram of the model is first analyzed within the mean-field approximation. This predicts at $T=0$ a phase transition from the ordered to the disordered state when the strength of quantum fluctuations, characterized by the rotational constant $\Theta$, exceeds a critical value $\Theta_{\rm c}^{MF}$. As a function of temperature, mean-field theory predicts a range of values of $\Theta$ where the system develops long-range order upon cooling, but enters again into a disordered state at sufficiently low temperatures (reentrance). The model is further studied by means of path integral Monte Carlo simulations in combination with finite-size scaling techniques, concentrating on the region of parameter space where reentrance is predicted to occur. The phase diagram determined from the simulations does not seem to exhibit reentrant behavior; at intermediate temperatures a pronounced increase of short-range order is observed rather than a genuine long-range order.Comment: 27 pages, 8 figures, RevTe

Fixed-parameter tractability of multicut parameterized by the size of the cutset

Given an undirected graph $G$, a collection $\{(s_1,t_1),..., (s_k,t_k)\}$ of pairs of vertices, and an integer $p$, the Edge Multicut problem ask if there is a set $S$ of at most $p$ edges such that the removal of $S$ disconnects every $s_i$ from the corresponding $t_i$. Vertex Multicut is the analogous problem where $S$ is a set of at most $p$ vertices. Our main result is that both problems can be solved in time $2^{O(p^3)}... n^{O(1)}$, i.e., fixed-parameter tractable parameterized by the size $p$ of the cutset in the solution. By contrast, it is unlikely that an algorithm with running time of the form $f(p)... n^{O(1)}$ exists for the directed version of the problem, as we show it to be W[1]-hard parameterized by the size of the cutset

Time-reversible Born-Oppenheimer molecular dynamics

We present a time-reversible Born-Oppenheimer molecular dynamics scheme, based on self-consistent Hartree-Fock or density functional theory, where both the nuclear and the electronic degrees of freedom are propagated in time. We show how a time-reversible adiabatic propagation of the electronic degrees of freedom is possible despite the non-linearity and incompleteness of the self-consistent field procedure. Time-reversal symmetry excludes a systematic long-term energy drift for a microcanonical ensemble and the number of self-consistency cycles can be kept low (often only 2-4 cycles per nuclear time step) thanks to a good initial guess given by the adiabatic propagation of the electronic degrees of freedom. The time-reversible Born-Oppenheimer molecular dynamics scheme therefore combines a low computational cost with a physically correct time-reversible representation of the dynamics, which preserves a detailed balance between propagation forwards and backwards in time.Comment: 4 pages, 4 figure

Epitaxial growth and transport properties of Sr2_2CrWO6_6 thin films

We report on the preparation and characterization of epitaxial thin films of the double-perovskite Sr$_2$CrWO$_6$ by Pulsed Laser Deposition (PLD). On substrates with low lattice mismatch like SrTiO$_3$, epitaxial Sr$_2$CrWO$_6$ films with high crystalline quality can be grown in a molecular layer-by-layer growth mode. Due to the similar ionic radii of Cr and W, these elements show no sublattice order. Nevertheless, the measured Curie temperature is well above 400 K. Due to the reducing growth atmosphere required for double perovskites, the SrTiO$_3$ substrate surface undergoes an insulator-metal transition impeding the separation of thin film and substrate electric transport properties.Comment: 3 pages, 5 figure

Parameterized Inapproximability of Target Set Selection and Generalizations

In this paper, we consider the Target Set Selection problem: given a graph and a threshold value $thr(v)$ for any vertex $v$ of the graph, find a minimum size vertex-subset to "activate" s.t. all the vertices of the graph are activated at the end of the propagation process. A vertex $v$ is activated during the propagation process if at least $thr(v)$ of its neighbors are activated. This problem models several practical issues like faults in distributed networks or word-to-mouth recommendations in social networks. We show that for any functions $f$ and $\rho$ this problem cannot be approximated within a factor of $\rho(k)$ in $f(k) \cdot n^{O(1)}$ time, unless FPT = W[P], even for restricted thresholds (namely constant and majority thresholds). We also study the cardinality constraint maximization and minimization versions of the problem for which we prove similar hardness results

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.Comment: 6 pages, 3 figure

Parameterized Approximation Schemes using Graph Widths

Combining the techniques of approximation algorithms and parameterized complexity has long been considered a promising research area, but relatively few results are currently known. In this paper we study the parameterized approximability of a number of problems which are known to be hard to solve exactly when parameterized by treewidth or clique-width. Our main contribution is to present a natural randomized rounding technique that extends well-known ideas and can be used for both of these widths. Applying this very generic technique we obtain approximation schemes for a number of problems, evading both polynomial-time inapproximability and parameterized intractability bounds