Quantum Monte Carlo study of small aluminum clusters Al-n (n=2-13)

Using fixed node diffusion quantum Monte Carlo (FN-DMC) simulations and density functional theory (DFT) within the generalized gradient approximations, we calculate the total energies of the relaxed and unrelaxed neutral, cationic, and anionic aluminum clusters, Al-n (n = 1-13). From the obtained total energies, we extract the ionization potential and electron detachment energy and compare with previous theoretical and experimental results. Our results for the electronic properties from both the FN-DMC and DFT calculations are in reasonably good agreement with the available experimental data. A comparison between the FN-DMC and DFT results reveals that their differences are a few tenths of electron volt for both the ionization potential and the electron detachment energy. We also observe two distinct behaviors in the electron correlation contribution to the total energies from smaller to larger clusters, which could be assigned to the structural transition of the clusters from planar to three-dimensional occurring at n = 4 to 5.

Quantum critical scaling at a Bose-glass/superfluid transition: Theory and experiment for a model quantum magnet

In this paper we investigate the quantum phase transition from magnetic Bose Glass to magnetic Bose-Einstein condensation induced by amagnetic field in NiCl2 center dot 4SC(NH2)(2) (dichloro-tetrakis-thiourea-nickel, or DTN), doped with Br (Br-DTN) or site diluted. Quantum Monte Carlo simulations for the quantum phase transition of the model Hamiltonian for Br-DTN, as well as for site-diluted DTN, are consistent with conventional scaling at the quantum critical point and with a critical exponent z verifying the prediction z = d; moreover the correlation length exponent is found to be nu = 0.75(10), and the order parameter exponent to be beta = 0.95(10). We investigate the low-temperature thermodynamics at the quantum critical field of Br-DTN both numerically and experimentally, and extract the power-law behavior of the magnetization and of the specific heat. Our results for the exponents of the power laws, as well as previous results for the scaling of the critical temperature to magnetic ordering with the applied field, are incompatible with the conventional crossover-scaling Ansatz proposed by Fisher et al. [Phys. Rev. B 40, 546 (1989)]. However they can all be reconciled within a phenomenological Ansatz in the presence of a dangerously irrelevant operator.

Characterization of an extrapolation chamber for low-energy X-rays: Experimental and Monte Carlo preliminary results

The extrapolation chamber is a parallel-plate ionization chamber that allows variation of its air-cavity volume. In this work, an experimental study and MCNP-4C Monte Carlo code simulations of an ionization chamber designed and constructed at the Calibration Laboratory at IFEN to be used as a secondary dosimetry standard for low-energy X-rays are reported. The results obtained were within the international recommendations, and the simulations showed that the components of the extrapolation chamber may influence its response up to 11.0%. (C) 2011 Elsevier Ltd. All rights reserved.

Optimization of x-ray spectra in digital mammography through Monte Carlo simulations

In this work, a Monte Carlo code was used to investigate the performance of different x-ray spectra in digital mammography, through a figure of merit (FOM), defined as FOM = CNR2/(D) over bar (g), with CNR being the contrast-to-noise ratio in image and (D) over bar (g) being the average glandular dose. The FOM was studied for breasts with different thicknesses t (2 cm <= t <= 8 cm) and glandular contents (25%, 50% and 75% glandularity). The anode/filter combinations evaluated were those traditionally employed in mammography (Mo/Mo, Mo/Rh, Rh/Rh), and a W anode combined with Al or K-edge filters (Zr, Mo, Rh, Pd, Ag, Cd, Sn), for tube potentials between 22 and 34 kVp. Results show that the W anode combined with K-edge filters provides higher values of FOM for all breast thicknesses investigated. Nevertheless, the most suitable filter and tube potential depend on the breast thickness, and for t >= 6 cm, they also depend on breast glandularity. Particularly for thick and dense breasts, a W anode combined with K-edge filters can greatly improve the digital technique, with the values of FOM up to 200% greater than that obtained with the anode/filter combinations and tube potentials traditionally employed in mammography. For breasts with t < 4 cm, a general good performance was obtained with the W anode combined with 60 mu m of the Mo filter at 24-25 kVp, while 60 mu m of the Pd filter provided a general good performance at 24-26 kVp for t = 4 cm, and at 28-30 and 29-31 kVp for t = 6 and 8 cm, respectively.

Low-temperature thermodynamic properties near the field-induced quantum critical point in NiCl2-4SC(NH2)(2)

We present a comprehensive experimental and theoretical investigation of the thermodynamic properties: specific heat, magnetization, and thermal expansion in the vicinity of the field-induced quantum critical point (QCP) around the lower critical field H-c1 approximate to 2 T in NiCl2-4SC(NH2)(2). A T-3/2 behavior in the specific heat and magnetization is observed at very low temperatures at H = H-c1, which is consistent with the universality class of Bose-Einstein condensation of magnons. The temperature dependence of the thermal expansion coefficient at H-c1 shows minor deviations from the expected T-1/2 behavior. Our experimental study is complemented by analytical calculations and quantum Monte Carlo simulations, which reproduce nicely the measured quantities. We analyze the thermal and the magnetic Gruneisen parameters, which are ideal quantities to identify QCPs. Both parameters diverge at H-c1 with the expected T-1 power law. By using the Ehrenfest relations at the second-order phase transition, we are able to estimate the pressure dependencies of the characteristic temperature and field scales.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

Monte Carlo Semi-Empirical Model for Si(Li) X-Ray Detector: Differences between Nominal and Fitted Parameters.

A detailed characterization of a X-ray Si(Li) detector was performed to obtain the energy dependence of efficiency in the photon energy range of 6.4 - 59.5 keV. which was measured and reproduced by Monte Carlo (MC) simulations. Significant discrepancies between MC and experimental values were found when lhe manufacturer parameters of lhe detector were used in lhe simulation. A complete Computerized Tomagraphy (CT) detector scan allowed to find the correct crystal dimensions and position inside the capsule. The computed efficiencies with the resulting detector model differed with the measured values no more than 10% in most of the energy range.

Hybrid density functional study of small Rh-n (n=2-15) clusters

The physical properties of small rhodium clusters, Rh-n, have been in debate due to the shortcomings of density functional theory (DFT). To help in the solution of those problems, we obtained a set of putative lowest energy structures for small Rh-n (n = 2-15) clusters employing hybrid-DFT and the generalized gradient approximation (GGA). For n = 2-6, both hybrid and GGA functionals yield similar ground-state structures (compact), however, hybrid favors compact structures for n = 7-15, while GGA favors open structures based on simple cubic motifs. Thus, experimental results are crucial to indicate the correct ground-state structures, however, we found that a unique set of structures (compact or open) is unable to explain all available experimental data. For example, the GGA structures (open) yield total magnetic moments in excellent agreement with experimental data, while hybrid structures (compact) have larger magnetic moments compared with experiments due to the increased localization of the 4d states. Thus, we would conclude that GGA provides a better description of the Rh-n clusters, however, a recent experimental-theoretical study [ Harding et al., J. Chem. Phys. 133, 214304 (2010)] found that only compact structures are able to explain experimental vibrational data, while open structures cannot. Therefore, it indicates that the study of Rh-n clusters is a challenging problem and further experimental studies are required to help in the solution of this conundrum, as well as a better description of the exchange and correlation effects on the Rh n clusters using theoretical methods such as the quantum Monte Carlo method.

Structure and electronic properties of hydrated mesityl oxide: a sequential quantum mechanics/molecular mechanics approach

The hydration of mesityl oxide (MOx) was investigated through a sequential quantum mechanics/molecular mechanics approach. Emphasis was placed on the analysis of the role played by water in the MOx syn-anti equilibrium and the electronic absorption spectrum. Results for the structure of the MOx-water solution, free energy of solvation and polarization effects are also reported. Our main conclusion was that in gas-phase and in low-polarity solvents, the MOx exists dominantly in syn-form and in aqueous solution in anti-form. This conclusion was supported by Gibbs free energy calculations in gas phase and in-water by quantum mechanical calculations with polarizable continuum model and thermodynamic perturbation theory in Monte Carlo simulations using a polarized MOx model. The consideration of the in-water polarization of the MOx is very important to correctly describe the solute-solvent electrostatic interaction. Our best estimate for the shift of the pi-pi* transition energy of MOx, when it changes from gas-phase to water solvent, shows a red-shift of -2,520 +/- 90 cm(-1), which is only 110 cm(-1) (0.014 eV) below the experimental extrapolation of -2,410 +/- 90 cm(-1). This red-shift of around -2,500 cm(-1) can be divided in two distinct and opposite contributions. One contribution is related to the syn -> anti conformational change leading to a blue-shift of similar to 1,700 cm(-1). Other contribution is the solvent effect on the electronic structure of the MOx leading to a red-shift of around -4,200 cm(-1). Additionally, this red-shift caused by the solvent effect on the electronic structure can by composed by approximately 60 % due to the electrostatic bulk effect, 10 % due to the explicit inclusion of the hydrogen-bonded water molecules and 30 % due to the explicit inclusion of the nearest water molecules.

Nuclear-modification factor for open-heavy-flavor production at forward rapidity in Cu plus Cu collisions at root s(NN)=200 GeV

Background: Heavy-flavor production in p + p collisions is a good test of perturbative-quantum-chromodynamics (pQCD) calculations. Modification of heavy-flavor production in heavy-ion collisions relative to binary-collision scaling from p + p results, quantified with the nuclear-modification factor (R-AA), provides information on both cold-and hot-nuclear-matter effects. Midrapidity heavy-flavor R-AA measurements at the Relativistic Heavy Ion Collider have challenged parton-energy-loss models and resulted in upper limits on the viscosity-entropy ratio that are near the quantum lower bound. Such measurements have not been made in the forward-rapidity region. Purpose: Determine transverse-momentum (p(T)) spectra and the corresponding R-AA for muons from heavy-flavor meson decay in p + p and Cu + Cu collisions at root s(NN) = 200 GeV and y = 1.65. Method: Results are obtained using the semileptonic decay of heavy-flavor mesons into negative muons. The PHENIX muon-arm spectrometers measure the p(T) spectra of inclusive muon candidates. Backgrounds, primarily due to light hadrons, are determined with a Monte Carlo calculation using a set of input hadron distributions tuned to match measured-hadron distributions in the same detector and statistically subtracted. Results: The charm-production cross section in p + p collisions at root s = 200 GeV, integrated over p(T) and in the rapidity range 1.4 < y < 1.9, is found to be d(sigma e (e) over bar)/dy = 0.139 +/- 0.029 (stat)(-0.058)(+0.051) (syst) mb. This result is consistent with a perturbative fixed-order-plus-next-to-leading-log calculation within scale uncertainties and is also consistent with expectations based on the corresponding midrapidity charm-production cross section measured by PHENIX. The R-AA for heavy-flavor muons in Cu + Cu collisions is measured in three centrality bins for 1 < p(T) < 4 GeV/c. Suppression relative to binary-collision scaling (R-AA < 1) increases with centrality. Conclusions: Within experimental and theoretical uncertainties, the measured charm yield in p + p collisions is consistent with state-of-the-art pQCD calculations. Suppression in central Cu + Cu collisions suggests the presence of significant cold-nuclear-matter effects and final-state energy loss.

Theoretical study of the absorption and nonradiative deactivation of 1-nitronaphthalene in the low-lying singlet and triplet excited states including methanol and ethanol solvent effects

The photophysics of the 1-nitronaphthalene molecular system, after the absorption transition to the first singlet excited state, is theoretically studied for investigating the ultrafast multiplicity change to the triplet manifold. The consecutive transient absorption spectra experimentally observed in this molecular system are also studied. To identify the electronic states involved in the nonradiative decay, the minimum energy path of the first singlet excited state is obtained using the complete active space self-consistent field//configurational second-order perturbation approach. A near degeneracy region was found between the first singlet and the second triplet excited states with large spin-orbit coupling between them. The intersystem crossing rate was also evaluated. To support the proposed deactivation model the transient absorption spectra observed in the experiments were also considered. For this, computer simulations using sequential quantum mechanic-molecular mechanic methodology was used to consider the solvent effect in the ground and excited states for proper comparison with the experimental results. The absorption transitions from the second triplet excited state in the relaxed geometry permit to describe the transient absorption band experimentally observed around 200 fs after the absorption transition. This indicates that the T-2 electronic state is populated through the intersystem crossing presented here. The two transient absorption bands experimentally observed between 2 and 45 ps after the absorption transition are described here as the T-1 -> T-3 and T-1 -> T-5 transitions, supporting that the intermediate triplet state (T-2) decays by internal conversion to T-1. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4738757]

Electronic Properties of Water in Liquid Environment. A Sequential QM/MM Study Using the Free Energy Gradient Method

There is a continuous search for theoretical methods that are able to describe the effects of the liquid environment on molecular systems. Different methods emphasize different aspects, and the treatment of both the local and bulk properties is still a great challenge. In this work, the electronic properties of a water molecule in liquid environment is studied by performing a relaxation of the geometry and electronic distribution using the free energy gradient method. This is made using a series of steps in each of which we run a purely molecular mechanical (MM) Monte Carlo Metropolis simulation of liquid water and subsequently perform a quantum mechanical/molecular mechanical (QM/MM) calculation of the ensemble averages of the charge distribution, atomic forces, and second derivatives. The MP2/aug-cc-pV5Z level is used to describe the electronic properties of the QM water. B3LYP with specially designed basis functions are used for the magnetic properties. Very good agreement is found for the local properties of water, such as geometry, vibrational frequencies, dipole moment, dipole polarizability, chemical shift, and spin-spin coupling constants. The very good performance of the free energy method combined with a QM/MM approach along with the possible limitations are briefly discussed.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

Ionization of chlorophyll-c(2) in liquid methanol

The ionization of chlorophyll-c(2) in liquid methanol was investigated by a sequential quantum mechanical/Monte Carlo approach. Focus was placed on the determination of the first ionization energy of chlorophyll-c(2). The results show that the first vertical ionization energy (IE) is red-shifted by 0.47 +/- 0.24 eV relative to the gas-phase value. The red-shift of the chlorophyll-c(2) IE in the liquid phase can be explained by Mg center dot center dot center dot OH hydrogen bonding and long-ranged electrostatic interactions in solution. The ionization threshold for chlorophyll-c2 in liquid methanol is close to 6 eV. (C) 2012 Elsevier B.V. All rights reserved.

Self-Assembly of Poly(4-vynil-N-alkylpyridinium) Bromide onto Silica: Effect of Side-Chain Length on Structural Aspects at a Molecular Level

Self-assembly of poly(4-vynil-N-alkyl)pyridinium bromide with alkyl side chains of 2, 5, 7, 10, or 16 carbons from ethanolic solutions onto flat silica surfaces was studied by means of ellipsometry, atomic force microscopy (AFM), contact angle measurements, and sum-frequency generation (SFG) vibrational spectroscopy in the CH3 and CH2 stretch region. Ab initio quantum-chemical calculations on the N-alkylpyridinium side-group with restricted Hartree-Fock (RHF) method and 6-311G (d,p) basis set were C one to estimate the charge distribution along the pyridinium ring and the alkyl side-chain. SFG results showed that longer side chains promote the disorientation of the alkyl groups at the surface, corroborating with the contact angle values. AFM images revealed film homogeneity, regardless the alkyl side group. However, after 24 h contact with water, ringlike structures appeared on the film surfaces, when the polycation alkyl side chain had 7 or less carbons, and as the alkyl chain increased to 10 or 16 carbons, the films dewetted because the hydrophobic interactions prevailed over the electrostatic interactions between the pyridinium charged groups and the negatively charged SiO2 surface. Under acid conditions (HCl 0.1 mol.L-1), the film mean thickness values decreased up to 50% of original values when the alkyl side chains were ethyl or pentyl groups due to ion-pair disruption, but for longer groups they remained unchanged. Quantum-chemical optimization and Mulliken electron population showed that (i) from C2 to C15 the positive charge at the headgroup (HG) decreased 0.025, while the charge at combined HG + alpha-CH2 increased 0.037; and (ii) for C6 or longer, the alkyl side group presents a tilt in the geometry, moving away from the plane. Such effects summed up over the whole polymer chain give support to suggest that when the side chains are longer than 7 carbons, the hydrophobic interaction decreases film stability and increases acid resistance.