MA160exam00pdf

Exam questions and solutions in PDF

Exercise sheet 3 pdf

Exercises and solutions in PDF

Exercise sheet 2 pdf

Exercises and solutions in PDF

Anharmonicity contributions to the vibrational second hyperpolarizability of conjugated oligomers

Restricted Hartree-Fock 6-31G calculations of electrical and mechanical anharmonicity contributions to the longitudinal vibrational second hyperpolarizability have been carried out for eight homologous series of conjugated oligomers - polyacetylene, polyyne, polydiacetylene, polybutatriene, polycumulene, polysilane, polymethineimine, and polypyrrole. To draw conclusions about the limiting infinite polymer behavior, chains containing up to 12 heavy atoms along the conjugated backbone were considered. In general, the vibrational hyperpolarizabilities are substantial in comparison with their static electronic counterparts for the dc-Kerr and degenerate four-wave mixing processes (as well as for static fields) but not for electric field-induced second harmonic generation or third harmonic generation. Anharmonicity terms due to nuclear relaxation are important for the dc-Kerr effect (and for the static hyperpolarizability) in the σ-conjugated polymer, polysilane, as well as the nonplanar π systems polymethineimine and polypyrrole. Restricting polypyrrole to be planar, as it is in the crystal phase, causes these anharmonic terms to become negligible. When the same restriction is applied to polymethineimine the effect is reduced but remains quantitatively significant due to the first-order contribution. We conclude that anharmonicity associated with nuclear relaxation can be ignored, for semiquantitative purposes, in planar π-conjugated polymers. The role of zero-point vibrational averaging remains to be evaluated

Effect of basis set superposition error on the electron density of molecular complexes

The effect of basis set superposition error (BSSE) on molecular complexes is analyzed. The BSSE causes artificial delocalizations which modify the first order electron density. The mechanism of this effect is assessed for the hydrogen fluoride dimer with several basis sets. The BSSE-corrected first-order electron density is obtained using the chemical Hamiltonian approach versions of the Roothaan and Kohn-Sham equations. The corrected densities are compared to uncorrected densities based on the charge density critical points. Contour difference maps between BSSE-corrected and uncorrected densities on the molecular plane are also plotted to gain insight into the effects of BSSE correction on the electron density

Initial convergence of the perturbation series expansion for vibrational nonlinear optical properties

Initial convergence of the perturbation series expansion for vibrational nonlinear optical (NLO) properties was analyzed. The zero-point vibrational average (ZPVA) was obtained through first-order in mechanical plus electrical anharmonicity. Results indicated that higher-order terms in electrical and mechanical anharmonicity can make substantial contributions to the pure vibrational polarizibility of typical NLO molecules

Estructura computacional i aplicacions de la semblança molecular quàntica

La tesis tracta diferents aspectes relacionats amb el càlcul de la semblança quàntica, així com la seva aplicació en la racionalització i predicció de l'activitat de fàrmacs. Es poden destacar dos progressos importants en el desenvolupament de noves metodologies que faciliten el càlcul de les mesures de semblança quàntica. En primer lloc, la descripció de les molècules mitjançant les funciones densitat aproximades PASA (Promolecular Atomic Shell Approximation) ha permès descriure amb suficient precisió la densitat electrònica dels sistemes moleculars analitzats, reduint substancialment el temps de càlcul de les mesures de semblança. En segon lloc, el desenvolupament de tècniques de superposició molecular específiques de les mesures de semblança quàntica ha permès resoldre el problema de l'alineament en l'espai dels compostos comparats. El perfeccionament d'aquests nous procediments i algoritmes matemàtics associats a les mesures de semblança molecular quàntica, ha estat essencial per poder progressar en diferents disciplines de la química computacional, sobretot les relacionades amb les anàlisis quantitatives entre les estructures moleculars i les seves activitats biològiques, conegudes amb les sigles angleses QSAR (Quantitative Structure-Activity Relationships). Precisament en l'àrea de les relacions estructura-activitat s'han presentat dues aproximacions fonamentades en la semblança molecular quàntica que s'originen a partir de dues representacions diferents de les molècules. La primera descripció considera la densitat electrònica global de les molècules i és important, entre altres, la disposició dels objectes comparats en l'espai i la seva conformació tridimensional. El resultat és una matriu de semblança amb les mesures de semblança de tots els parells de compostos que formen el conjunt estudiat. La segona descripció es fonamenta en la partició de la densitat global de les molècules en fragments. S'utilitzen mesures d'autosemblança per analitzar els requeriments bàsics d'una determinada activitat des del punt de vista de la semblança quàntica. El procés permet la detecció de les regions moleculars que són responsables d'una alta resposta biològica. Això permet obtenir un patró amb les regions actives que és d'evident interès per als propòsits del disseny de fàrmacs. En definitiva, s'ha comprovat que mitjançant la simulació i manipulació informàtica de les molècules en tres dimensions es pot obtenir una informació essencial en l'estudi de la interacció entre els fàrmacs i els seus receptors macromoleculars.

New mono- and dinuclear ruthenium complexes containing the 3,5-bis(2-pyridyl)pyrazole ligand. Synthesis, characterization and applications

Se han sintetizado dos nuevos complejos mononucleares de Ru, con formula [RuCl2(Hbpp)(dmso)2], a partir de la reacción entre [RuCl2(dmso)4] y Hbpp (3,5-bis(2-piridil)pirazola). El hecho que sólo tres de los seis posibles estereoisómeros se obtengan a partir de esta reacción, se ha racionalizado en base a factores estructurales y electrónicos. Estos complejos se han caracterizado de forma estructural, espectroscópica y electroquímica. En acetonitrilo en medio básico, el isómero trans,cis-[RuCl2(Hbpp)(dmso)2] da lugar a procesos de isomerización de enlace de un ligando dmso cuando el Ru(II) se oxida a Ru(III). Las constantes termodinámicas y cinéticas para el proceso se han determinado por voltametria cíclica. La irradiación de trans,cis-[RuCl2(Hbpp)(dmso)2] y cis(out),cis-[RuCl2(Hbpp)(dmso)2] con luz UV o solar da lugar a reacciones de fotosustitución de un ligando dmso por una molécula de acetonitrilo para dar un nuevo compuesto el cual ha sido caracterizado en solución por técnicas espectroscópicas y electroquímicas. Ambos complejos resultan catalizadores útiles en la transferencia de hidrógeno de isopropanol a acetofenona, obteniéndose 1-feniletanol como único producto y un 42.1% de conversión (36.1 ciclos metálicos) a 80ºC con el isómero trans,cis-[RuCl2(Hbpp)(dmso)2], que resulta significativamente más eficaz que el complejo cis(out),cis-[RuCl2(Hbpp)(dmso)2]. La reacción de cis(out),cis-[RuCl2(Hbpp)(dmso)2] con trpy (2,2':6',2"-terpiridina) da lugar a los dos isómeros geométricos del complejo [Ru(Hbpp)(trpy)(Cl)]+, el in y el out. Estos complejos se han aislado y caracterizado por técnicas estructurales, espectroscópicas y electroquímicas. Estos cloro complejos han sido utilizados como precursores para la síntesis de los complejos análogos con ligandos aqua (in,out-[Ru(Hbpp)(trpy)(H2O)]2+) y piridina (in,out-[Ru(Hbpp)(trpy)(py)]2+), los cuales también han sido aislados y caracterizados. Las propiedades ácido-base de los aqua complejos, y del complejo out-py se han estudiado detalladamente por voltametria cíclica y mediante valoraciones espectrofotométricas ácido-base. El tratamiento matemático de los datos así obtenidos nos ha permitido determinar los valores de pKa para los distintos equilibrios de protonación de los complejos en los estados de oxidación II y III. El complejo out-aqua ha demostrado ser un buen catalizador para la oxidación electroquímica del alcohol benzílico, presumiblemente a benzaldehido. La constante de velocidad de segundo orden para el proceso ha sido determinada como 17.1 M-1 s-1, por simulación matemática. El dímero con un puente cloro, [Ru2Cl(bpp)(trpy)2]2+ ha sido preparado por dos rutas sintéticas diferentes. El dímero análogo con un puente acetato se ha obtenido por reacción del cloro dímero con un exceso de acetato sódico. El dímero con dos ligandos aqua [Ru2(bpp)(trpy)2(OH2)2]3+ puede obtenerse por hidrólisis ácida del complejo con un acetato puente o por hidrólisis básica del complejo con un puente cloro. Estos complejos han sido caracterizados por técnicas estructurales, espectroscópicas y electroquímicas. Las soluciones del dímero con dos ligandos aqua en medio ácido resultan inestables a la coordinación de aniones de la solución con el tiempo. Las propiedades ácido-base del dímero con dos aguas coordinadas han sido estudiadas por voltametria cíclica y mediante experimentos de electrólisis a potencial controlado. El pKa para la desprotonación de uno de los ligandos aqua ha sido determinado mediante una valoración espectrofotométrica ácido-base como 6.7. Este valor tan bajo de pKa se atribuye a la formación de la entidad {Ru2O2H3}, favorable termodinámicamente. Los espectros UV-vis para los distintos estados de oxidación del aqua dímero, de RuIIRuII a RuIIIRuIV, han sido obtenidos por oxidación química y electroquímica del complejo. Se han llevado a cabo estudios cinéticos de la oxidación, paso a paso, de RuII,II a RuIV,IV , y se han determinado las constantes de oxidación de segundo orden para los distintos procesos de oxidación. La capacidad del aqua dímero en la oxidación del agua a oxígeno molecular ha sido investigada en solución homogénea utilizando CeIV como oxidante. La evolución de oxígeno se ha demostrado por cromatografia de gases. Se ha obtenido una eficiencia del 73% y 18.6 ciclos catalíticos, cuando 1.83 x 10-6 moles de dímero se han mezclado con un exceso de 100 equivalentes de cerio. El dímero con dos aguas cataliza también la oxidación del agua de forma heterogénea, con el complejo adsorbido sobre una membrana de nafion, aunque la eficiencia es menor. Se ha propuesto un mecanismo intramolecular para la reacción de oxidación del agua. Consiste en la oxidación a 4 electrones del dímero, de RuII,II a RuIV,IV, el cual reacciona con el agua para formar oxígeno y revierte nuevamente al estado de oxidación II,II. Este modelo es consistente con estudios catalíticos de la evolución de oxígeno en función de las concentraciones de cerio y catalizador, llevados a cabo en solución ácida homogénea, que demuestran que la oxidación a 4 electrones del agua se encuentra catalizada por una sola molécula de complejo bajo concentraciones elevadas de cerio. La constante de pseudo-primer-orden para la evolución de oxígeno tiene un valor de 1.4 x 10-2 s-1, que es uno de los valores de constante más elevados obtenidos hasta la fecha. Desafortunadamente, el aqua dímero se desactiva durante el proceso de catálisis dando lugar a una especie naranja, la cual estamos actualmente tratando de caracterizar.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Can the frequency of blocking be described by a red noise process?

The frequency of persistent atmospheric blocking events in the 40-yr ECMWF Re-Analysis (ERA-40) is compared with the blocking frequency produced by a simple first-order Markov model designed to predict the time evolution of a blocking index [defined by the meridional contrast of potential temperature on the 2-PVU surface (1 PVU ≡ 1 × 10−6 K m2 kg−1 s−1)]. With the observed spatial coherence built into the model, it is able to reproduce the main regions of blocking occurrence and the frequencies of sector blocking very well. This underlines the importance of the climatological background flow in determining the locations of high blocking occurrence as being the regions where the mean midlatitude meridional potential vorticity (PV) gradient is weak. However, when only persistent blocking episodes are considered, the model is unable to simulate the observed frequencies. It is proposed that this persistence beyond that given by a red noise model is due to the self-sustaining nature of the blocking phenomenon.

Non-invasive imaging of roots with high resolution X-ray micro-tomography

ray micro-tomography is a well-established technique for non-invasive imaging and evaluation of heterogeneous materials. An inexpensive X-ray micro-tomography system has been designed and built for the specific purposes of examining root growth and root/soil interactions. The system uses a silver target X-ray source with a focal spot diameter of 80 mum, an X-ray image intensifier with a sampling aperture of about 100 mum, and a sample with a diameter of 25 mm. Pre-germinated wheat and rape seeds were grown for up to 8-10 days in plastic containers in a sandy loam soil sieved to < 250 μm, and imaged with the X-ray system at regular intervals. The quality of 3 D image obtained was good allowing the development and growth of both root axes and some first-order laterals to be observed. The satisfactory discrimination between soil and roots enabled measurements of root diameter (wheat values were 0.48-1.22 mm) in individual tomographic slices and, by tracking from slice to slice, root lengths were also measured. The measurements obtained were generally within 10% of those obtained from destructive samples measured manually and with a flat-bed scanner. Further developments of the system will allow more detailed examination of the root: soil interface.

The acoustic signature of fluid flow in complex porous media

Effective medium approximations for the frequency-dependent and complex-valued effective stiffness tensors of cracked/ porous rocks with multiple solid constituents are developed on the basis of the T-matrix approach (based on integral equation methods for quasi-static composites), the elastic - viscoelastic correspondence principle, and a unified treatment of the local and global flow mechanisms, which is consistent with the principle of fluid mass conservation. The main advantage of using the T-matrix approach, rather than the first-order approach of Eshelby or the second-order approach of Hudson, is that it produces physically plausible results even when the volume concentrations of inclusions or cavities are no longer small. The new formulae, which operates with an arbitrary homogeneous (anisotropic) reference medium and contains terms of all order in the volume concentrations of solid particles and communicating cavities, take explicitly account of inclusion shape and spatial distribution independently. We show analytically that an expansion of the T-matrix formulae to first order in the volume concentration of cavities (in agreement with the dilute estimate of Eshelby) has the correct dependence on the properties of the saturating fluid, in the sense that it is consistent with the Brown-Korringa relation, when the frequency is sufficiently low. We present numerical results for the (anisotropic) effective viscoelastic properties of a cracked permeable medium with finite storage porosity, indicating that the complete T-matrix formulae (including the higher-order terms) are generally consistent with the Brown-Korringa relation, at least if we assume the spatial distribution of cavities to be the same for all cavity pairs. We have found an efficient way to treat statistical correlations in the shapes and orientations of the communicating cavities, and also obtained a reasonable match between theoretical predictions (based on a dual porosity model for quartz-clay mixtures, involving relatively flat clay-related pores and more rounded quartz-related pores) and laboratory results for the ultrasonic velocity and attenuation spectra of a suite of typical reservoir rocks. (C) 2003 Elsevier B.V. All rights reserved.

Modelling impacts of pollution in river systems: a new dispersion model and a case study of mine discharges in the Abrud, Aries and Mures river system in Transylvania, Romania

A new model of dispersion has been developed to simulate the impact of pollutant discharges on river systems. The model accounts for the main dispersion processes operating in rivers as well as the dilution from incoming tributaries and first-order kinetic decay processes. The model is dynamic and simulates the hourly behaviour of river flow and pollutants along river systems. The model has been applied to the Aries and Mures River System in Romania and has been used to assess the impacts of potential dam releases from the Roia Montan Mine in Transylvania, Romania. The question of mine water release is investigated under a range of scenarios. The impacts on pollution levels downstream at key sites and at the border with Hungary are investigated.

Monte Carlo phase diagram for diblock copolymer melts

A partial phase diagram is constructed for diblock copolymer melts using lattice-based Monte Carlo simulations. This is done by locating the order-disorder transition (ODT) with the aid of a recently proposed order parameter and identifying the ordered phase over a wide range of copolymer compositions (0.2 <= f <= 0.8). Consistent with experiments, the disordered phase is found to exhibit direct first-order transitions to each of the ordered morphologies. This includes the spontaneous formation of a perforated-lamellar phase, which presumably forms in place of the gyroid morphology due to finite-size and/or nonequilibrium effects. Also included in our study is a detailed examination of disordered cylinder-forming (f=0.3) diblock copolymers, revealing a substantial degree of pretransitional chain stretching and short-range order that set in well before the ODT, as observed previously in analogous studies on lamellar-forming (f=0.5) molecules. (c) 2006 American Institute of Physics.

The influence of hilly terrain on canopy-atmosphere carbon dioxide exchange

Topography influences many aspects of forest-atmosphere carbon exchange; yet only a small number of studies have considered the role of topography on the structure of turbulence within and above vegetation and its effect on canopy photosynthesis and the measurement of net ecosystem exchange of CO2 (N-ee) using flux towers. Here, we focus on the interplay between radiative transfer, flow dynamics for neutral stratification, and ecophysiological controls on CO2 sources and sinks within a canopy on a gentle cosine hill. We examine how topography alters the forest-atmosphere CO2 exchange rate when compared to uniform flat terrain using a newly developed first-order closure model that explicitly accounts for the flow dynamics, radiative transfer, and nonlinear eco physiological processes within a plant canopy. We show that variation in radiation and airflow due to topography causes only a minor departure in horizontally averaged and vertically integrated photosynthesis from their flat terrain values. However, topography perturbs the airflow and concentration fields in and above plant canopies, leading to significant horizontal and vertical advection of CO2. Advection terms in the conservation equation may be neglected in flow over homogeneous, flat terrain, and then N-ee = F-c, the vertical turbulent flux of CO2. Model results suggest that vertical and horizontal advection terms are generally of opposite sign and of the same order as the biological sources and sinks. We show that, close to the hilltop, F-c departs by a factor of three compared to its flat terrain counterpart and that the horizontally averaged F-c-at canopy top differs by more than 20% compared to the flat-terrain case.