SUSYFLAVOR v2: A Computational tool for FCNC and CP-violating processes in the MSSM

We present SUSY_FLAVOR version 2 — a Fortran 77 program that calculates low-energy flavor observables in the general R-parity conserving MSSM. For a set of MSSM parameters as input, the code gives predictions for: 1. Electric dipole moments of the leptons and the neutron. 2. Anomalous magnetic moments (i.e. g − 2) of the leptons. 3. Radiative lepton decays (μ → eγ and τ → μγ , eγ ). 4. Rare Kaon decays (K0 L → π0 ¯νν and K+ → π+ ¯νν). 5. Leptonic B decays (Bs,d → l+l−, B → τ ν and B → Dτ ν). 6. Radiative B decays (B → ¯ Xsγ ). 7. ΔF = 2 processes ( ¯ K0–K0, ¯D–D, ¯Bd–Bd and ¯Bs–Bs mixing). Comparing to SUSY_FLAVOR v1, where the matching conditions were calculated strictly at one-loop level, SUSY_FLAVOR v2 performs the resummation of all chirally enhanced corrections, i.e. takes into account the enhanced effects from tan β and/or large trilinear soft mixing terms to all orders in perturbation theory. Also, in SUSY_FLAVOR v2 new routines calculation of B → (D)τ ν, g − 2, radiative lepton decays and Br(l → l′γ ) were added. All calculations are done using exact diagonalization of the sfermion mass matrices. The program can be obtained from http://www.fuw.edu.pl/susy_flavor.

A stochastic principle behind polar properties of condensed molecular matter

A statistical mechanics view leads to the conclusion that polar molecules allowed to populate a degree of freedom for orientational disorder in a condensed phase thermalize into a bi-polar state featuring zero net polarity. In cases of orientational disorder polar order of condensed molecular matter can only exist in corresponding sectors of opposite average polarities. Channel type inclusion compounds, single component molecular crystals, solid solutions, optically anomalous crystals, inorganic ionic crystals, biomimetic crystals and biological tissues investigated by scanning pyroelectric and phase sensitive second harmonic generation microscopy all showed domains of opposite polarities in their final grown state. For reported polar molecular crystal structures it is assumed that kinetic hindrance along one direction of the polar axis is preventing the formation of a bi-polar state, thus allowing for a kinetically controlled mono-domain state. In this review we summarize theoretical and experimental findings leading to far reaching conclusions on the polar state of solid molecular matter. “… no stationary state … of a system has an electrical dipole moment.” P. W. Anderson, Science, 1972, 177, 393.

Pyroelectric and piezoelectric scanning microscopy applied to reveal the bipolar state of 4-iodo-4′-nitrobiphenyl (INBP)

Two recent scanning probe techniques were applied to investigate the bipolar twin state of 4-iodo-4′-nitrobiphenyl (INBP) crystals. Solution grown crystals of INBP show typically a morphology which does not express that of a mono-domain polar structure (Fdd2, mm2). From previous X-ray diffraction a twinning volume ratio of [similar]70 : 30 is now explained by two unipolar domains (Flack parameter: 0.075(29)) of opposite orientation of the molecular dipoles, joined by a transition zone showing a width of [similar]140 μm. Scanning pyroelectric microscopy (SPEM) demonstrates a continuous transition of the polarization P from +P into −P across the zone. Application of piezoelectric force microscopy (PFM) confirms unipolar alignment of INBP molecules down to a resolution of [similar]20 nm. A previously proposed real structure for INBP crystals built from lamellae with antiparallel alignment is thus rejected. Anomalous X-ray scattering was used to determine the absolute molecular orientation in the two domains. End faces of the polar axis 2 are thus made up by NO2 groups. Using a previously determined negative pyroelectric coefficient pc leads to a confirmation also by a SPEM analysis. Calculated values for functional group interactions (DA), (AA), (DD) and the stochastic theory of polarity formation allow us to predict that NO2 groups should terminate corresponding faces. Following the present analysis, INBP may represent a first example undergoing dipole reversal upon growth to end up in a bipolar state.

Changing correlation structures of the Northern Hemisphere atmospheric circulation from 1000 to 2100 AD

Atmospheric circulation modes are important concepts in understanding the variability of atmospheric dynamics. Assuming their spatial patterns to be fixed, such modes are often described by simple indices from rather short observational data sets. The increasing length of reanalysis products allows these concepts and assumptions to be scrutinised. Here we investigate the stability of spatial patterns of Northern Hemisphere teleconnections by using the Twentieth Century Reanalysis as well as several control and transient millennium-scale simulations with coupled models. The observed and simulated centre of action of the two major teleconnection patterns, the North Atlantic Oscillation (NAO) and to some extent the Pacific North American (PNA), are not stable in time. The currently observed dipole pattern of the NAO, its centre of action over Iceland and the Azores, split into a north–south dipole pattern in the western Atlantic with a wave train pattern in the eastern part, connecting the British Isles with West Greenland and the eastern Mediterranean during the period 1940–1969 AD. The PNA centres of action over Canada are shifted southwards and over Florida into the Gulf of Mexico during the period 1915–1944 AD. The analysis further shows that shifts in the centres of action of either teleconnection pattern are not related to changes in the external forcing applied in transient simulations of the last millennium. Such shifts in their centres of action are accompanied by changes in the relation of local precipitation and temperature with the overlying atmospheric mode. These findings further undermine the assumption of stationarity between local climate/proxy variability and large-scale dynamics inherent when using proxy-based reconstructions of atmospheric modes, and call for a more robust understanding of atmospheric variability on decadal timescales.

Solute transport in crystalline rocks at Äspö — I: Geological basis and model calibration

Water-conducting faults and fractures were studied in the granite-hosted A¨ spo¨ Hard Rock Laboratory (SE Sweden). On a scale of decametres and larger, steeply dipping faults dominate and contain a variety of different fault rocks (mylonites, cataclasites, fault gouges). On a smaller scale, somewhat less regular fracture patterns were found. Conceptual models of the fault and fracture geometries and of the properties of rock types adjacent to fractures were derived and used as input for the modelling of in situ dipole tracer tests that were conducted in the framework of the Tracer Retention Understanding Experiment (TRUE-1) on a scale of metres. After the identification of all relevant transport and retardation processes, blind predictions of the breakthroughs of conservative to moderately sorbing tracers were calculated and then compared with the experimental data. This paper provides the geological basis and model calibration, while the predictive and inverse modelling work is the topic of the companion paper [J. Contam. Hydrol. 61 (2003) 175]. The TRUE-1 experimental volume is highly fractured and contains the same types of fault rocks and alterations as on the decametric scale. The experimental flow field was modelled on the basis of a 2D-streamtube formalism with an underlying homogeneous and isotropic transmissivity field. Tracer transport was modelled using the dual porosity medium approach, which is linked to the flow model by the flow porosity. Given the substantial pumping rates in the extraction borehole, the transport domain has a maximum width of a few centimetres only. It is concluded that both the uncertainty with regard to the length of individual fractures and the detailed geometry of the network along the flowpath between injection and extraction boreholes are not critical because flow is largely one-dimensional, whether through a single fracture or a network. Process identification and model calibration were based on a single uranine breakthrough (test PDT3), which clearly showed that matrix diffusion had to be included in the model even over the short experimental time scales, evidenced by a characteristic shape of the trailing edge of the breakthrough curve. Using the geological information and therefore considering limited matrix diffusion into a thin fault gouge horizon resulted in a good fit to the experiment. On the other hand, fresh granite was found not to interact noticeably with the tracers over the time scales of the experiments. While fracture-filling gouge materials are very efficient in retarding tracers over short periods of time (hours–days), their volume is very small and, with time progressing, retardation will be dominated by altered wall rock and, finally, by fresh granite. In such rocks, both porosity (and therefore the effective diffusion coefficient) and sorption Kds are more than one order of magnitude smaller compared to fault gouge, thus indicating that long-term retardation is expected to occur but to be less pronounced.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

A reconstruction of radiocarbon production and total solar irradiance from the Holocene ¹⁴C and CO₂ records: implications of data and model uncertainties

Radiocarbon production, solar activity, total solar irradiance (TSI) and solar-induced climate change are reconstructed for the Holocene (10 to 0 kyr BP), and TSI is predicted for the next centuries. The IntCal09/SHCal04 radiocarbon and ice core CO2 records, reconstructions of the geomagnetic dipole, and instrumental data of solar activity are applied in the Bern3D-LPJ, a fully featured Earth system model of intermediate complexity including a 3-D dynamic ocean, ocean sediments, and a dynamic vegetation model, and in formulations linking radiocarbon production, the solar modulation potential, and TSI. Uncertainties are assessed using Monte Carlo simulations and bounding scenarios. Transient climate simulations span the past 21 thousand years, thereby considering the time lags and uncertainties associated with the last glacial termination. Our carbon-cycle-based modern estimate of radiocarbon production of 1.7 atoms cm−2 s−1 is lower than previously reported for the cosmogenic nuclide production model by Masarik and Beer (2009) and is more in-line with Kovaltsov et al. (2012). In contrast to earlier studies, periods of high solar activity were quite common not only in recent millennia, but throughout the Holocene. Notable deviations compared to earlier reconstructions are also found on decadal to centennial timescales. We show that earlier Holocene reconstructions, not accounting for the interhemispheric gradients in radiocarbon, are biased low. Solar activity is during 28% of the time higher than the modern average (650 MeV), but the absolute values remain weakly constrained due to uncertainties in the normalisation of the solar modulation to instrumental data. A recently published solar activity–TSI relationship yields small changes in Holocene TSI of the order of 1 W m−2 with a Maunder Minimum irradiance reduction of 0.85 ± 0.16 W m−2. Related solar-induced variations in global mean surface air temperature are simulated to be within 0.1 K. Autoregressive modelling suggests a declining trend of solar activity in the 21st century towards average Holocene conditions.

Constraints on the effective fluid theory of stationary branes

We develop further the effective fluid theory of stationary branes. This formalism applies to stationary blackfolds as well as to other equilibrium brane systems at finite temperature. The effective theory is described by a Lagrangian containing the information about the elastic dynamics of the brane embedding as well as the hydrodynamics of the effective fluid living on the brane. The Lagrangian is corrected order-by-order in a derivative expansion, where we take into account the dipole moment of the brane which encompasses finite-thickness corrections, including transverse spin. We describe how to extract the thermodynamics from the Lagrangian and we obtain constraints on the higher-derivative terms with one and two derivatives. These constraints follow by comparing the brane thermodynamics with the conserved currents associated with background Killing vector fields. In particular, we fix uniquely the one- and two-derivative terms describing the coupling of the transverse spin to the background space-time. Finally, we apply our formalism to two blackfold examples, the black tori and charged black rings and compare the latter to a numerically generated solution.

Lepton Flavor Violation in the Standard Model with general Dimension-Six Operators

We study lepton flavor observables in the Standard Model (SM) extended with all dimension-6 operators which are invariant under the SM gauge group. We calculate the complete one-loop predictions to the radiative lepton decays μ → eγ, τ → μγ and τ → eγ as well as to the closely related anomalous magnetic moments and electric dipole moments of charged leptons, taking into account all dimension-6 operators which can generate lepton flavor violation. Also the 3-body flavor violating charged lepton decays τ ± → μ ± μ + μ −, τ ± → e ± e + e −, τ ± → e ± μ + μ −, τ ± → μ ± e + e −, τ ± → e ∓ μ ± μ ±, τ ± → μ ∓ e ± e ± and μ ± → e ± e + e − and the Z 0 decays Z 0 → ℓ+iℓ−j are considered, taking into account all tree-level contributions.

Contrast Formation in Kelvin Probe Force Microscopy of Single π-Conjugated Molecules

We report the contrast formation in the local contact potential difference (LCPD) measured by Kelvin probe force microscopy (KPFM) on single charge-transfer complexes (CTCs) on a NaCl bilayer on Cu(111). At different tip heights, we found quantitatively different LCPD contrasts that characterize different properties of the molecule. In the small distance regime, the tip penetrates the electron density of the molecule, and the contrast is related to the size and topography of the electron shell of the molecule. For larger distances, the LCPD contrast corresponds to the electrostatic field above the molecule. However, in the medium-distance regime, that is, for tip heights similar to the size of the molecule, the nonspherical distribution of π- and σ-electrons often conceals the effect of the partial charges within the molecule. Only for large distances does the LCPD map converge toward the simple field of a dipole for a polar molecule.

Excitonic Splitting, Delocalization, and Vibronic Quenching in the Benzonitrile Dimer

The excitonic S1/S2 state splitting and the localization/delocalization of the S1 and S2 electronic states are investigated in the benzonitrile dimer (BN)2 and its 13C and d5 isotopomers by mass-resolved two-color resonant two-photon ionization spectroscopy in a supersonic jet, complemented by calculations. The doubly hydrogen-bonded (BN-h5)2 and (BN-d5)2 dimers are C2h symmetric with equivalent BN moieties. Only the S0 → S2 electronic origin is observed, while the S0 → S1 excitonic component is electric-dipole forbidden. A single 12C/13C or 5-fold h5/d5 isotopic substitution reduce the dimer symmetry to Cs, so that the heteroisotopic dimers (BN)2-(h5 – h513C), (BN)2-(h5 – d5), and (BN)2-(h5 – h513C) exhibit both S0 → S1 and S0 → S2 origins. Isotope-dependent contributions Δiso to the excitonic splittings arise from the changes of the BN monomer zero-point vibrational energies; these range from Δiso(12C/13C) = 3.3 cm–1 to Δiso(h5/d5) = 155.6 cm–1. The analysis of the experimental S1/S2 splittings of six different isotopomeric dimers yields the S1/S2 exciton splitting Δexc = 2.1 ± 0.1 cm–1. Since Δiso(h5/d5) ≫ Δexc and Δiso(12C/13C) > Δexc, complete and near-complete exciton localization occurs upon 12C/13C and h5/d5 substitutions, respectively, as diagnosed by the relative S0 → S1 and S0 → S2 origin band intensities. The S1/S2 electronic energy gap of (BN)2 calculated by the spin-component scaled approximate second-order coupled-cluster (SCS-CC2) method is Δelcalc = 10 cm–1. This electronic splitting is reduced by the vibronic quenching factor Γ. The vibronically quenched exciton splitting Δelcalc·Γ = Δvibroncalc = 2.13 cm–1 is in excellent agreement with the observed splitting Δexc = 2.1 cm–1. The excitonic splittings can be converted to semiclassical exciton hopping times; the shortest hopping time is 8 ps for the homodimer (BN-h5)2, the longest is 600 ps for the (BN)2(h5 – d5) heterodimer.

Excited-State Structure, Vibrations, and Nonradiative Relaxation of Jet-Cooled 5-Fluorocytosine

The S0 → S1 vibronic spectrum and S1 state nonradiative relaxation of jet-cooled keto-amino 5-fluorocytosine (5FCyt) are investigated by two-color resonant two-photon ionization spectroscopy at 0.3 and 0.05 cm–1 resolution. The 000 rotational band contour is polarized in-plane, implying that the electronic transition is 1ππ*. The electronic transition dipole moment orientation and the changes of rotational constants agree closely with the SCS-CC2 calculated values for the 1ππ* (S1) transition of 5FCyt. The spectral region from 0 to 300 cm–1 is dominated by overtone and combination bands of the out-of-plane ν1′ (boat), ν2′ (butterfly), and ν3′ (HN–C6H twist) vibrations, implying that the pyrimidinone frame is distorted out-of-plane by the 1ππ* excitation, in agreement with SCS-CC2 calculations. The number of vibronic bands rises strongly around +350 cm–1; this is attributed to the 1ππ* state barrier to planarity that corresponds to the central maximum of the double-minimum out-of-plane vibrational potentials along the ν1′, ν2′, and ν3′ coordinates, which gives rise to a high density of vibronic excitations. At +1200 cm–1, rapid nonradiative relaxation (knr ≥ 1012 s–1) sets in, which we interpret as the height of the 1ππ* state barrier in front of the lowest S1/S0 conical intersection. This barrier in 5FCyt is 3 times higher than that in cytosine. The lifetimes of the ν′ = 0, 2ν1′, 2ν2′, 2ν1′ + 2ν2′, 4ν2′, and 2ν1′ + 4ν2′ levels are determined from Lorentzian widths fitted to the rotational band contours and are τ ≥ 75 ps for ν′ = 0, decreasing to τ ≥ 55 ps at the 2ν1′ + 4ν2′ level at +234 cm–1. These gas-phase lifetimes are twice those of S1 state cytosine and 10–100 times those of the other canonical nucleobases in the gas phase. On the other hand, the 5FCyt gas-phase lifetime is close to the 73 ps lifetime in room-temperature solvents. This lack of dependence on temperature and on the surrounding medium implies that the 5FCyt nonradiative relaxation from its S1 (1ππ*) state is essentially controlled by the same ∼1200 cm–1 barrier and conical intersection both in the gas phase and in solution.

On vertical electric fields at lunar magnetic anomalies

We study the interaction between a magnetic dipole mimicking the Gerasimovich magnetic anomaly on the lunar surface and the solar wind in a self-consistent 3-D quasi-neutral hybrid simulation where ions are modeled as particles and electrons as a charge-neutralizing fluid. Especially, we consider the origin of the recently observed electric potentials at lunar magnetic anomalies. An antimoonward Hall electric field forms in our simulation resulting in a potential difference of <300V on the lunar surface, in which the value is similar to observations. Since the hybrid model assumes charge neutrality, our results suggest that the electric potentials at lunar magnetic anomalies can be formed by decoupling of ion and electron motion even without charge separation.

Luminescence and Energy Transfer of [Ru(bpy) 3 ] 2+ , [Cr(ox) 3 ] 3- , and [Os(bpy) 3 ] 2+ in Three-Dimensional Oxalato-Networks

Luminescence and energy transfer in [Zn1-xRux(bpy)3][NaAl1-yCry(ox)3] (x ≈ 0.01, y = 0.006 − 0.22; bpy = 2,2‘-bipyridine, ox = C2O42-) and [Zn1-x-yRuxOsy(bpy)3][NaAl(ox)3] (x ≈ 0.01, y = 0.012) are presented and discussed. Surprisingly, the luminescence of the isolated luminophores [Ru(bpy)3]2+ and [Os(bpy)3]2+ in [Zn(bpy)3][NaAl(ox)3] is hardly quenched at room temperature. Steady-state luminescence spectra and decay curves show that energy transfer occurs between [Ru(bpy)3]2+ and [Cr(ox)3]3- and between [Ru(bpy)3]2+ and [Os(bpy)3]2+ in [Zn1-xRux(bpy)3][NaAl1-yCry(ox)3] and [Zn1-x-yRuxOsy(bpy)3] [NaAl(ox)3], respectively. For a quantitative investigation of the energy transfer, a shell type model is developed, using a Monte Carlo procedure and the structural parameters of the systems. A good description of the experimental data is obtained assuming electric dipole−electric dipole interaction between donors and acceptors, with a critical distance Rc for [Ru(bpy)3]2+ to [Cr(ox)3]3- energy transfer of 15 Å and for [Ru(bpy)3]2+ to [Os(bpy)3]2+ energy transfer of 33 Å. These values are in good agreement with those derived using the Förster−Dexter theory.

Zeeman effect in atmospheric O₂ measured by ground-based microwave radiometry

In this work we study the Zeeman effect on stratospheric O₂ using ground-based microwave radiometer measurements. The interaction of the Earth magnetic field with the oxygen dipole leads to a splitting of O₂ energy states, which polarizes the emission spectra. A special campaign was carried out in order to measure this effect in the oxygen emission line centered at 53.07 GHz. Both a fixed and a rotating mirror were incorporated into the TEMPERA (TEMPERature RAdiometer) in order to be able to measure under different observational angles. This new configuration allowed us to change the angle between the observational path and the Earth magnetic field direction. Moreover, a high-resolution spectrometer (1 kHz) was used in order to measure for the first time the polarization state of the radiation due to the Zeeman effect in the main isotopologue of oxygen from ground-based microwave measurements. The measured spectra showed a clear polarized signature when the observational angles were changed, evidencing the Zeeman effect in the oxygen molecule. In addition, simulations carried out with the Atmospheric Radiative Transfer Simulator (ARTS) allowed us to verify the microwave measurements showing a very good agreement between model and measurements. The results suggest some interesting new aspects for research of the upper atmosphere.