Monte Carlo thermodynamic and structural properties of the TIP4P water model: dependence on the computational conditions

Classical Monte Carlo simulations were carried out on the NPT ensemble at 25°C and 1 atm, aiming to investigate the ability of the TIP4P water model [Jorgensen, Chandrasekhar, Madura, Impey and Klein; J. Chem. Phys., 79 (1983) 926] to reproduce the newest structural picture of liquid water. The results were compared with recent neutron diffraction data [Soper; Bruni and Ricci; J. Chem. Phys., 106 (1997) 247]. The influence of the computational conditions on the thermodynamic and structural results obtained with this model was also analyzed. The findings were compared with the original ones from Jorgensen et al [above-cited reference plus Mol. Phys., 56 (1985) 1381]. It is notice that the thermodynamic results are dependent on the boundary conditions used, whereas the usual radial distribution functions g(O/O(r)) and g(O/H(r)) do not depend on them.

Termodinâmica estatística de líquidos com o método de Monte Carlo. II. O líquido clorofórmio

Thermodynamic properties and radial distribution functions for liquid chloroform were calculated using the Monte Carlo method implemented with Metropolis algorithm in the NpT ensemble at 298 K and 1 atm. A five site model was developed to represent the chloroform molecules. A force field composed by Lennard-Jones and Coulomb potential functions was used to calculate the intermolecular energy. The partial charges needed to represent the Coulombic interactions were obtained from quantum chemical ab initio calculations. The Lennard-Jones parameters were adjusted to reproduce experimental values for density and enthalpy of vaporization for pure liquid. All thermodynamic results are in excelent agreement with experimental data. The correlation functions calculated are in good accordance with theoretical results avaliable in the literature. The free energy for solvating one chloroform molecule into its own liquid at 298 K and 1 atm was computed as an additional test of the potential model. The result obtained compares well with the experimental value. The medium effects on cis/trans convertion of a hypotetical solute in water TIP4P and chloroform solvents were also accomplished. The results obtained from this investigation are in agreement with estimates of the continuous theory of solvation.

Redes neurais e suas aplicações em calibração multivariada

Neural Networks are a set of mathematical methods and computer programs designed to simulate the information process and the knowledge acquisition of the human brain. In last years its application in chemistry is increasing significantly, due the special characteristics for model complex systems. The basic principles of two types of neural networks, the multi-layer perceptrons and radial basis functions, are introduced, as well as, a pruning approach to architecture optimization. Two analytical applications based on near infrared spectroscopy are presented, the first one for determination of nitrogen content in wheat leaves using multi-layer perceptrons networks and second one for determination of BRIX in sugar cane juices using radial basis functions networks.

Fragmentation of massive dense cores down to <~ 1000 AU: relation between fragmentation and density structure

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 μm (or 1.2 mm in two cases) and the spectral energy distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. Even though the estimated fragmentation level is strictly speaking a lower limit, its relative value is significant and several trends could be explored with our data. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

3D kinematics of the near-IR HH 223 outflow in L723

In this work, we derive the full 3D kinematics of the near-infrared outflow HH 223, located in the dark cloud Lynds 723 (L723), where a well-defined quadrupolar CO outflow is found. HH 223 appears projected on to the two lobes of the eastwest CO outflow. The radio continuum source VLA 2, towards the centre of the CO outflow, harbours a multiple system of low-mass young stellar objects. One of the components has been proposed to be the exciting source of the eastwest CO outflow. From the analysis of the kinematics, we get further evidence on the relationship between the near-infrared and CO outflows and on the location of their exciting source. The proper motions were derived using multi-epoch, narrow-band H2 (2.122 μm line) images. Radial velocities were derived from the 2.122 μm line of the spectra. Because of the extended (∼5 arcmin), S-shaped morphology of the target, the spectra were obtained with the multi-object-spectroscopy (MOS) observing mode using the instrument Long-Slit Intermediate Resolution Infrared Spectrograph (LIRIS) at the 4.2 m William Herschel Telescope. To our knowledge, this work is the first time that MOS observing mode has been successfully used in the near-infrared range for an extended target.

Aplicação de alguns modelos quimiométricos à espectroscopia de fluorescência de raios-X de energia dispersiva

The objective of this work was to accomplish the simultaneous determination of some chemical elements by Energy Dispersive X-ray Fluorescence (EDXRF) Spectroscopy through multivariate calibration in several sample types. The multivariate calibration models were: Back Propagation neural network, Levemberg-Marquardt neural network and Radial Basis Function neural network, fuzzy modeling and Partial Least Squares Regression. The samples were soil standards, plant standards, and mixtures of lead and sulfur salts diluted in silica. The smallest Root Mean Square errors (RMS) were obtained with Back Propagation neural networks, which solved main EDXRF problems in a better way.

SSPBE: um programa para solução numérica da equação de Poisson-Boltzmann em simetria esférica com modelo de adsorção

A Fortran77 program, SSPBE, designed to solve the spherically symmetric Poisson-Boltzmann equation using cell model for ionic macromolecular aggregates or macroions is presented. The program includes an adsorption model for ions at the aggregate surface. The working algorithm solves the Poisson-Boltzmann equation in the integral representation using the Picard iteration method. Input parameters are introduced via an ASCII file, sspbe.txt. Output files yield the radial distances versus mean field potentials and average molar ion concentrations, the molar concentration of ions at the cell boundary, the self-consistent degree of ion adsorption from the surface and other related data. Ion binding to ionic, zwitterionic and reverse micelles are presented as representative examples of the applications of the SSPBE program.

Radiaalikompressorin juoksupyörän virtauskanavan suunnittelu

Diplomityössä on tutkittu radiaalikompressorin juoksupyörän virtauskanavan suunnittelua ja esitelty yleisesti käytetyt suunnittelumenetelmät. Virtauskanavan muodon vaikutukset virtauskenttään on esitetty. Työssä on kehitetty uusi tietokoneohjelma juoksupyörän kolmiulotteisen geometrian luomiseen, yksiulotteisesta suunnittelusta saatujen päämittojen pohjalta. Ohjelmalla luotua juoksupyörän geometriaa voidaan käyttää virtaus- ja lujuusanalyyseissä. Ohjelmaan on myös liitetty kaksiulotteinen kokoonpuristuvan viskoosittoman virtauksen ratkaisija.

Estudo da estrutura e da solvatação do HNP-3, um antibiótico natural, por dinâmica molecular

The structure and hydration of the HNP-3 have been derived from molecular dynamics data using root mean square deviation, radial and energy distributions. Three antiparallel beta sheets were found to be preserved. 15 intramolecular hydrogen bonds were identified together with 36 hydrogen bonds on the backbone and 35 on the side chain atoms. From the point of view of the hydration dynamics, the analysis shows a high solvent accessibility of the monomer and attractive interactions with water molecules.

Functional diversity of macrophyte communities within and between Pyrenean lakes

Submersed vegetation is a common feature in about 70% Pyrene an high mountain (>1500 m a.s.l.) lakes. Isoetids and soft-water elodeids are common elements of this underw ater flora and can form distinct vegetation units (i.e. patches of vegetation dominated by different species) within complex mosaics of vegetation in shallow waters (<7 m). Since is oetids exert a strong influence on sediment biogeochemistry due to high radial oxygen loss, we examined the small scale characteristics of the lake environment (water and sediment) associated to vegetation patches in order to ascertain potential functional differences among them. To do so, we characterised the species composition and biomass of the main vegetation units from 11 lakes, defined plant communities based on biomass data, and then related each community with sediment properties (redox and dissolved nutrient concentration in the pore water) and water nutrient concentration within plant canopy. We also characterised lake water and sediment in areas without vegetation as a reference. A total of twenty-one vegetation units were identified, ranging from one to five per lake. A cluster analysis on biomass species composition suggested seven different macrophyte communities that were named after the most dominant species:Nitella sp.,Potamogeton praelongus, Myriophyllum alterniflorum, Sparganium angustifolium , Isoetes echinospora,Isoetes lacustris and Carex rostrata . Coupling between macrophyte communities and their immediate environment (overlying water and sediment) was manifested mainly as variation in sediment redox conditions and the dominant form of inorganic nitrogen in pore-water. These effects depended on the specific compositi on of the community, and on the allocation between above- and belowground biomass, and could be predicted with a model relating the average and standard deviation of sediment redox potential from 0 down to -20 cm, across macrophyte communities. Differences in pore-water total dissolved phosphorus were related to the trophic state of the lakes. There was no correlation between sediment and water column dissolved nutrients. However, nitrate concentrations tended to be lower in the water overlaying isoetid communities, in apparent contradiction to the patterns of dissolved nitrates in the pore-water. These tendencies were robust even when comparing the water over laying communities within the same lake, thus pointing towards a potential effect of isoetids in reducing dissolved nitrogen in the lakes.

El fotorreactor elíptico

In this paper, the mathematical model of the elliptical photoreactor , an special type of reactor that uses ultraviolet radiation, is presented. In the elliptical photoreactor the cylindrical reactor is irradiated from the outside by placing the lamp and the reactor at the foci of an elliptical reflector. The two main models of radiation -radial and difusse- are studied, an finally the general method of resolution of the mathematical model and its resolution in certain simple cases is shown.

A Practical Micro Mechanical Model Based Local Approach Methodology for the Analysis of Ductile Fracture of Welded T-Joints

This thesis concentrates on developing a practical local approach methodology based on micro mechanical models for the analysis of ductile fracture of welded joints. Two major problems involved in the local approach, namely the dilational constitutive relation reflecting the softening behaviour of material, and the failure criterion associated with the constitutive equation, have been studied in detail. Firstly, considerable efforts were made on the numerical integration and computer implementation for the non trivial dilational Gurson Tvergaard model. Considering the weaknesses of the widely used Euler forward integration algorithms, a family of generalized mid point algorithms is proposed for the Gurson Tvergaard model. Correspondingly, based on the decomposition of stresses into hydrostatic and deviatoric parts, an explicit seven parameter expression for the consistent tangent moduli of the algorithms is presented. This explicit formula avoids any matrix inversion during numerical iteration and thus greatly facilitates the computer implementation of the algorithms and increase the efficiency of the code. The accuracy of the proposed algorithms and other conventional algorithms has been assessed in a systematic manner in order to highlight the best algorithm for this study. The accurate and efficient performance of present finite element implementation of the proposed algorithms has been demonstrated by various numerical examples. It has been found that the true mid point algorithm (a = 0.5) is the most accurate one when the deviatoric strain increment is radial to the yield surface and it is very important to use the consistent tangent moduli in the Newton iteration procedure. Secondly, an assessment of the consistency of current local failure criteria for ductile fracture, the critical void growth criterion, the constant critical void volume fraction criterion and Thomason's plastic limit load failure criterion, has been made. Significant differences in the predictions of ductility by the three criteria were found. By assuming the void grows spherically and using the void volume fraction from the Gurson Tvergaard model to calculate the current void matrix geometry, Thomason's failure criterion has been modified and a new failure criterion for the Gurson Tvergaard model is presented. Comparison with Koplik and Needleman's finite element results shows that the new failure criterion is fairly accurate indeed. A novel feature of the new failure criterion is that a mechanism for void coalescence is incorporated into the constitutive model. Hence the material failure is a natural result of the development of macroscopic plastic flow and the microscopic internal necking mechanism. By the new failure criterion, the critical void volume fraction is not a material constant and the initial void volume fraction and/or void nucleation parameters essentially control the material failure. This feature is very desirable and makes the numerical calibration of void nucleation parameters(s) possible and physically sound. Thirdly, a local approach methodology based on the above two major contributions has been built up in ABAQUS via the user material subroutine UMAT and applied to welded T joints. By using the void nucleation parameters calibrated from simple smooth and notched specimens, it was found that the fracture behaviour of the welded T joints can be well predicted using present methodology. This application has shown how the damage parameters of both base material and heat affected zone (HAZ) material can be obtained in a step by step manner and how useful and capable the local approach methodology is in the analysis of fracture behaviour and crack development as well as structural integrity assessment of practical problems where non homogeneous materials are involved. Finally, a procedure for the possible engineering application of the present methodology is suggested and discussed.

Color ad spectral image assessment using novel quality and fidelity techniques

The ongoing development of the digital media has brought a new set of challenges with it. As images containing more than three wavelength bands, often called spectral images, are becoming a more integral part of everyday life, problems in the quality of the RGB reproduction from the spectral images have turned into an important area of research. The notion of image quality is often thought to comprise two distinctive areas – image quality itself and image fidelity, both dealing with similar questions, image quality being the degree of excellence of the image, and image fidelity the measure of the match of the image under study to the original. In this thesis, both image fidelity and image quality are considered, with an emphasis on the influence of color and spectral image features on both. There are very few works dedicated to the quality and fidelity of spectral images. Several novel image fidelity measures were developed in this study, which include kernel similarity measures and 3D-SSIM (structural similarity index). The kernel measures incorporate the polynomial, Gaussian radial basis function (RBF) and sigmoid kernels. The 3D-SSIM is an extension of a traditional gray-scale SSIM measure developed to incorporate spectral data. The novel image quality model presented in this study is based on the assumption that the statistical parameters of the spectra of an image influence the overall appearance. The spectral image quality model comprises three parameters of quality: colorfulness, vividness and naturalness. The quality prediction is done by modeling the preference function expressed in JNDs (just noticeable difference). Both image fidelity measures and the image quality model have proven to be effective in the respective experiments.

Geoestatística como metodologia para estudar a dinâmica espaço-temporal de doenças associadas a Colletotrichum spp. transmitidos por sementes

Objetivou-se no presente trabalho caracterizar o progresso espaço-temporal da severidade da antracnose do feijoeiro comum e da ramulose do algodoeiro por meio da Geoestatística. Os experimentos foram conduzidos no campo, durante o período das águas, em diferentes épocas. Sementes inoculadas pelo método da restrição hídrica foram semeadas no centro de parcelas experimentais constituindo fonte de inóculo do tipo ponto. Foram realizadas semanalmente seis avaliações da severidade das doenças com base em escalas de notas. Pelos modelos de semivariograma isotrópicos esféricos e gaussianos ajustados aos dados, verificou-se o padrão de distribuição agregado e a dependência espacial de ambos os patossistemas. Com o mapeamento da severidade das doenças pelo método da krigagem ordinária em blocos, verificou-se menor severidade das doenças nos primeiros estádios, com aumento gradual ao longo do tempo. A antracnose apresentou formação inicial de focos de inóculo secundário, além do foco com inóculo inicial, que coalesceram com o tempo. A severidade da ramulose diminuiu gradualmente de forma radial e contínua do centro da parcela para as extremidades, com maior capacidade de infecção de plantas vizinhas, quando comparada à antracnose, possivelmente pela maior agressividade do patógeno e hábito de crescimento arbustivo do algodoeiro.