A hybrid approach to Fermi operator expansion

In a recent paper we have suggested that the finite temperature density matrix can be computed efficiently by a combination of polynomial expansion and iterative inversion techniques. We present here significant improvements over this scheme. The original complex-valued formalism is turned into a purely real one. In addition, we use Chebyshev polynomials expansion and fast summation techniques. This drastically reduces the scaling of the algorithm with the width of the Hamiltonian spectrum, which is now of the order of the cubic root of such parameter. This makes our method very competitive for applications to ab-initio simulations, when high energy resolution is required.Comment: preprint of ICCMSE08 proceeding

Evolution of the structure of amorphous ice - from low-density amorphous (LDA) through high-density amorphous (HDA) to very high-density amorphous (VHDA) ice

We report results of molecular dynamics simulations of amorphous ice for pressures up to 22.5 kbar. The high-density amorphous ice (HDA) as prepared by pressure-induced amorphization of Ih ice at T=80 K is annealed to T=170 K at various pressures to allow for relaxation. Upon increase of pressure, relaxed amorphous ice undergoes a pronounced change of structure, ranging from the low-density amorphous ice (LDA) at p=0, through a continuum of HDA states to the limiting very high-density amorphous ice (VHDA) regime above 10 kbar. The main part of the overall structural change takes place within the HDA megabasin, which includes a variety of structures with quite different local and medium-range order as well as network topology and spans a broad range of densities. The VHDA represents the limit to densification by adapting the hydrogen-bonded network topology, without creating interpenetrating networks. The connection between structure and metastability of various forms upon decompression and heating is studied and discussed. We also discuss the analogy with amorphous and crystalline silica. Finally, some conclusions concerning the relation between amorphous ice and supercooled water are drawn.Comment: 11 pages, 12 postscript figures. To be published in The Journal of Chemical Physic

An Efficient and Accurate Car-Parrinello-like Approach to Born-Oppenheimer Molecular Dynamics

We present a new method which combines Car-Parrinello and Born-Oppenheimer molecular dynamics in order to accelerate density functional theory based ab-initio simulations. Depending on the system a gain in efficiency of one to two orders of magnitude has been observed, which allows ab-initio molecular dynamics of much larger time and length scales than previously thought feasible. It will be demonstrated that the dynamics is correctly reproduced and that high accuracy can be maintained throughout for systems ranging from insulators to semiconductors and even to metals in condensed phases. This development considerably extends the scope of ab-initio simulations.Comment: 4 pages, 3 figures; Accepted by Phys. Rev. Lett. for publicatio

Coexistence of tetrahedral and octahedral-like sites in amorphous phase change materials

Chalcogenide alloys are materials of interest for optical recording and non-volatile memories. We perform ab-initio molecular dynamics simulations aiming at shading light onto the structure of amorphous Ge2Sb2Te5 (GST), the prototypical material in this class. First principles simulations show that amorphous GST obtained by quenching from the liquid phase displays two types of short range order. One third of Ge atoms are in a tetrahedral environment while the remaining Ge, Sb and Te atoms display a defective octahedral environment, reminiscent of cubic crystalline GST.Comment: 6 pages, 10 figures including supplementary materia

Freezing of a Lennard-Jones fluid: from nucleation to spinodal regine

Using molecular dynamics, we investigate the crystal nucleation in a Lennard-Jones fluid as a function of the degree of supercooling. At moderate supercooling, a nucleation picture applies, while for deeper quenches, the phenomenon progressively acquires a spinodal character. We show that in the nucleation regime, the freezing is a two-step process. The formation of the critical nucleus is indeed preceded by the abrupt formation of a precritical crystallite from a density fluctuation in the fluid. In contrast, as the degree of supercooling is increased, crystallization proceeds in a more continuous and collective fashion and becomes more spatially diffuse, indicating that the liquid is unstable and crystallizes by a spinodal mechanism

The Ciona intestinalis immune-related galectin genes (CiLgals-a and CiLgals-b) are expressed by the gastric epithelium

The transcription of two Ciona intestinalis galectin genes (CiLgals-a and CiLgalseb) is uparegulated by LPS in the pharynxis (hemocytes, vessel epithelium, endostilar zones) which is retained the main organ of the immunity. In this ascidian, for the first time we show, by immunohistochemistry and in situ hybridization methods, that these two immune-related genes are expressed in the gastric epithelium of na\uefve ascidians, whereas the galectins appear to be only contained in the intestine columnar epithelium. In addition, according to previous results on the pharynx, the genes are also expressed and galectins produced by hemocytes scattered in the connective tissue surrounding the gut. The genes expression and galectin localization in several tissues, including the previous findings on the transcription upregulation, the constitutive expression of these genes by endostylar zones and by the gastric epithelium suggest a potential multifunctional role of these galectins. In this respect, it is of interest to define where the CiLgals are normally found as related to the tissue functions. Such an approach should be a starting point for further investigations

Upregulated transcription of phenoloxidase genes in the pharynx and endostyle of Ciona intestinalis in response to LPS

We investigated the role of phenoloxidases (POs) in ascidians inflammatory reaction, a components of a copper-containing protein family involved in invertebrate immune system. In Ciona intestinalis two phenoloxidases (CinPO-1, CinPO-2) have been sequenced. In the present study, real time PCR analysis showed that both CinPO-1 and CinPO-2 genes were modulated by LPS inoculation suggesting that they are inducible and highly expressed in the inflamed pharynx. In situ hybridization disclosed CinPO-1 and CinPO-2 transcripts in pharynx hemocytes (granulocytes) and, mainly, in unilocular refractile granulocytes (URG) which mainly populated the inflamed tunic matrix. Interestingly, the genes are also upregulated by LPS in the endostyle (zones 7, 8 and 9) that is considered homolog to the vertebrate thyroid