2,256 research outputs found

    Theoretical study of the finite temperature spectroscopy in van der Waals clusters. III Solvated Chromophore as an effective diatomics

    Full text link
    The absorption spectroscopy of calcium-doped argon clusters is described in terms of an effective diatomics molecule Ca-(Ar_n), in the framework of semiclassical vertical transitions. We show how, upon choosing a suitable reaction coordinate, the effective finite-temperature equilibrium properties can be obtained for the ground- and excited-surfaces from the potential of mean force (PMF). An extension of the recent multiple range random-walk method is used to calculate the PMF over continuous intervals of distances. The absorption spectra calculated using this single-coordinate description are found to be in good agreement with the spectra obtained from high-statistics Monte Carlo data, in various situations. For CaAr13_{13}, we compare the performances of two different choices of the reaction coordinate. For CaAr_37, the method is seen to be accurate enough to distinguish between different low-energy structures. Finally, the idea of casting the initial many-body problem into a single degree of freedom problem is tested on the spectroscopy of calcium in bulk solid argon.Comment: 8 pages, 9 figure

    Stacked clusters of polycyclic aromatic hydrocarbon molecules

    Full text link
    Clusters of polycyclic aromatic hydrocarbon (PAH) molecules are modelled using explicit all-atom potentials using a rigid body approximation. The PAH's considered range from pyrene (C10H8) to circumcoronene (C54H18), and clusters containing between 2 and 32 molecules are investigated. In addition to the usual repulsion-dispersion interactions, electrostatic point-charge interactions are incorporated, as obtained from density functional theory calculations. The general electrostatic distribution in neutral or singly charged PAH's is reproduced well using a fluctuating charges analysis, which provides an adequate description of the multipolar distribution. Global optimization is performed using a variety of methods, including basin-hopping and parallel tempering Monte Carlo. We find evidence that stacking the PAH molecules generally yields the most stable motif. A structural transition between one-dimensional stacks and three-dimensional shapes built from mutiple stacks is observed at larger sizes, and the threshold for this transition increases with the size of the monomer. Larger aggregates seem to evolve toward the packing observed for benzene in bulk.Difficulties met in optimizing these clusters are analysed in terms of the strong anisotropy of the molecules. We also discuss segregation in heterogeneous clusters and vibrational properties in the context of astrophysical observations.Comment: 12 pages, 7 figure

    Channelized melt flow in downwelling mantle: Implications for 226Ra-210Pb disequilibria in arc magmas

    Get PDF
    We present the results of an analytical model of porous flow of viscous melt into a steadily dilating ‘‘channel’’ (defined as a cluster of smaller veins) in downwelling subarc mantle. The model predicts the pressure drop in the mantle wedge matrix surrounding the channel needed to drive melt flow as a function of position and time. Melt is sucked toward the dilatant region at a near-constant velocity (105 s1) until veins comprising the channel stop opening (t = t). Fluid elements that complete their journey within the time span t < t arrive at a channel. Our results make it possible to calculate the region of influence sampled by melt that surrounds the channel. This region is large compared to the model size of the channelized region driving flow. For a baseline dilation time of 1 year and channel half width of 2 m, melt can be sampled over an 80-m radius and has the opportunity to sample matrix material with potentially contrasting chemistry on geologically short timescales. Our mechanical results are consistent with a downgoing arc mantle wedge source region where melting and melt extraction by porous flow to a channel network are sufficiently rapid to preserve source-derived 238U-230Th-226Ra, and potentially also 226 Ra-210Pb, disequilibria, prior to magma ascent to the surface. Since this is the rate-determining step in the overall process, it allows the possibility that such short-lived disequilibria measured in arc rocks at the surface are derived from deep in the mantle wedge. Stresses due to partial melting do not appear capable of producing the desired sucking effect, while the order of magnitude rate of shear required to drive dilation of 107 s1 is much larger than values resulting from steady state subduction. We conclude that local deformation rates in excess of background plate tectonic rates are needed to ‘‘switch on’’ the dilatant channel network and to initiate the sucking effect

    Theoretical study of finite temperature spectroscopy in van der Waals clusters. II Time-dependent absorption spectra

    Full text link
    Using approximate partition functions and a master equation approach, we investigate the statistical relaxation toward equilibrium in selected CaArn_n clusters. The Gaussian theory of absorption (previous article) is employed to calculate the average photoabsorption intensity associated with the 4s^2-> 4s^14p^1 transition of calcium as a function of time during relaxation. In CaAr_6 and CaAr_10 simple relaxation is observed with a single time scale. CaAr_13 exhibits much slower dynamics and the relaxation occurs over two distinct time scales. CaAr_37 shows much slower relaxation with multiple transients, reminiscent of glassy behavior due to competition between different low-energy structures. We interpret these results in terms of the underlying potential energy surfaces for these clusters.Comment: 10 pages, 9 figure

    Response to the Letter to the Editor

    Full text link
    This paper has attracted interest around the world from the media (both TV and newspapers). In addition, we have received letters, emails and telephone calls. One of our favorites was a voicemail message asking us to return a call to Australia at which point we would learn who really killed JFK. We welcome the opportunity to respond to the letter to the editor from Mr. Fiorentino. Mr. Fiorentino claims that our ``statement relating to the likelihood of a second assassin based on the premise of three or more separate bullets is demonstrably false.'' In response we would like to simply quote from page 327 of Gerald Posner's book Case Closed, one of the most well known works supporting the single assassin theory: ``If Connally was hit by another bullet, it had to be fired from a second shooter, since the Warren Commission's own reconstructions showed that Oswald could not have operated the bolt and refired in 1.4 seconds.'' Mr. Fiorentino also claims that the ``second fatal flaw is the use of a rather uncomplicated formula based on Bayes Theorem.'' Let EE denote the evidence and TT denote the theory that there were just two bullets (and hence a single shooter). We used Bayes Theorem to hypothetically calculate P(T∣E)P(T|E) from P(E∣T)P(E|T) and the prior probability P(T)P(T). In order to make P(T∣E)P(T|E) ten times more likely than P(Tˉ∣E)P(\bar{T}|E), the ratio of the prior probabilities [i.e., P(T)/P(Tˉ)P(T) / P(\bar{T})] would have to be greater than 15. Thus, we again conclude that this casts serious doubt on Dr. Guinn's conclusion that the evidence supported just two bullets. Sadly, this is far from the first time that probability has been misunderstood and/or misapplied in a case of public interest. A notable British example is the Clark case. See Nobles and Schiff (2005) for details. Finally, we welcome and, in fact, encourage members of the scientific community to provide alternative analyses of the data.Comment: Published in at http://dx.doi.org/10.1214/07-AOAS154 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org
    • …
    corecore