Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.

Thermal error evaluation and modelling of a CNC cylindrical grinding machine

With the relentless quest for improved performance driving ever tighter tolerances for manufacturing, machine tools are sometimes unable to meet the desired requirements. One option to improve the tolerances of machine tools is to compensate for their errors. Among all possible sources of machine tool error, thermally induced errors are, in general for newer machines, the most important. The present work demonstrates the evaluation and modelling of the behaviour of the thermal errors of a CNC cylindrical grinding machine during its warm-up period.

Comparison between detailed model simulation and artificial neural network for forecasting building energy consumption

There are several ways to attempt to model a building and its heat gains from external sources as well as internal ones in order to evaluate a proper operation, audit retrofit actions, and forecast energy consumption. Different techniques, varying from simple regression to models that are based on physical principles, can be used for simulation. A frequent hypothesis for all these models is that the input variables should be based on realistic data when they are available, otherwise the evaluation of energy consumption might be highly under or over estimated. In this paper, a comparison is made between a simple model based on artificial neural network (ANN) and a model that is based on physical principles (EnergyPlus) as an auditing and predicting tool in order to forecast building energy consumption. The Administration Building of the University of Sao Paulo is used as a case study. The building energy consumption profiles are collected as well as the campus meteorological data. Results show that both models are suitable for energy consumption forecast. Additionally, a parametric analysis is carried out for the considered building on EnergyPlus in order to evaluate the influence of several parameters such as the building profile occupation and weather data on such forecasting. (C) 2008 Elsevier B.V. All rights reserved.

Thermal/mechanical simulation and laboratory fatigue testing of an alternative yttria tetragonal zirconia polycrystal core-veneer all-ceramic layered crown design

This study evaluated the stress levels at the core layer and the veneer layer of zirconia crowns (comprising an alternative core design vs. a standard core design) under mechanical/thermal simulation, and subjected simulated models to laboratory mouth-motion fatigue. The dimensions of a mandibular first molar were imported into computer-aided design (CAD) software and a tooth preparation was modeled. A crown was designed using the space between the original tooth and the prepared tooth. The alternative core presented an additional lingual shoulder that lowered the veneer bulk of the cusps. Finite element analyses evaluated the residual maximum principal stresses fields at the core and veneer of both designs under loading and when cooled from 900 degrees C to 25 degrees C. Crowns were fabricated and mouth-motion fatigued, generating master Weibull curves and reliability data. Thermal modeling showed low residual stress fields throughout the bulk of the cusps for both groups. Mechanical simulation depicted a shift in stress levels to the core of the alternative design compared with the standard design. Significantly higher reliability was found for the alternative core. Regardless of the alternative configuration, thermal and mechanical computer simulations showed stress in the alternative core design comparable and higher to that of the standard configuration, respectively. Such a mechanical scenario probably led to the higher reliability of the alternative design under fatigue.

Study of the optical and magnetic properties of pyrimidine in water combining PCM and QM/MM methodologies

The solvent effects on the low-lying absorption spectrum and on the (15)N chemical shielding of pyrimidine in water are calculated using the combined and sequential Monte Carlo simulation and quantum mechanical calculations. Special attention is devoted to the solute polarization. This is included by an iterative procedure previously developed where the solute is electrostatically equilibrated with the solvent. In addition, we verify the simple yet unexplored alternative of combining the polarizable continuum model (PCM) and the hybrid QM/MM method. We use PCM to obtain the average solute polarization and include this in the MM part of the sequential QM/MM methodology, PCM-MM/QM. These procedures are compared and further used in the discrete and the explicit solvent models. The use of the PCM polarization implemented in the MM part seems to generate a very good description of the average solute polarization leading to very good results for the n-pi* excitation energy and the (15)N nuclear chemical shield of pyrimidine in aqueous environment. The best results obtained here using the solute pyrimidine surrounded by 28 explicit water molecules embedded in the electrostatic field of the remaining 472 molecules give the statistically converged values for the low lying n-pi* absorption transition in water of 36 900 +/- 100 (PCM polarization) and 36 950 +/- 100 cm(-1) (iterative polarization), in excellent agreement among one another and with the experimental value observed with a band maximum at 36 900 cm(-1). For the nuclear shielding (15)N the corresponding gas-water chemical shift obtained using the solute pyrimidine surrounded by 9 explicit water molecules embedded in the electrostatic field of the remaining 491 molecules give the statistically converged values of 24.4 +/- 0.8 and 28.5 +/- 0.8 ppm, compared with the inferred experimental value of 19 +/- 2 ppm. Considering the simplicity of the PCM over the iterative polarization this is an important aspect and the computational savings point to the possibility of dealing with larger solute molecules. This PCM-MM/QM approach reconciles the simplicity of the PCM model with the reliability of the combined QM/MM approaches.

Precipitation diurnal cycle and summer climatology assessment over South America: An evaluation of Regional Climate Model version 3 simulations

Regional Climate Model version 3 (RegCM3) simulations of 17 summers (1988-2004) over part of South America south of 5 degrees S were evaluated to identify model systematic errors. Model results were compared to different rainfall data sets (Climate Research Unit (CRU), Climate Prediction Center (CPC), Global Precipitation Climatology Project (GPCP), and National Centers for Environmental Prediction (NCEP) reanalysis), including the five summers mean (1998-2002) precipitation diurnal cycle observed by the Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR). In spite of regional differences, the RegCM3 simulates the main observed aspects of summer climatology associated with the precipitation (northwest-southeast band of South Atlantic Convergence Zone (SACZ)) and air temperature (warmer air in the central part of the continent and colder in eastern Brazil and the Andes Mountains). At a regional scale, the main RegCM3 failures are the underestimation of the precipitation in the northern branch of the SACZ and some unrealistic intense precipitation around the Andes Mountains. However, the RegCM3 seasonal precipitation is closer to the fine-scale analyses (CPC, CRU, and TRMM-PR) than is the NCEP reanalysis, which presents an incorrect north-south orientation of SACZ and an overestimation of its intensity. The precipitation diurnal cycle observed by TRMM-PR shows pronounced contrasts between Tropics and Extratropics and land and ocean, where most of these features are simulated by RegCM3. The major similarities between the simulation and observation, especially the diurnal cycle phase, are found over the continental tropical and subtropical SACZ regions, which present afternoon maximum (1500-1800 UTC) and morning minimum (0900-1200 UTC). More specifically, over the core of SACZ, the phase and amplitude of the simulated precipitation diurnal cycle are very close to the TRMM-PR observations. Although there are amplitude differences, the RegCM3 simulates the observed nighttime rainfall in the eastern Andes Mountains, over the Atlantic Ocean, and also over northern Argentina. The main simulation deficiencies are found in the Atlantic Ocean and near the Andes Mountains. Over the Atlantic Ocean the convective scheme is not triggered; thus the rainfall arises from the grid-scale scheme and therefore differs from the TRMM-PR. Near the Andes, intense (nighttime and daytime) simulated precipitation could be a response of an incorrect circulation and topographic uplift. Finally, it is important to note that unlike most reported bias of global models, RegCM3 does not trigger the moist convection just after sunrise over the southern part of the Amazon.

Dipole polarizability and Rayleigh light scattering by the hydrated electron

Assuming the existence of a confined state of the electron in bulk water the polarizability of the hydrated electron is analyzed. Statistically uncorrelated supermolecular structures composed of seven water molecules (first solvation shell) with an extra electron were extracted from classical Monte Carlo simulation and used in quantum mechanical second-order Moller-Plesset calculations. It is found that the bound excess electron contributes with 274 a.u. to the total dipole polarizability of 345 a.u. for (H(2)O)(7)(-). From the calculated polarizabilities the Rayleigh elastic light scattering properties are inferred and found to considerably enhance activity and light depolarization. (C) 2009 Elsevier B.V. All rights reserved.

Polarization and Spectral Shift of Benzophenone in Supercritical Water

Monte Carlo simulation and quantum mechanics calculations based on the INDO/CIS and TD-DFT methods were utilized to study the solvatochromic shift of benzophenone when changing the environment from normal water to supercritical (P = 340.2 atm and T = 673 K) condition. Solute polarization increases the dipole moment of benzophenone, compared to gas phase, by 88 and 35% in normal and supercritical conditions, giving the in-solvent dipole value of 5.8 and 4.2 D, respectively. The average number of solute-solvent hydrogen bonds was analyzed, and a large decrease of 2.3 in normal water to only 0.8 in the supercritical environment was found. By using these polarized models of benzophenone in the two different conditions of water, we performed MC simulations to generate statistically uncorrelated configurations of the solute surrounded by the solvent molecules and subsequent quantum mechanics calculations on these configurations. When changing from normal to supercritical water environment, INDO/CIS calculations explicitly considering all valence electrons of the 235 solvent water molecules resulted in a solvatochromic shift of 1425 cm(-1) for the most intense transition of benzophenone, that is, slightly underestimated in comparison with the experimentally inferred result of 1700 cm(-1). TD-B3LYP/6-311+G(2d,p) calculations on the same configurations but with benzophenone electrostatically embedded in the 320 water molecules resulted in a solvatochromic shift of 1715 cm(-1) for this transition, in very good agreement with the experimental result. When using the unpolarized model of the benzophenone, this calculated solvatochromic shift was only 640 cm(-1). Additional calculations were also made by using BHandHLYP/6-311+G(2d,p) to analyze the effect of the asymptotic decay of the exchange functional. This study indicates that, contrary to the general expectation, there is a sizable solute polarization even in the low-density regime of supercritical condition and that the inclusion of this polarization is important for a reliable description of the spectral shifts considered here.

Self-assembling of phenothiazine compounds investigated by small-angle X-ray scattering and electron paramagnetic resonance spectroscopy

Small-angle X-ray scattering (SAXS) and electron paramagnetic resonance (EPR) have been carried out to investigate the structure of the self-aggregates of two phenothiazine drugs, chlorpromazine (CPZ) and trifluoperazine (TFP), in aqueous solution. In the SAXS studies, drug solutions of 20 and 60 mM, at pH 4.0 and 7.0, were investigated and the best data fittings were achieved assuming several different particle form factors with a homogeneous electron density distribution in respect to the water environment. Because of the limitation of scattering intensity in the q range above 0.15 angstrom(-1), precise determination of the aggregate shape was not possible and all of the tested models for ellipsoids, cylinders, or parallelepipeds fitted the experimental data equally well. The SAXS data allows inferring, however, that CPZ molecules might self-assemble in a basis set of an orthorhombic cell, remaining as nanocrystallites in solution. Such nanocrystals are composed of a small number of unit cells (up to 10, in c-direction), with CPZ aggregation numbers of 60-80. EPR spectra of 5- and 16-doxyl stearic acids bound to the aggregates were analyzed through simulation, and the dynamic and magnetic parameters were obtained. The phenothiazine concentration in EPR experiments was in the range of 5-60 mM. Critical aggregation concentration of TFP is lower than that for CPZ, consistent with a higher hydrophobicity of TFP. At acidic pH 4.0 a significant residual motion of the nitroxide relative to the aggregate is observed, and the EPR spectra and corresponding parameters are similar to those reported for aqueous surfactant micelles. However, at pH 6.5 a significant motional restriction is observed, and the nitroxide rotational correlation times correlate very well with those estimated for the whole aggregated particle from SAXS data. This implies that the aggregate is densely packed at this pH and that the nitroxide is tightly bound to it producing a strongly immobilized EPR spectrum. Besides that, at pH 6.5 the differences in motional restriction observed between 5- and 16-DSA are small, which is different from that observed for aqueous surfactant micelles.

Performance of GEANT4 in dosimetry applications: Calculation of X-ray spectra and kerma-to-dose equivalent conversion coefficients

In order to validate the Geant4 toolkit for dosimetry applications, simulations were performed to calculate conversion coefficients h(10, alpha) from air kerma free-in-air to personal dose equivalent Hp(10, a). The simulations consisted of two parts: the production of X-rays with radiation qualities of narrow and wide spectra, and the interaction of radiation with ICRU tissue-equivalent and ISO water slab phantoms. The half-value layers of the X-ray spectra obtained by simulation were compared with experimental results. Mean energy, spectral resolution, half-value layers and conversion coefficients were compared with ISO reference values. The good agreement between results from simulation and reference data shows that the Geant4 is suitable for dosimetry applications which involve photons with energies in the range of ten to a few hundreds of keV. (C) 2008 Elsevier Ltd. All rights reserved.