Additional compact formulas for vibrational dynamic dipole polarizabilities and hyperpolarizabilities

Compact expressions, complete through second order in electrical and/or mechanical anharmonicity, are given for the dynamic dipole vibrational polarizability and dynamic first and second vibrational hyperpolarizabilities. Certain contributions not previously formulated are now included

Calculation of Franck-Condon factors including anharmonicity: simulation of the C2 H4+ X̃2B 3u←C2 H4X̃1 Ag band in the photoelectron spectrum of ethylene

Our new simple method for calculating accurate Franck-Condon factors including nondiagonal (i.e., mode-mode) anharmonic coupling is used to simulate the C2H4+X2B 3u←C2H4X̃1 Ag band in the photoelectron spectrum. An improved vibrational basis set truncation algorithm, which permits very efficient computations, is employed. Because the torsional mode is highly anharmonic it is separated from the other modes and treated exactly. All other modes are treated through the second-order perturbation theory. The perturbation-theory corrections are significant and lead to a good agreement with experiment, although the separability assumption for torsion causes the C2 D4 results to be not as good as those for C2 H4. A variational formulation to overcome this circumstance, and deal with large anharmonicities in general, is suggested

The effect of counterpoise correction and relaxation energy term to the internal rotation barriers: application to the BF3 ...NH3 and C2H4 ...SO2 dimers

The relevance of the fragment relaxation energy term and the effect of the basis set superposition error on the geometry of the BF3⋯NH3 and C2H4⋯SO2 van der Waals dimers have been analyzed. Second-order Møller-Plesset perturbation theory calculations with the d95(d,p) basis set have been used to calculate the counterpoise-corrected barrier height for the internal rotations. These barriers have been obtained by relocating the stationary points on the counterpoise-corrected potential energy surface of the processes involved. The fragment relaxation energy can have a large influence on both the intermolecular parameters and barrier height. The counterpoise correction has proved to be important for these systems

C-H···O H-bonded complexes: how does basis set superposition error change their potential-energy surfaces?

Geometries, vibrational frequencies, and interaction energies of the CNH⋯O3 and HCCH⋯O3 complexes are calculated in a counterpoise-corrected (CP-corrected) potential-energy surface (PES) that corrects for the basis set superposition error (BSSE). Ab initio calculations are performed at the Hartree-Fock (HF) and second-order Møller-Plesset (MP2) levels, using the 6-31G(d,p) and D95++(d,p) basis sets. Interaction energies are presented including corrections for zero-point vibrational energy (ZPVE) and thermal correction to enthalpy at 298 K. The CP-corrected and conventional PES are compared; the unconnected PES obtained using the larger basis set including diffuse functions exhibits a double well shape, whereas use of the 6-31G(d,p) basis set leads to a flat single-well profile. The CP-corrected PES has always a multiple-well shape. In particular, it is shown that the CP-corrected PES using the smaller basis set is qualitatively analogous to that obtained with the larger basis sets, so the CP method becomes useful to correctly describe large systems, where the use of small basis sets may be necessary

How does basis set superposition error change the potential surfaces for hydrogen-bonded dimers?

We describe a simple method to automate the geometric optimization of molecular orbital calculations of supermolecules on potential surfaces that are corrected for basis set superposition error using the counterpoise (CP) method. This method is applied to the H-bonding complexes HF/HCN, HF/H2O, and HCCH/H2O using the 6-31G(d,p) and D95 + + (d,p) basis sets at both the Hartree-Fock and second-order Møller-Plesset levels. We report the interaction energies, geometries, and vibrational frequencies of these complexes on the CP-optimized surfaces; and compare them with similar values calculated using traditional methods, including the (more traditional) single point CP correction. Upon optimization on the CP-corrected surface, the interaction energies become more negative (before vibrational corrections) and the H-bonding stretching vibrations decrease in all cases. The extent of the effects vary from extremely small to quite large depending on the complex and the calculational method. The relative magnitudes of the vibrational corrections cannot be predicted from the H-bond stretching frequencies alone

A comparative analysis by means of quantum molecular similarity measures of density distributions derived from conventional ab initio and density functional methods

A procedure based on quantum molecular similarity measures (QMSM) has been used to compare electron densities obtained from conventional ab initio and density functional methodologies at their respective optimized geometries. This method has been applied to a series of small molecules which have experimentally known properties and molecular bonds of diverse degrees of ionicity and covalency. Results show that in most cases the electron densities obtained from density functional methodologies are of a similar quality than post-Hartree-Fock generalized densities. For molecules where Hartree-Fock methodology yields erroneous results, the density functional methodology is shown to yield usually more accurate densities than those provided by the second order Møller-Plesset perturbation theory

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Relacions entre la identitat catalana i la percepció de vitalitat etnolingüística en una mostra d’estudiants universitaris

This article analyses, on the one hand, the perception of ethno-linguistic vitality in a sample of university students in the province of Girona, and, on the other, links this perceived vitality to the feeling of identity in the group to which they belong. 112 Catalan-speaking students who described themselves as Catalan (Catalan identity) responded to the questionnaire on subjective vitality, in terms of adaptation to Catalan. This instrument assesses the beliefs about one’s own group’s and its language’s vitality (in this case Catalan) and that of the exocentric group and language (in this case Spanish). The results show that the students see their language (Catalan) as having greater status and institutional support than Spanish, but a lower demographic impact (birth rate, immigration, mixed marriages, emigration). Identity is also positively correlated with the perception of ethno-linguistic perception. These results are discussed in relation to the role of vitality in the choice of the language for education (Catalan versus Spanish) and day-to-day communication, and the level at which language policy is defined in order to promote communication in Catalan in the fields of business, politics, education and the media

Aprofitament de residus vegetals per a la concentració i separació d'ions metàl·lics de solucions aquoses

La adsorció sobre carbó actiu és un tractament habitual para a l'eliminació de metalls contaminants d'afluents aquosos. Amb la finalitat de reduir el cost dels tractament actuals s'està estudiant la utilització d'adsorbents de baix cost com a substituts dels sorbents convencionals que s'utilitzen usualment en els tractaments d'aigües. En aquesta tesi s'ha avaluat la possibilitat d'aprofitar diversos residus vegetals procedents de processos agrícoles o residus industrials, com a sorbents per a l'eliminació de metalls de solucions aquoses. En concret, el treball s'ha centrat en assolir dos grans objectius, primer de tot, conèixer la possibilitat d'utilitzar el pinyol d'oliva com a adsorbent de metalls divalents de solucions aquoses, com són el Cd(II), Cu(II), Ni(II) i Pb(II), i en segon lloc, avaluar la utilització de quatre residus vegetals (escorça de yohimbe, rapa de raïm, suro i pinyol d'oliva) per a l'extracció de Cr(VI) de solucions aquoses i utilitzar el residu vegetal més apropiat com a base per tal d'obtenir un nou sorbent, amb més bona morfologia per ser utilitzat en sistemes d'eliminació de Cr(VI) en continu. Les dades cinètiques d'adsorció es van tractar amb el model cinètic de pseudo-segon ordre, amb el que es van determinar la velocitat inicial d'adsorció dels metalls. Els models d'isotermes d'equilibri de Langmuir i Freundlich es van utilitzar per descriure els equilibris d'adsorció i calcular els paràmetres de les isotermes. En els assajos amb mescles binàries de metalls, els resultats es van tractar amb el model ampliat de Langmuir. En aquest treball s'ha determinat, primerament, que el pinyol d'oliva és un bon sorbent dels metalls divalents estudiats. S'ha determinat que l'adsorció dels metalls varia amb el pH de la solució i que la presència de sals en la solució provoca una davallada en l'adsorció. També es va determinar la competència en l'adsorció entre els metalls al realitzar assajos en mescles binàries. Amb els estudis d'extracció de Cr(VI) amb els quatre materials s'ha determinat que aquests vegetals poden adsorbir el Cr(VI) en solució, tot i que la rapa i la yohimbe són els materials que presenten unes capacitats d'adsorció més elevades. Les anàlisis amb espectroscòpia electrònica de Raig-X realitzats en la rapa de raïm i la yohimbe en contacte amb Cr(VI) varen demostrar la presència del metall en els dos estats d'oxidació, Cr(VI) i Cr(III), evidenciant la participació del mecanisme de reducció de Cr(VI) en l'adsorció del metall sobre els materials. Amb l'encapsulació de la rapa de raïm es va incrementar significativament la capacitat d'adsorció del metall respecta a la seva utilització en la forma original. Les dades experimentals es van tractar amb un model de difusió i es va determinar el coeficient de difusió en l'adsorció en el sorbent amb diferent percentatge de rapa en l'interior. Finalment s'han realitzat estudis preliminars d'adsorció de crom hexavalent en continu utilitzant columnes de petites dimensions farcides amb el sorbent format per boles de rapa encapsulada en alginat càlcic. En aquests assajos s'ha determinat l'efecte en l'adsorció de metall de diversos paràmetres d'operació, com la concentració inicial de metall, el cabal i l'alçada de llit. El model de transferència de matèria i difusió en columna emprat pel tractament de les dades experimentals va mostrar l'efecte de la concentració inicial i del cabal en els coeficients de difusió i transferència de matèria. El model BDST no va predir correctament el temps de servei en els diferents paràmetres d'operació experimentats a causa de la complexitat del procés d'adsorció. Amb aquests assajos preliminars s'han establert les bases de l'adsorció de Cr(VI) amb el nou sorbent en sistemes en continu amb la finalitat de realitzar estudis posteriors d'eliminació del metall en columnes de majors dimensions.

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.

Biblioterapia : a Função Terapêutica do Livro em Ambiente Prisional

Este trabalho de investigação consiste num estudo empírico sobre uma nova ciência - a Biblioterapia, e a situação mais específica, a sua aplicação em contexto prisional. Para alcançar este objetivo formulámos a seguinte hipótese: “o livro exerce sobre as pessoas em geral e, de uma forma mais particular, sobre indivíduos em contexto de detenção institucional uma função terapêutica”. A metodologia adotada para esta investigação assenta essencialmente em fontes documentais textuais e digitais, sobre autores especialistas na matéria, que de uma forma mais precisa se identificam e citam especificamente no capítulo referente à revisão da literatura. Para atingir o nosso propósito começámos por definir os conceitos de leitura, de leitura terapêutica, de fenomenologia da linguagem, de terapia e de diálogo assim como procurámos saber quais os fundamentos filosóficos e os componentes biblioterapêuticos, público-alvo e áreas com maiores potencialidades de aplicação, benefícios e limitações que norteiam à aplicação da Biblioterapia, nomeadamente, para fins de desenvolvimento pessoal e de alterações comportamentais. Fez-se igualmente uma análise da evolução do conceito, desde Aristóteles até aos nossos dias. Só com o conhecimento efetivo desta realidade seria possível prosseguir o nosso objeto de estudo. Entendemos ainda ser importante perceber qual é a situação da leitura em Portugal e as políticas desenvolvidas nesta área pelo Estado, para promover a leitura junto dos estudantes e da população em geral. Apenas assente nestas duas realidades que se entrecruzam, no processo de leitor/livro é possível identificar e relacionar esta problemática em estudo. Constatamos a existência de dois tipos de Biblioterapia, a Biblioterapia-arte e a Biblioterapia-ciência, identificando algumas das suas características, principalmente o que as aproxima e o que as separa. Apresentamos um exemplo de uma área reconhecidamente bem sucedida na aplicação da Biblioterapia não-clínica, em instituições prisionais, efetuando uma abordagem às questões mais pertinentes neste campo: os seus aspetos institucionais, os psicossociais e os de reinserção social. Considerámos ainda as recomendações da IFLA para este setor populacional – Guidelines for Library Services to Prisoners. Ao longo do trabalho e de uma forma transversal realçamos as interrelações biblioterapêuticas entre pacientes, bibliotecários e outros intervenientes no processo biblioterapêutico, tais como médicos e orientadores de leitura. Finalizamos o trabalho concluindo que a prática de leitura de temas específicos, em qualquer campo de atuação e sob orientação de um elemento profissional com conhecimentos aprofundados tanto em relação às personalidades e problemas dos destinatários como aos materiais que tem à sua disposição – o livro, em sentido lato - pode efetivamente desencadear benefícios terapêuticos nos leitores.

Approximate Gauss-Newton methods for nonlinear least squares problems

The Gauss–Newton algorithm is an iterative method regularly used for solving nonlinear least squares problems. It is particularly well suited to the treatment of very large scale variational data assimilation problems that arise in atmosphere and ocean forecasting. The procedure consists of a sequence of linear least squares approximations to the nonlinear problem, each of which is solved by an “inner” direct or iterative process. In comparison with Newton’s method and its variants, the algorithm is attractive because it does not require the evaluation of second-order derivatives in the Hessian of the objective function. In practice the exact Gauss–Newton method is too expensive to apply operationally in meteorological forecasting, and various approximations are made in order to reduce computational costs and to solve the problems in real time. Here we investigate the effects on the convergence of the Gauss–Newton method of two types of approximation used commonly in data assimilation. First, we examine “truncated” Gauss–Newton methods where the inner linear least squares problem is not solved exactly, and second, we examine “perturbed” Gauss–Newton methods where the true linearized inner problem is approximated by a simplified, or perturbed, linear least squares problem. We give conditions ensuring that the truncated and perturbed Gauss–Newton methods converge and also derive rates of convergence for the iterations. The results are illustrated by a simple numerical example. A practical application to the problem of data assimilation in a typical meteorological system is presented.

Relationships between lightning and properties of convective cloud clusters

Satellite observations of convective system properties and lightning flash rate are used to investigate the ability of potential lightning parameterizations to capture both the dominant land-ocean contrast in lightning occurrence and regional differences between Africa, the Amazon and the islands of the maritime continent. As found in previous studies, the radar storm height is tightly correlated with the lightning flash rate. A roughly second order power-law fit to the mean radar echo top height above the 0C isotherm is shown to capture both regional and land-ocean contrasts in lightning occurrence and flash rate using a single set of parameters. Recent developments should soon make it possible to implement a parameterization of this kind in global models. Parameterizations based on cloud top height, convective rain rate and convective rain fraction all require the use of separate fits over land and ocean and fail to capture observed differences between continental regions.

Structure of turbulent flow over arrays of cubical roughness

The structure of turbulent flow over large roughness consisting of regular arrays of cubical obstacles is investigated numerically under constant pressure gradient conditions. Results are analysed in terms of first- and second-order statistics, by visualization of instantaneous flow fields and by conditional averaging. The accuracy of the simulations is established by detailed comparisons of first- and second-order statistics with wind-tunnel measurements. Coherent structures in the log region are investigated. Structure angles are computed from two-point correlations, and quadrant analysis is performed to determine the relative importance of Q2 and Q4 events (ejections and sweeps) as a function of height above the roughness. Flow visualization shows the existence of low-momentum regions (LMRs) as well as vortical structures throughout the log layer. Filtering techniques are used to reveal instantaneous examples of the association of the vortices with the LMRs, and linear stochastic estimation and conditional averaging are employed to deduce their statistical properties. The conditional averaging results reveal the presence of LMRs and regions of Q2 and Q4 events that appear to be associated with hairpin-like vortices, but a quantitative correspondence between the sizes of the vortices and those of the LMRs is difficult to establish; a simple estimate of the ratio of the vortex width to the LMR width gives a value that is several times larger than the corresponding ratio over smooth walls. The shape and inclination of the vortices and their spatial organization are compared to recent findings over smooth walls. Characteristic length scales are shown to scale linearly with height in the log region. Whilst there are striking qualitative similarities with smooth walls, there are also important differences in detail regarding: (i) structure angles and sizes and their dependence on distance from the rough surface; (ii) the flow structure close to the roughness; (iii) the roles of inflows into and outflows from cavities within the roughness; (iv) larger vortices on the rough wall compared to the smooth wall; (v) the effect of the different generation mechanism at the wall in setting the scales of structures.

Application of the Lighthill-Ford theory of spontaneous imbalance to clear-air turbulence forecasting

A new method of clear-air turbulence (CAT) forecasting based on the Lighthill–Ford theory of spontaneous imbalance and emission of inertia–gravity waves has been derived and applied on episodic and seasonal time scales. A scale analysis of this shallow-water theory for midlatitude synoptic-scale flows identifies advection of relative vorticity as the leading-order source term. Examination of leading- and second-order terms elucidates previous, more empirically inspired CAT forecast diagnostics. Application of the Lighthill–Ford theory to the Upper Mississippi and Ohio Valleys CAT outbreak of 9 March 2006 results in good agreement with pilot reports of turbulence. Application of Lighthill–Ford theory to CAT forecasting for the 3 November 2005–26 March 2006 period using 1-h forecasts of the Rapid Update Cycle (RUC) 2 1500 UTC model run leads to superior forecasts compared to the current operational version of the Graphical Turbulence Guidance (GTG1) algorithm, the most skillful operational CAT forecasting method in existence. The results suggest that major improvements in CAT forecasting could result if the methods presented herein become operational.