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.

Liquid crystal necklaces: cholesteric drops threaded by thin cellulose fibres

Liquid crystals in confined geometries exhibit numerous complex structures often including topological defects that are controlled by the nematic elasticity, chirality and surface anchoring. In this work, we study the structures of cholesteric droplets pierced by cellulose fibres with planar anchoring at droplet and fibre surfaces. By varying the temperature we demonstrate the role of twisting power and droplet diameter on the equilibrium structures. The observed structures are complemented by detailed numerical simulations of possible director fields decorated by defects. Three distinct structures, a bipolar and two ring configurations, are identified experimentally and numerically. Designing cholesteric liquid crystal microdroplets on thin long threads opens new routes to produce fibre waveguides decorated with complex microresonators.

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.

Störungstheoretische Behandlung relativistischer Korrekturen zu Gleichgewichtsstrukturen im Rahmen der Coupled-Cluster-Theorie

Quantenchemische Untersuchungen von Atomen und Molekülen haben in den letzten Jahren durch die systematische Erweiterung der Methoden und Computerresourcen zunehmend für die Interpretation und Vorhersage experimenteller Ergebnisse an Bedeutung gewonnen. Relativistische Effekte in der Chemie werden zum Beispiel für die gelbe Farbe von Gold und den flüssigen Aggregatzustand von Quecksilber verantwortlich gemacht und müssen daher in quantenchemischen Rechnungen berücksichtigt werden. Relativistische Effekte sind bei leichten Elementen oft so klein, daß sie in vielen quantenchemischen Betrachtungen vernachlässigt werden. Dennoch sind es gerade diese Beiträge, die verbleibende Abweichungen von noch so genauen nichtrelativistischen Rechnungen von ebenso genauen experimentellen Ergebnissen ausmachen können. Relativistische Effekte können auf viele Arten in quantenchemischen Rechnungen berücksichtigt werden. Eine Möglichkeit ist die Störungstheorie. Ein derartiger Ansatz ist die Mass-velocity-Darwin-Näherung, ein anderer die Direkte Störungstheorie. Hier entspricht die relativistische Energiekorrektur erster Ordnung der ersten Ableitung der Energie nach einem relativistischen Störparameter. Für eine Bestimmung der Gleichgewichtsstruktur eines Moleküls müssen die Kräfte auf die Atomkerne bestimmt werden. Diese entsprechen einer ersten Ableitung der Gesamtenergie nach den Kernkoordinaten. Eine Einbeziehung der relativistischen Effekte auf diese Kräfte erfordert daher die gemischte zweite Ableitung der Energie nach dem relativistischen Störparameter und den Kernkoordinaten. Diese relativistischen Korrekturen wurden in dem quantenchemischen Programmpaket ACES2 implementiert. Ein Resultat dieser Arbeit ist, daß nun erstmalig eine Implementierung analytischer Gradienten für die Berechnung relativistischer Korrekturen zu Strukturparametern mit Hilfe der relativistischen Störungstheorie für den Coupled-Cluster-Ansatz bereit steht. Die Coupled-Cluster-Theorie eignet sich besonders gut für die hochgenaue Vorhersage von molekularen Eigenschaften, wie der Gleichgewichtsstruktur. Im Rahmen dieser Arbeit wurde die Basissatzabhängigkeit der relativistischen Beiträge zu Energien, Strukturparametern und harmonischen Schwingungsfrequenzen im Detail untersucht. Für die hier untersuchten Moleküle sind die relativistischen Effekte und Effekte aufgrund der Elektronenkorrelation nicht additiv, so verkürzt die Berücksichtigung relativistischer Effekte bei Hartree-Fock-Rechnungen die Bindung in den Hydrogenhalogeniden, während die Einbeziehung der Elektronenkorrelation durch CCSD(T)-Rechnungen zu einer verlängerten Bindung im Fluorwasserstoff und weniger stark ausgeprägten Korrekturen im Chlor- und Bromwasserstoff führt. Für die anderen hier untersuchten mehratomigen Moleküle findet sich kein einheitlicher Trend; dies unterstreicht die Notwendigkeit expliziter Rechnungen. Damit steht ein leistungsfähiges und vielseitiges Werkzeug für die Berechnung relativistischer Korrekturen auf verschiedenste molekulare Eigenschaften zur Verfügung, das mit modernen, systematisch verbesserbaren quantenchemischen Methoden verknüpft ist. Hiermit ist es möglich, hochgenaue Rechnungen zur Vorhersage und Interpretation von Experimenten durchzuführen.

A soft-tetramer model for diblock copolymer melts: structure formation and geometric characterization

The ability of block copolymers to spontaneously self-assemble into a variety of ordered nano-structures not only makes them a scientifically interesting system for the investigation of order-disorder phase transitions, but also offers a wide range of nano-technological applications. The architecture of a diblock is the most simple among the block copolymer systems, hence it is often used as a model system in both experiment and theory. We introduce a new soft-tetramer model for efficient computer simulations of diblock copolymer melts. The instantaneous non-spherical shape of polymer chains in molten state is incorporated by modeling each of the two blocks as two soft spheres. The interactions between the spheres are modeled in a way that the diblock melt tends to microphase separate with decreasing temperature. Using Monte Carlo simulations, we determine the equilibrium structures at variable values of the two relevant control parameters, the diblock composition and the incompatibility of unlike components. The simplicity of the model allows us to scan the control parameter space in a completeness that has not been reached in previous molecular simulations.The resulting phase diagram shows clear similarities with the phase diagram found in experiments. Moreover, we show that structural details of block copolymer chains can be reproduced by our simple model.We develop a novel method for the identification of the observed diblock copolymer mesophases that formalizes the usual approach of direct visual observation,using the characteristic geometry of the structures. A cluster analysis algorithm is used to determine clusters of each component of the diblock, and the number and shape of the clusters can be used to determine the mesophase.We also employ methods from integral geometry for the identification of mesophases and compare their usefulness to the cluster analysis approach.To probe the properties of our model in confinement, we perform molecular dynamics simulations of atomistic polyethylene melts confined between graphite surfaces. The results from these simulations are used as an input for an iterative coarse-graining procedure that yields a surface interaction potential for the soft-tetramer model. Using the interaction potential derived in that way, we perform an initial study on the behavior of the soft-tetramer model in confinement. Comparing with experimental studies, we find that our model can reflect basic features of confined diblock copolymer melts.

Higher-order molecular properties and excitation energies in single-reference and multireference coupled-cluster theory

Coupled-cluster (CC) theory is one of the most successful approaches in high-accuracy quantum chemistry. The present thesis makes a number of contributions to the determination of molecular properties and excitation energies within the CC framework. The multireference CC (MRCC) method proposed by Mukherjee and coworkers (Mk-MRCC) has been benchmarked within the singles and doubles approximation (Mk-MRCCSD) for molecular equilibrium structures. It is demonstrated that Mk-MRCCSD yields reliable results for multireference cases where single-reference CC methods fail. At the same time, the present work also illustrates that Mk-MRCC still suffers from a number of theoretical problems and sometimes gives rise to results of unsatisfactory accuracy. To determine polarizability tensors and excitation spectra in the MRCC framework, the Mk-MRCC linear-response function has been derived together with the corresponding linear-response equations. Pilot applications show that Mk-MRCC linear-response theory suffers from a severe problem when applied to the calculation of dynamic properties and excitation energies: The Mk-MRCC sufficiency conditions give rise to a redundancy in the Mk-MRCC Jacobian matrix, which entails an artificial splitting of certain excited states. This finding has established a new paradigm in MRCC theory, namely that a convincing method should not only yield accurate energies, but ought to allow for the reliable calculation of dynamic properties as well. In the context of single-reference CC theory, an analytic expression for the dipole Hessian matrix, a third-order quantity relevant to infrared spectroscopy, has been derived and implemented within the CC singles and doubles approximation. The advantages of analytic derivatives over numerical differentiation schemes are demonstrated in some pilot applications.

Conformational transitions of glycine induced by vibrational excitation of the O-H stretch

Vibrational energy flow and conformational transitions following excitation of the OH stretching mode of the most stable conformer of glycine are studied by classical trajectories. "On the fly" simulations with the PM3 semiempirical electronic structure method for the potential surface are used. Initial conditions are selected to correspond to the v = 1 excitation of the OH stretch. The main findings are: (1) An an equilibrium-like ratio is established between the populations of the 3 lowest-lying conformers after about 10 picoseconds. (2) There is a high probability throughout the 150 ps of the simulations for finding the molecule in geometries far from the equilibrium structures of the lowest-energy conformers. (3) Energy from the initial excited OH (v = 1) stretch flows preferentially to 5 other vibrational modes, including the bending motion of the H atom. (4) RRK theory yields conformational transition rates that deviate substantially from the classical trajectory results. Possible implication of these results for vibrational energy flow and conformational transitions in small biological molecules are discussed.

Intermolecular Clamping by Hydrogen Bonds: 2-Pyridone center dot NH3

A combined spectroscopic and ab initio theoretical study of the doubly hydrogen-bonded complex of 2-pyridone (2PY) with NH3 has been performed. The S-1 <- S-0 spectrum extends up to approximate to 1200 cm(-1) above the 0(0)(0) band, close to twice the range observed for 2PY. The S-1 state nonradiative decay for vibrations above approximate to 300 cm(-1) in the NH3 complex is dramatically slowed down relative to bare 2PY. Also, the Delta v=2,4,... overtone bands of the v(1)' and v(2)' out-of-plane vibrations that dominate the low-energy spectral region of 2PY are much weaker or missing for 2PY center dot NH3, which implies that the bridging (2PY)NH center dot center dot center dot NH3 and H2NH center dot center dot center dot O=C H-bonds clamp the 2PY at a planar geometry in the S-1 state. The mass-resolved UV vibronic spectra of jet-cooled 2PY center dot NH3 and its H/D mixed isotopomers are measured using two-color resonant two-photon ionization spectroscopy. The S-0 and S-1 equilibrium structures and normal-mode frequencies are calculated by density functional (B3LYP) and correlated ab initio methods (MP2 and approximate second-order coupled-cluster, CC2). The S-1 <- S-0 vibronic assignments are based on configuration interaction singles (CIS) and CC2 calculations. A doubly H-bonded bridged structure of C-S symmetry is predicted, in agreement with that of Held and Pratt [J. Am. Chem. Soc. 1993, 115, 9718]. While the B3LYP and MP2 calculated rotational constants are in very good agreement with experiment, the calculated H2NH center dot center dot center dot O=C H-bond distance is approximate to 0.7 angstrom shorter than that derived by Held and Pratt. On the other hand, this underlines their observation that ammonia can act as a strong H-bond donor when built into an H-bonded bridge. The CC2 calculations predict the H2NH center dot center dot center dot O distance to increase by 0.2 angstrom upon S-1 <- S-0 electronic excitation, while the (2PY)NH center dot center dot center dot NH3 H-bond remains nearly unchanged. Thus, the expansion of the doubly H-bonded bridge in the excited state is asymmetric and almost wholly due to the weakening of the interaction of ammonia with the keto acceptor group.

The antimicrobial action of pseudin-2

The amphibian antimicrobial peptide pseudin-2 is a peptide derived from the skin of the South-American frog Pseudis paradoxa (Olson et al., 2001). This peptide possesses tremendous potential as a therapeutic lead since it has been shown to possess both antimicrobial as well insulin-releasing properties (Olson et al., 2001; Abdel-Wahab et al., 2008). This study aimed to develop pseudin-2’s potential by understanding and improving its properties as an antimicrobial agent. The structure-function relationships of pseudin-2 were explored using a combination of in-vitro and in-silico techniques, with an aim to predict how the structure of the peptide may be altered in order to improve its efficacy. A library of pseudin-2 mutants was generated by randomizing codons at positions 10, 14 and 18 of a synthetic gene, using NNK saturation mutagenesis. Analysis of these novel peptides broadly confirmed, in line with literature precedent, that anti-microbial activity increases with increased positive charge. Specifically, 2 positively-charged residues at positions 10 and 14 and a hydrophobic at position 18 are preferred. However, substitution at position 14 with some polar, non-charged residues also created peptides with antimicrobial activity. Interestingly, the pseudin-2 analogue [10-E, 14-Q, 18-L] which is identical to pseudin-2, except that the residues at positions 10 and 14 are switched, showed no anti-microbial activity at all. Molecular dynamics simulations of pseudin-2 showed that the peptide possesses two equilibrium structures in a membrane environment: a linear and a kinked a-helix which both embed into the membrane at an angle. Biophysical characterization using circular dichroism spectroscopy confirmed that the peptide is helical within the membrane environment whilst linear dichroism established that the peptide has no defined orientation within the membrane. Collectively, these data indicate that Pseudin-2 exerts its antimicrobial activity via the carpet model.