Lead-free piezoelectric system (Na0.5Bi0.5)TiO3-BaTiO3: Equilibrium structures and irreversible structural transformations driven by electric field and mechanical impact

The structure-property correlation in the lead-free piezoelectric (1 - x)(Na0.5Bi0.5)TiO3-(x)BaTiO3 has been systematically investigated in detail as a function of composition (0 < x <= 0.11), temperature, electric field, and mechanical impact by Raman scattering, ferroelectric, piezoelectric measurement, x-ray, and neutron powder diffraction methods. Although x-ray diffraction study revealed three distinct composition ranges characterizing different structural features in the equilibrium state at room temperature: (i) monoclinic (Cc) + rhombohedral (R3c) for the precritical compositions, 0 <= x <= 0.05, (ii) cubiclike for 0.06 <= x <= 0.0675, and (iii) morphotropic phase boundary (MPB) like for 0.07 <= x < 0.10, Raman and neutron powder diffraction studies revealed identical symmetry for the cubiclike and the MPB compositions. The cubiclike structure undergoes irreversible phase separation by electric poling as well as by pure mechanical impact. This cubiclike phase exhibits relaxor ferroelectricity in its equilibrium state. The short coherence length (similar to 50A degrees) of the out-of-phase octahedral tilts does not allow the normal ferroelectric state to develop below the dipolar freezing temperature, forcing the system to remain in a dipolar glass state at room temperature. Electric poling helps the dipolar glass state to transform to a normal ferroelectric state with a concomitant enhancement in the correlation length of the out-of-phase octahedral tilt.

Comparison of the experimental, semi-experimental and ab initio equilibrium structures of acetylene: Influence of relativisitic effects and of the diagonal Born-Oppenheimer corrections

The equilibrium structure of acetylene (also named ethyne) has been reinvestigated to resolve the small discrepancies noted between different determinations. The size of the system as well as the large amount of available experimental data provides the quite unique opportunity to check the magnitude and relevance of various contributions to equilibrium structure as well as to verify the accuracy of experimental results. With respect to pure theoretical investigation, quantum-chemical calculations at the coupled-cluster level have been employed together with extrapolation to the basis set limit, consideration of higher excitations in the cluster operator, inclusion of core correlation effects as well as relativistic and diagonal Born-Oppenheimer corrections. In particular, it is found that the extrapolation to the complete basis set limit, the inclusion of higher excitations in the electronic-correlation treatment and the relativistic corrections are of the same order of magnitude. It also appears that a basis set as large as a core-valence quintuple-zeta set is required for accurately accounting for the inner-shell correlation contribution. From a pure experimental point of view, the equilibrium structure has been determined using very accurate rotational constants recently obtained by a global analysis (that is to say that all non-negligible interactions are explicitely included in the Hamiltonian matrix) of rovibrational spectra. Finally, a semi-experimental equilibrium structure (where the equilibrium rotational constants are obtained from the experimental ground state rotational constants and computed rovibrational corrections) has been obtained from the available experimental ground-state rotational constants for ten isotopic species corrected for computed vibrational corrections. Such a determination led to the revision of the ground-state rotational constants of two isotopologues, thus showing that structural determination is a good method to identify errors in experimental rotational constants. The three structures are found in a very good agreement, and our recommended values are rCC 120.2958(7) pm and rCH 106.164(1) pm. © 2011 American Institute of Physics.

Prediction of the equilibrium structures and photomagnetic properties of the Prussian blue analogue RbMn[Fe(CN)6] by density functional theory

A periodic density functional theory method using the B3LYP hybrid exchange-correlation potential is applied to the Prussian blue analogue RbMn[Fe(CN)6] to evaluate the suitability of the method for studying, and predicting, the photomagnetic behavior of Prussian blue analogues and related materials. The method allows correct description of the equilibrium structures of the different electronic configurations with regard to the cell parameters and bond distances. In agreement with the experimental data, the calculations have shown that the low-temperature phase (LT; Fe(2+)(t(6)2g, S = 0)-CN-Mn(3+)(t(3)2g e(1)g, S = 2)) is the stable phase at low temperature instead of the high-temperature phase (HT; Fe(3+)(t(5)2g, S = 1/2)-CN-Mn(2+)(t(3)2g e(2)g, S = 5/2)). Additionally, the method gives an estimation for the enthalpy difference (HT LT) with a value of 143 J mol(-1) K(-1). The comparison of our calculations with experimental data from the literature and from our calorimetric and X-ray photoelectron spectroscopy measurements on the Rb0.97Mn[Fe(CN)6]0.98 x 1.03 H2O compound is analyzed, and in general, a satisfactory agreement is obtained. The method also predicts the metastable nature of the electronic configuration of the high-temperature phase, a necessary condition to photoinduce that phase at low temperatures. It gives a photoactivation energy of 2.36 eV, which is in agreement with photoinduced demagnetization produced by a green laser.

Ultrafast diffraction and structural dynamics: The nature of complex molecules far from equilibrium

Studies of molecular structures at or near their equilibrium configurations have long provided information on their geometry in terms of bond distances and angles. Far-from-equilibrium structures are relatively unknown—especially for complex systems—and generally, neither their dynamics nor their average geometries can be extrapolated from equilibrium values. For such nonequilibrium structures, vibrational amplitudes and bond distances play a central role in phenomena such as energy redistribution and chemical reactivity. Ultrafast electron diffraction, which was developed to study transient molecular structures, provides a direct method for probing the nature of complex molecules far from equilibrium. Here we present our ultrafast electron diffraction observations of transient structures for two cyclic hydrocarbons. At high internal energies of ≈4 eV, these molecules display markedly different behavior. For 1,3,5-cycloheptatriene, excitation results in the formation of hot ground-state structures with bond distances similar to those of the initial structure, but with nearly three times the average vibrational amplitude. Energy is redistributed within 5 ps, but with a negative temperature characterizing the nonequilibrium population. In contrast, the ring-opening reaction of 1,3-cyclohexadiene is shown to result in hot structures with a C—C bond distance of over 1.7 Å, which is 0.2 Å away from any expected equilibrium value. Even up to 400 ps, energy remains trapped in large-amplitude motions comprised of torsion and asymmetric stretching. These studies promise a new direction for studying structural dynamics in nonequilibrium complex systems.

Effect of constraints by threonine on proline containing αa-helix—A molecular dynamics approach

Proline plays an important role in the secondary structure of proteins. In the pursuit of understanding its structural role, Proline containing helices with constraints have been studied by employing molecular dynamics (MD) technique. In the present study, the constraint introduced is a threonine residue, whose sidechain has intramolecular hydrogen bond interaction with the backbone oxygen atom. The three systems that have been chosen for characterization are: (1) Ace-(Ala)12−Thr-Pro-(Ala)10−NHMe, (2) Ace-(Ala)13-Pro-Ala-Thr- (Ala)8-NHMe and (3) Ace-(Ala)13-Pro-(Ala)3-Thr-(Ala)6-NHMe. The equilibrium structures and structural transitions have been identified by monitoring the backbone dihedral angles, bend related parameters and the hydrogen bond interactions. The MD averages and root mean square (r.m.s.) fluctuations are compared and discussed. Energy minimization has been carried out on selected MD simulated points in order to analyze the characteristics of different conformations.

Protein Structure and Function: Looking through the Network of Side-Chain Interactions

Network theory has become an excellent method of choice through which biological data are smoothly integrated to gain insights into complex biological problems. Understanding protein structure, folding, and function has been an important problem, which is being extensively investigated by the network approach. Since the sequence uniquely determines the structure, this review focuses on the networks of non-covalently connected amino acid side chains in proteins. Questions in structural biology are addressed within the framework of such a formalism. While general applications are mentioned in this review, challenging problems which have demanded the attention of scientific community for a long time, such as allostery and protein folding, are considered in greater detail. Our aim has been to explore these important problems through the eyes of networks. Various methods of constructing protein structure networks (PSN) are consolidated. They include the methods based on geometry, edges weighted by different schemes, and also bipartite network of protein-nucleic acid complexes. A number of network metrics that elegantly capture the general features as well as specific features related to phenomena, such as allostery and protein model validation, are described. Additionally, an integration of network theory with ensembles of equilibrium structures of a single protein or that of a large number of structures from the data bank has been presented to perceive complex phenomena from network perspective. Finally, we discuss briefly the capabilities, limitations, and the scope for further explorations of protein structure networks.

Density functional theory embedding for correlated wavefunctions

Methods that exploit the intrinsic locality of molecular interactions show significant promise in making tractable the electronic structure calculation of large-scale systems. In particular, embedded density functional theory (e-DFT) offers a formally exact approach to electronic structure calculations in which the interactions between subsystems are evaluated in terms of their electronic density. In the following dissertation, methodological advances of embedded density functional theory are described, numerically tested, and applied to real chemical systems.

First, we describe an e-DFT protocol in which the non-additive kinetic energy component of the embedding potential is treated exactly. Then, we present a general implementation of the exact calculation of the non-additive kinetic potential (NAKP) and apply it to molecular systems. We demonstrate that the implementation using the exact NAKP is in excellent agreement with reference Kohn-Sham calculations, whereas the approximate functionals lead to qualitative failures in the calculated energies and equilibrium structures.

Next, we introduce density-embedding techniques to enable the accurate and stable calculation of correlated wavefunction (CW) in complex environments. Embedding potentials calculated using e-DFT introduce the effect of the environment on a subsystem for CW calculations (WFT-in-DFT). We demonstrate that WFT-in-DFT calculations are in good agreement with CW calculations performed on the full complex.

We significantly improve the numerics of the algorithm by enforcing orthogonality between subsystems by introduction of a projection operator. Utilizing the projection-based embedding scheme, we rigorously analyze the sources of error in quantum embedding calculations in which an active subsystem is treated using CWs, and the remainder using density functional theory. We show that the embedding potential felt by the electrons in the active subsystem makes only a small contribution to the error of the method, whereas the error in the nonadditive exchange-correlation energy dominates. We develop an algorithm which corrects this term and demonstrate the accuracy of this corrected embedding scheme.

The limits of flexoelectricity in liquid crystals

The flexoelectric conversion of mechanical to electrical energy in nematic liquid crystals is investigated using continuum theory. Since the electrical energy produced cannot exceed the mechanical energy supplied, and vice-versa, upper bounds are imposed on the magnitudes of the flexoelectric coefficients in terms of the elastic and dielectric coefficients. For conventional values of the elastic and dielectric coefficients, it is shown that the flexoelectric coefficients may not be larger than a few tens of pC/m. This has important consequences for the future use of such flexoelectric materials in devices and the related energetics of distorted equilibrium structures. © 2011 Author(s).

Micromagnetic simulation of domain walls in iron films

The micromagnetic structure and energy of 180° domain walls spanning laminar crystals of iron having (100) or (110) surfaces and ranging in thickness from 145 to 580 nm have been investigated by numerical integration of the Landau-Lifshitz-Gilbert equation. Stable equilibrium structures with two flux symmetries were obtained for both crystal orientations at all thicknesses studied.

On Ti-18Si-6B and Ti-7.5Si-22.5B powder alloys prepared by high-energy ball milling and sintering

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

High-pressure assisted sintering of Ti-20Si-10B and Ti-10Si-5B powders

This work presents the structural characterization of Ti-10Si-5B and Ti-20Si-10B (at-%) alloys produced by high-pressure assisted sintering. Sintering was performed in air at 1100 and 1200°C for 60 s using pressure levels of 5 GPa. Structural evaluation of sintered samples was conducted by means of scanning electron microscopy and energy dispersive spectrometry. Samples were successfully consolidated after sintering, which presented theoretical density values higher than 99%. The microstructures of the sintered Ti-10Si-5B and Ti-20Si-10B alloys revealed the presence of the TiSS, TiB, TiB2, Ti5Si3, Ti5Si4, TiSi, and TiSi2.phases. A small amount of Ti6Si2B was formed after high-pressure assisted sintering of the Ti-20Si-10B alloy (5GPa, 1100°C for 60 s) indicating that equilibrium structures were not achieved during short sintering times. No oxygen and carbon contamination was detected in structures of Ti-Si-B alloys after high-pressure sintering at 1100 and 1200°C without controlled atmosphere. © (2012) Trans Tech Publications, Switzerland.

On the Mechanisms of Triplet Excited State Population in 8-Azaadenine

The photophysics of 8-azaadenine (8-AA) has been studied with the CASPT2//CASSCF protocol and ANO-L double-zeta basis sets. Stationary equilibrium structures, surface crossings, minimum energy paths, and linear interpolations have been used to study possible mechanisms to populate the lowest triplet state, T-1 (3)(pi pi*), capable of sensitizing molecular oxygen. Our results show that two main mechanisms can occur after photoexcitation to the S-2 (1)(pi pi*) state. The first one is through the S-2/S-1 conical intersection (((1)pi pi*/(1)n pi*)(Cl)), leading to the S-1 ((1)n pi*) state minimum, (S-1 ((1)n pi*))(min), where a singlet-triplet crossing, ((1)n pi*/(3)pi pi*)(STC), is accessible. The second one starts with the ((1)pi pi*/(3)n pi*)(STC) at the (S-2((1)pi pi*))(min), from which the system can evolve to the (T-2 ((3)n pi*))(min), with subsequent population of the T-1 excited electronic state, due to the ((3)n pi*/(3)pi pi*)(Cl) conical intersection.

On Ti-18Si-6B and Ti-7.5Si-22.5B powder alloys prepared by high-energy ball milling and sintering

Recently, a new ternary phase was discovered in the Ti-Si-B system, located near the Ti6Si2B composition. The present study concerns the preparation of titanium alloys that contain such phase mixed with α-titanium and other intermetallic phases. High-purity powders were initially processed in a planetary ball-mill under argon atmosphere with Ti-18Si-6B and Ti-7.5Si-22.5B at. (%) initial compositions. Variation of parameters such as rotary speed, time, and ball diameters were adopted. The as-milled powders were pressureless sintered and hot pressed. Both the as-milled and sintered materials were characterized by X-ray diffraction, scanning electron microscopy and energy-dispersive spectrometry. Sintered samples have presented equilibrium structures formed mainly by the α-Ti+Ti6Si2B+Ti5Si3+TiB phases. Silicon and boron peaks disappear throughout the milling processes, as observed in the powder diffraction data. Furthermore, an iron contamination of up to 10 at. (%) is measured by X-ray spectroscopy analysis on some regions of the sintered samples. Density, hardness and tribological results for these two compositions are also presented here.

Interpolyelektrolytkomplexe zylindrischer Bürstenpolymere

Die elektrostatische Wechselwirkung zwischen entgegengesetzt geladenen Polyelektrolyten führt zur spontanen Bildung von Interpolyelektrolytkomplexen. Besonders im Fokus des akademischen und biotechnologischen Interesses stehen derzeit Komplexe aus DNA und synthetischen Polykationen, da eine Anwendung dieser speziellen Interpolyelektrolytkomplexe in der nicht-viralen Gentherapie vielfältig diskutiert wird. Ziel der vorliegenden Arbeit war es, den Einfluss von Kettensteifheit auf die Bildung von Interpolyelektrolytkomplexen zu untersuchen und dabei Wege aufzuzeigen, die eine Kontrolle der Topologie von Interpolyelektrolytkomplexen ermöglichen. Neben dem topologischen Einfluss wurde untersucht, wie durch Komplexierung hochmolekularer Polyelektrolyte equilibrierbare Strukturen erhalten werden können. Als Modellsystem für diese Untersuchungen wurden zylindrische Bürstenpolymere verwendet, denen als topologischer „Kontrast“ das Komplexierungsverhalten kommerzieller PAMAM-G5-Dendrimere mit kugelförmiger Topologie gegenüber gestellt wurde. Um den Ladungsgrad des Bürstenpolymers beliebig variieren zu können, wurden Bürstenpolymere mit Poly(ethylenimin)-Seitenketten synthetisiert, deren Ladungsdichte über den Protonierungsgrad einstellbar ist. Die vorliegende Arbeit zeigt, wie diese mit Hilfe der Makromonomermethode hergestellt werden konnten. Die Komplexbildung von DNA mit semiflexiblen zylindrischen Bürstenpolymeren mit unterschiedlichen Seitenketten und Ladungsdichten in wässriger Lösung hat gezeigt, dass diese in allen untersuchten Fällen unter kinetischer Kontrolle verläuft und Nicht-Gleichgewichtsstrukturen gebildet werden. Sehr überraschend wurde festgestellt, dass die Größen der mit vorgelegter DNA gebildeten Komplexe ungeachtet des verwendeten Polykations identisch sind und DNA-Komplexe mit einem Radius von 30 bis 50 nm und einer kugelförmigen Topologie resultieren. Diese kinetisch kontrollierte Komplexbildung konnte in nicht-wässriger Lösung durch starke Reduktion der Anzahl wechselwirkender Ladungen verhindert werden, so dass eine thermodynamische Kontrolle möglich wurde. Unter diesen Bedingungen ist es gelungen, aus hochgeladenen Poly(styrolsulfonat)-Bürsten mit modifizierten Poly(ethylenimin)-Bürsten oder PAMAM-Dendrimeren Komplexe zylindrischer Topologie herzu-stellen. Für letztere konnte darüber hinaus postuliert werden, dass diese Komplexe eine dichteste Packung der PAMAM-Dendrimere darstellen, für deren Bildung das Polyanion mit seiner größeren Konturlänge und seiner zylindrischen Topologie als Templat dient.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.