We use our recent electric dipole moment (EDM) measurement data to constrain the possibility that the HfF$^+$ EDM oscillates in time due to interactions with candidate dark matter axion-like particles (ALPs). We employ a Bayesian analysis method which accounts for both the look-elsewhere effect and the uncertainties associated with stochastic density fluctuations in the ALP field. We find no evidence of an oscillating EDM over a range spanning from 27 nHz to 400 mHz, and we use this result to constrain the ALP-gluon coupling over the mass range $10^{-22}-10^{-15}$ eV. This is the first laboratory constraint on the ALP-gluon coupling in the $10^{-17}-10^{-15}$ eV range, and the first laboratory constraint to properly account for the stochastic nature of the ALP field.

We advocate an idea that some mysterious explosions, the so-called sky-quakes, which have been known for centuries could be a manifestation of the dark matter Axion Quark Nuggets (AQN) when they propagate in the Earth's atmosphere. We specifically study the event which occurred on July 31-st 2008 and was properly recorded by the dedicated Elginfield Infrasound Array (ELFO) near London, Ontario, Canada. The infrasound detection was accompanied by non-observation of any meteors by an all-sky camera network. Our interpretation is based on the AQN dark matter model which was originally invented with a completely different purpose---to explain the similarity of the dark and visible cosmological matter densities $\Omega_{\rm dark}\sim \Omega_{\rm visible}$. Our estimates for the infrasonic frequency $\nu\simeq 5$ Hz and overpressure $\delta p\sim 0.3$ Pa are consistent with the ELFO record. We propose a detection strategy for a systematic study to search for such explosions originating from AQNs by using Distributed Acoustic Sensing and briefly mention other possible detection methods. Specific signals from AQN tracks may also be detected by an existing network of seismic stations.

We advocate an idea that the presence of the daily and annual modulations of the axion flux on the Earth surface may dramatically change the strategy of the axion searches. Our computations are based on the so-called Axion Quark Nugget (AQN) dark matter model which was originally put forward to explain the similarity of the dark and visible cosmological matter densities $\Omega_{\rm dark}\sim \Omega_{\rm visible}$. In our framework, the population of galactic axions with mass $ 10^{-6} {\rm eV}\lesssim m_a\lesssim 10^{-3}{\rm eV}$ and velocity $< v_a>\sim 10^{-3} c$ will be always accompanied by the axions with typical velocities $<v_a>\sim 0.6 c$ emitted by AQNs. We formulate the broadband detection strategy to search for such relativistic axions by studying the daily and annual modulations. We describe several tests which could effectively discriminate a true signal from noise. These AQN-originated axions can be observed as correlated events which could be recorded by synchronized stations in the global network. The correlations can be effectively studied if the detectors are positioned at distances shorter than a few hundred kilometres.

We discuss the sensitivity of the present and near-future axion dark matter experiments to a halo of axions or axion-like particles gravitationally bound to the Earth or the Sun. The existence of such halos, assuming they are formed, renders a significant gain in the sensitivity of axion searches while satisfying all the present experimental bounds. The structure and coherence properties of these halos also imply novel signals, which can depend on the latitude or orientation of the detector. We demonstrate this by analysing the sensitivity of several distinct types of axion dark matter experiments.

We investigate the possibility for the direct detection of low mass (GeV scale) WIMP dark matter in scintillation experiments. Such WIMPs are typically too light to leave appreciable nuclear recoils, but may be detected via their scattering off atomic electrons. In particular, the DAMA Collaboration [R. Bernabei et al., Nucl. Phys. At. Energy 19, 307 (2018)] has recently presented strong evidence of an annual modulation in the scintillation rate observed at energies as low as 1 keV. Despite a strong enhancement in the calculated event rate at low energies, we find that an interpretation in terms of electron-interacting WIMPs cannot be consistent with existing constraints. We also demonstrate the importance of correct treatment of the atomic wavefunctions, and show the resulting event rate is very sensitive to the low-energy performance of the detectors, meaning it is crucial that the detector uncertainties be taken into account. Finally, we demonstrate that the potential scintillation event rate can be much larger than may otherwise be expected, meaning that competitive searches can be performed for GeV scale WIMPs using the conventional prompt S1 scintillation signals. This is important given the recent and upcoming very large liquid xenon detectors.

Interaction with the axion dark matter (DM) field generates an oscillating nuclear electric dipole moment (EDM) with a frequency corresponding to the axion's Compton frequency. Within an atom, an oscillating EDM can drive electric dipole transitions in the electronic shell. In the absence of radiation, and if the axion frequency matches a dipole transition, it can promote the electron into the excited state. The excitation events can be detected, for example, via subsequent uorescence or photoionization. Here we calculate the rates of such transitions. For a single light atom and an axion Compton frequency resonant with a transition energy corresponding to 1 eV, the rate is on the order of 10^(-22) per year, so a macroscopic atomic sample would be needed. A fundamental challenge is discriminating against background processes that may lead to the excitation of the same electric dipole transition. The ways to enhance the signals to potentially observable levels exceeding backgrounds and to search for axions in an extended frequency range are discussed.

Dark photon is a massive vector field which interacts only with the physical photon through the kinetic mixing. This coupling is assumed to be weak so that the dark photon becomes almost unobservable in processes with elementary particles, but can serve as a dark matter particle. We argue that in very early Universe ($z>3000$) this vector field may have the equation of state of radiation ($w=1/3$) but later behaves as cold dark matter ($w=0$). This may slightly change the expansion rate of the Universe at early time and reduce the value of the sound horizon of baryon acoustic oscillations (standard ruler). As a result, in this model the value of the Hubble constant appears to be larger than that in the standard $\Lambda$CDM model. In particular, it is sufficient to have the dark photon mass of order $m\sim 10^{-27}-10^{-25}$ eV to fit the value of the Hubble constant to $H_0 = 73$ km$\cdot$s$^{-1}$Mpc$^{-1}$ thus resolving the Hubble tension.

Astrophysics today faces a number of mysteries which defiant their resolutions in spite of drastic improvements in instrumental design, better technique being developed, and gradual improvements in theoretical and computation methods over last decades. Primordial Lithium Puzzle is known to stay with us for at least two decades, and it is very likely that its final resolution will require some fundamentally new ideas, novel frameworks and a non-conventional paradigms. We propose that Primordial Lithium Puzzle finds its natural resolution within the so-called Axion Quark Nugget (AQN) dark matter model. This model was invented long ago as a natural explanation of the observed ratio $\Omega_{\rm dark} \sim \Omega_{\rm visible}$ without any references to BBN physics. In this new paradigm, in contrast with conventional WIMP (Weakly Interacting Massive Particle) framework the dark matter (DM) takes the form of the macroscopically large quark nuggets without requiring any new fields beyond the standard model physics, except for the axion. The time evolution of these AQNs in primordial soup at $T\sim 20$ KeV suggests a strong suppression of the abundances of nuclei with high charges $Z\geq 3$. This suppression mechanism represents the resolution of the primordial lithium abundance within AQN dark matter scenario.

White dwarf atmospheres are subjected to gravitational potentials around $10^5$ times larger than occur on Earth. They provide a unique environment in which to search for any possible variation in fundamental physics in the presence of strong gravitational fields. However, a sufficiently strong magnetic field will alter absorption line profiles and introduce additional uncertainties in measurements of the fine structure constant. Estimating the magnetic field strength is thus essential in this context. Here we model the absorption profiles of a large number of atomic transitions in the white dwarf photosphere, including first-order Zeeman effects in the line profiles, varying the magnetic field as a free parameter. We apply the method to a high signal-to-noise, high-resolution, far-ultraviolet HST/STIS spectrum of the white dwarf G191-B2B. The method yields a sensitive upper limit on its magnetic field of $B < 2300$ Gauss at the $3\sigma$ level. Using this upper limit we find that the potential impact of quadratic Zeeman shifts on measurements of the fine structure constant in G191-B2B is 4 orders of magnitude below laboratory wavelength uncertainties.

The Maxwell-Proca electrodynamics corresponding to a finite photon mass causes a substantial change of the Maxwell stress tensor and, under certain circumstances, may cause the electromagnetic stresses to act effectively as "negative pressure." The paper describes a model where this negative pressure imitates gravitational pull and may produce forces comparable to gravity and even become dominant. The effect is associated with the random magnetic fields in the galactic disk with a scale exceeding the photon Compton wavelength. The presence of a weaker regular field does not affect the forces under consideration. The stresses act predominantly on the interstellar gas and cause an additional force pulling the gas towards the center and towards the galactic plane. The stars do not experience any significant direct force but get involved in this process via a "recycling loop" where rapidly evolving massive stars are formed from the gas undergoing galactic rotation and then lose their masses back to the gas within a time shorter than roughly 1/6 of the rotation period. This makes their dynamics inseparable from that of the rotating gas. The lighter, slowly evolving stars, as soon as they are formed, lose connection to the gas and are confined within the galaxy only gravitationally. Numerical examples based on the parameters of our galaxy reveal both opportunities and challenges of this model and motivate further analysis. The critical issue is the plausibility of formation of the irregular magnetic field that would be force free. Another challenge is developing a predictive model of the evolution of the gaseous and stellar population of the galaxy under the aforementioned scenario. It may be interesting to also explore possible broader cosmological implications of the negative-pressure model.

Detection schemes for the quantum chromodynamics axions and other axion-like particles in light-shining-through-a-wall (LSW) experiments are based on the conversion of these particles into photons in a magnetic field. An alternative scheme may involve the detection via a resonant atomic or molecular transition induced by resonant axion absorption. The signal obtained in this process is second order in the axion-electron interaction constant but may become first order if we allow interference between the axion-induced transition amplitude and the transition amplitude induced by the electromagnetic radiation that produces the axions.

Various fundamental-physics experiments such as measurement of the birefringence of the vacuum, searches for ultralight dark matter (e.g., axions), and precision spectroscopy of complex systems (including exotic atoms containing antimatter constituents) are enabled by high-field magnets. We give an overview of current and future experiments and discuss the state-of-the-art DC- and pulsed-magnet technologies and prospects for future developments.

We derive a mean-field relativistic formula for the isotope shift of an electronic energy level for arbitrary angular momentum; we then use it to predict the spectra of superheavy metastable neutron-rich isotopes belonging to the hypothetical island of stability. Our results may be applied to the search for superheavy atoms in astrophysical spectra using the known values of the transition frequencies for the neutron deficient isotopes produced in the laboratory. An example of a relevant astrophysical system may be the spectra of the Przybylski's star where superheavy elements up to Z=99 have been possibly identified. In addition, it has been recently suggested to use the measurements of King plot non-linearity in a search for hypothetical new light bosons. On the other hand, one can find the non-linear corrections to the King-plot arising already in the Standard Model framework. We investigate contributions to the non-linearity arising from relativistic effects in the isotope field-shift, the nuclear polarizability and many-body effects. It is found that the nuclear polarizability contribution can lead to the significant deviation of the King plot from linearity. Therefore, the measurements of the non-linearity of King plots may be applied to obtain the nuclear polarizability change between individual isotopes. We then proceed with providing a rough analytical estimate of the non-linearity arising solely from the effect of a hypothetical scalar boson. Our predictions give theoretical limitations on the sensitivity of the search for new interactions and should help to identify the most suitable atoms for corresponding experiments.

We report on a search for ultra-low-mass axion-like dark matter by analysing the ratio of the spin-precession frequencies of stored ultracold neutrons and $^{199}$Hg atoms for an axion-induced oscillating electric dipole moment of the neutron and an axion-wind spin-precession effect. No signal consistent with dark matter is observed for the axion mass range $10^{-24}~\textrm{eV} \le m_a \le 10^{-17}~\textrm{eV}$. Our null result sets the first laboratory constraints on the coupling of axion dark matter to gluons, which improve on astrophysical limits by up to 3 orders of magnitude, and also improves on previous laboratory constraints on the axion coupling to nucleons by up to a factor of 40.

Hot white dwarf stars are the ideal probe for a relationship between the fine-structure constant and strong gravitational fields, providing us with an opportunity for a direct observational test. We study a sample of hot white dwarf stars, combining far-UV spectroscopic observations, atomic physics, atmospheric modelling and fundamental physics, in the search for variation in the fine structure constant. This variation manifests as shifts in the observed wavelengths of absorption lines, such as quadruply ionized iron (FeV) and quadruply ionized nickel (NiV), when compared to laboratory wavelengths. Berengut et al. (Phys. Rev. Lett. 2013, 111, 010801) demonstrated the validity of such an analysis using high-resolution Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) spectra of G191-B2B. We have made three important improvements by: (a) using three new independent sets of laboratory wavelengths, (b) analysing a sample of objects, and (c) improving the methodology by incorporating robust techniques from previous studies towards quasars (the Many Multiplet method). A successful detection would be the first direct measurement of a gravitational field effect on a bare constant of nature. Here we describe our approach and present preliminary results from nine objects using both FeV and NiV.

Spectral lines belonging to the short-lifetime heavy radioactive elements up to Es ($Z$=99) have been found in the spectra of the Przybylski's star. We suggest that these unstable elements may be decay products of a "magic" metastable nucleus belonging to the the island of stability where the nuclei have a magic number of neutrons $N=184$. The laboratory-produced nuclei have a significantly smaller number of neutrons. To identify spectra of the $N=184$ isotopes of these nuclei and their neutron-reach superheavy decay products in astrophysical data we calculate the isotope shift which should be added to the laboratory - measured wavelenghs. The results for the isotopic shifts in the strongest optical electromagnetic transitions in No, Lr, Nh, Fl,and $Z$=120 elements are presented.

Low-mass (sub-eV) spin-0 dark matter particles, which form a coherently oscillating classical field $\phi = \phi_0 \cos(m_\phi t)$, can induce oscillating variations in the fundamental constants through their interactions with the Standard Model sector. We calculate the effects of such possible interactions, which may include the linear interaction of $\phi$ with the Higgs boson, on atomic and molecular transitions. Using recent atomic clock spectroscopy measurements, we derive new limits on the linear interaction of $\phi$ with the Higgs boson, as well as its quadratic interactions with the photon and light quarks. For the linear interaction of $\phi$ with the Higgs boson, our derived limits improve on existing constraints by up to $2-3$ orders of magnitude.

In the comment of Avelino, Sousa and Lobo [arXiv:1506.06028], it is argued, by comparing the kinetic energy of a topological defect with the overall energy of a pulsar, that the origin of the pulsar glitch phenomenon due to the passage of networks of topological defects through pulsars is faced with serious difficulties. Here, we point out that topological defects may trigger pulsar glitches within traditional scenarios, such as vortex unpinning. If the energy transfer from a topological defect exceeds the activation energy for a single pinned vortex, this may lead to an avalanche of unpinning of vortices and consequently a pulsar glitch, and therefore the source of angular momentum and energy required for a glitch event is provided by the pulsar itself. Indeed, the activation energy for such a process can be very small (essentially zero compared with the observed increase in the pulsar's rotational kinetic energy at the onset of a glitch). The unpinning of a vortex by a topological defect may occur through the passage of the defect into the core of the pulsar.

We revisit the WIMP-type dark matter scattering on electrons that results in atomic ionization, and can manifest itself in a variety of existing direct-detection experiments. Unlike the WIMP-nucleon scattering, where current experiments probe typical interaction strengths much smaller than the Fermi constant, the scattering on electrons requires a much stronger interaction to be detectable, which in turn requires new light force carriers. We account for such new forces explicitly, by introducing a mediator particle with scalar or vector couplings to dark matter and to electrons. We then perform state of the art numerical calculations of atomic ionization relevant to the existing experiments. Our goals are to consistently take into account the atomic physics aspect of the problem (e.g., the relativistic effects, which can be quite significant), and to scan the parameter space: the dark matter mass, the mediator mass, and the effective coupling strength, to see if there is any part of the parameter space that could potentially explain the DAMA modulation signal. While we find that the modulation fraction of all events with energy deposition above 2 keV in NaI can be quite significant, reaching ~50%, the relevant parts of the parameter space are excluded by the XENON10 and XENON100 experiments.

We present an analysis of 23 absorption systems along the lines of sight towards 18 quasars in the redshift range of $0.4 \leq z_{abs} \leq 2.3$ observed on the Very Large Telescope (VLT) using the Ultraviolet and Visual Echelle Spectrograph (UVES). Considering both statistical and systematic error contributions we find a robust estimate of the weighted mean deviation of the fine-structure constant from its current, laboratory value of $\Delta\alpha/\alpha=\left(0.22\pm0.23\right)\times10^{-5}$, consistent with the dipole variation reported in Webb et al. and King et al. This paper also examines modelling methodologies and systematic effects. In particular we focus on the consequences of fitting quasar absorption systems with too few absorbing components and of selectively fitting only the stronger components in an absorption complex. We show that using insufficient continuum regions around an absorption complex causes a significant increase in the scatter of a sample of $\Delta\alpha/\alpha$ measurements, thus unnecessarily reducing the overall precision. We further show that fitting absorption systems with too few velocity components also results in a significant increase in the scatter of $\Delta\alpha/\alpha$ measurements, and in addition causes $\Delta\alpha/\alpha$ error estimates to be systematically underestimated. These results thus identify some of the potential pitfalls in analysis techniques and provide a guide for future analyses.

Atoms and molecules can become ionized during the scattering of a slow, heavy particle off a bound electron. Such an interaction involving leptophilic weakly interacting massive particles (WIMPs) is a promising possible explanation for the anomalous 9 sigma annual modulation in the DAMA dark matter direct detection experiment [R. Bernabei et al., Eur. Phys. J. C 73, 2648 (2013)]. We demonstrate the applicability of the Born approximation for such an interaction by showing its equivalence to the semiclassical adiabatic treatment of atomic ionization by slow-moving WIMPs. Conventional wisdom has it that the ionization probability for such a process should be exponentially small. We show, however, that due to nonanalytic, cusp-like behaviour of Coulomb functions close to the nucleus this suppression is removed, leading to an effective atomic structure enhancement. Crucially, we also show that electron relativistic effects actually give the dominant contribution to such a process, enhancing the differential cross section by up to 1000 times.

We propose to search for scalar dark matter via its effects on the electromagnetic fine-structure constant and particle masses. Scalar dark matter that forms an oscillating classical field produces `slow' linear-in-time drifts and oscillating variations of the fundamental constants, while scalar dark matter that forms topological defects produces transient-in-time variations of the constants of Nature. These variations can be sought for with atomic clock, laser interferometer and pulsar timing measurements. Atomic spectroscopy and Big Bang nucleosynthesis measurements already give improved bounds on the quadratic interaction parameters of scalar dark matter with the photon, electron, and light quarks by up to 15 orders of magnitude, while Big Bang nucleosynthesis measurements provide the first such constraints on the interaction parameters of scalar dark matter with the massive vector bosons.

We demonstrate that massive fields, such as dark matter, can directly produce a cosmological evolution of the fundamental constants of Nature. We show that a scalar or pseudoscalar (axion-like) dark matter field $\phi$, which forms a coherently oscillating classical field and interacts with Standard Model particles via quadratic couplings in $\phi$, produces `slow' cosmological evolution and oscillating variations of the fundamental constants. We derive limits on the quadratic interactions of $\phi$ with the photon, electron and light quarks from measurements of the primordial $^4$He abundance produced during Big Bang nucleosynthesis and recent atomic dysprosium spectroscopy measurements. These limits improve on existing constraints by up to 15 orders of magnitude. We also derive limits on the previously unconstrained linear and quadratic interactions of $\phi$ with the massive vector bosons from measurements of the primordial $^4$He abundance.

In the recent work [Phys. Rev. D 91, 111702(R) (2015)], C. Hill concludes that the axion electromagnetic anomaly induces an oscillating electron electric dipole moment of frequency $m_a$ and strength $\sim 10^{-32}~e$ cm, in the limit $v/c \to 0$ for the axion field. Here, we demonstrate that a proper treatment of this problem in the lowest order yields $\textit{no}$ electric dipole moment of the electron in the same limit. Instead, oscillating electric dipole moments of atoms and molecules are produced by different mechanisms.

We present an overview of recent developments in the detection of light bosonic dark matter, including axion, pseudoscalar axion-like and scalar dark matter, which form either a coherently oscillating classical field or topological defects (solitons). We emphasise new high-precision laboratory and astrophysical measurements, in which the sought effects are linear in the underlying interaction strength between dark matter and ordinary matter, in contrast to traditional detection schemes for dark matter, where the effects are quadratic or higher order in the underlying interaction parameters and are extremely small. New terrestrial experiments include measurements with atomic clocks, spectroscopy, atomic and solid-state magnetometry, torsion pendula, ultracold neutrons, and laser interferometry. New astrophysical observations include pulsar timing, cosmic radiation lensing, Big Bang nucleosynthesis and cosmic microwave background measurements. We also discuss various recently proposed mechanisms for the induction of slow `drifts', oscillating variations and transient-in-time variations of the fundamental constants of Nature by dark matter, which offer a more natural means of producing a cosmological evolution of the fundamental constants compared with traditional dark energy-type theories, which invoke a (nearly) massless underlying field. Thus, measurements of variation of the fundamental constants gives us a new tool in dark matter searches.

We present a brief overview of a new generation of high-precision laboratory and astrophysical measurements to search for ultralight (sub-eV) axion, axion-like pseudoscalar and scalar dark matter, which form either a coherent condensate or topological defects (solitons). In these new detection methods, the sought effects are linear in the interaction constant between dark matter and ordinary matter, which is in stark contrast to traditional searches for dark matter, where the sought effects are quadratic or higher order in the underlying interaction constants (which are extremely small).

Scalar dark matter can interact with Standard Model (SM) particles, altering the fundamental constants of Nature in the process. Changes in the fundamental constants during and prior to Big Bang nucleosynthesis (BBN) produce changes in the primordial abundances of the light elements. By comparing the measured and calculated (within the SM) primordial abundance of $^{4}$He, which is predominantly determined by the ratio of the neutron-proton mass difference to freeze-out temperature at the time of weak interaction freeze-out prior to BBN, we are able to derive stringent constraints on the mass of a scalar dark matter particle $\phi$ together with its interactions with the photon, light quarks and massive vector bosons via quadratic couplings in $\phi$, as well as its interactions with massive vector bosons via linear couplings in $\phi$. We also derive a stringent constraint on the quadratic interaction of $\phi$ with the photon from recent atomic dysprosium spectroscopy measurements.

Any slight variations in the fundamental constants of Nature, which may be induced by dark matter or some yet-to-be-discovered cosmic field, would characteristically alter the phase of a light beam inside an interferometer, which can be measured extremely precisely. Laser and maser interferometry may be applied to searches for the linear-in-time drift of the fundamental constants, detection of topological defect dark matter through transient-in-time effects and for a relic, coherently oscillating condensate, which consists of scalar dark matter fields, through oscillating effects. Our proposed experiments require either minor or no modifications of existing apparatus, and offer extensive reach into important and unconstrained spaces of physical parameters.

We show that the interaction of an axion field, or in general a pseudoscalar field, with the axial-vector current generated by an electron through a derivative-type coupling can give rise to a time-dependent mixing of opposite-parity states in atomic and molecular systems. Likewise, the analogous interaction of an axion field with the axial-vector current generated by a nucleon can give rise to time-dependent mixing of opposite-parity states in nuclear systems. This mixing can induce oscillating electric dipole moments, oscillating parity non-conservation effects and oscillating anapole moments in such systems. By adjusting the energy separation between the opposite-parity states of interest to match the axion mass energy, axion-induced experimental observables can be enhanced by many orders of magnitude. Oscillating atomic electric dipole moments can also be generated by axions through hadronic mechanisms, namely the P,T-violating nucleon-nucleon interaction and through the axion-induced electric dipole moments of valence nucleons, which comprise the nuclei. The axion field is modified by the Earth's gravitational field. The interaction of the spin of either an electron or nucleon with this modified axion field leads to axion-induced observable effects. These effects, which are of the form $\mathbf{g} \cdot \mathbf{\sigma}$, differ from the axion-wind effect, which has the form $\mathbf{p}_{\textrm{a}} \cdot \mathbf{\sigma}$.

We calculate the proton and neutron spin contributions for nuclei using semi-empirical methods, as well as a novel hybrid \emph{ab initio}/semi-empirical method, for interpretation of experimental data. We demonstrate that core-polarisation corrections to \emph{ab initio} nuclear shell model calculations generally reduce discrepancies in proton and neutron spin expectation values from different calculations. We derive constraints on the spin-dependent P,T-violating interaction of a bound proton with nucleons, which for certain ranges of exchanged pseudoscalar boson masses improve on the most stringent laboratory limits by several orders of magnitude. We derive a limit on the CPT and Lorentz-invariance-violating parameter $|\tilde{b}_{\perp}^p| < 7.6 \times 10^{-33}$ GeV, which improves on the most stringent existing limit by a factor of 8, and demonstrate sensitivities to the parameters $\tilde{d}_{\perp}^p$ and $\tilde{g}_{ D\perp}^p$ at the level $\sim 10^{-29} - 10^{-28}$ GeV, which is a one order of magnitude improvement compared to the corresponding existing sensitivities. We extend previous analysis of nuclear anapole moment data for Cs to obtain new limits on several other CPT and Lorentz-invariance-violating parameters: $\left|b_0^p \right| < 7 \times 10^{-8}$ GeV, $\left|d_{00}^p \right| < 8 \times 10^{-8}$, $\left|b_0^n \right| < 3 \times 10^{-7}$ GeV and $\left|d_{00}^n \right| < 3 \times 10^{-7}$.

We propose methods and present calculations that can be used to search for evidence of cosmic fields by investigating the parity-violating effects, including parity nonconservation amplitudes and electric dipole moments, that they induce in atoms. The results are used to constrain important fundamental parameters describing the strength of the interaction of various cosmic fields with electrons, protons, and neutrons. Candidates for such fields are dark matter (including axions) and dark energy, as well as several more exotic sources described by standard-model extensions. Existing parity nonconservation experiments in Cs, Dy, Yb, and Tl are combined with our calculations to directly place limits on the interaction strength between the temporal component, b_0, of a static pseudovector cosmic field and the atomic electrons, with the most stringent limit of |b_0^e| < 7*10^(-15) GeV, in the laboratory frame of reference, coming from Dy. From a measurement of the nuclear anapole moment of Cs, and a limit on its value for Tl, we also extract limits on the interaction strength between the temporal component of this cosmic field, as well as a related tensor cosmic-field component d_00, with protons and neutrons. The most stringent limits of |b_0^p| < 4*10^(-8) GeV and |d_00^p| < 5*10^(-8) for protons, and |b_0^n| < 2*10^(-7) GeV and |d_00^n| < 2*10^(-7) for neutrons (in the laboratory frame) come from the results using Cs. Axions may induce oscillating P- and T-violating effects in atoms and molecules through the generation of oscillating nuclear magnetic quadrupole and Schiff moments, which arise from P- and T-odd intranuclear forces and from the electric dipole moments of constituent nucleons. Nuclear-spin-independent parity nonconservation effects may be enhanced in diatomic molecules possessing close pairs of opposite-parity levels in the presence of time-dependent interactions.

We present a brief overview of recently proposed detection schemes for axion, axion-like pseudoscalar particle and topological defect dark matter. We focus mainly on the possibility of using atomic and molecular systems for dark matter detection. For axions and axion-like particles, these methods are complementary probes to ongoing photon-axion interconversion experiments and astrophysical observations. For topological defects, these methods are complementary to conventional astrophysical search schemes based on gravitational signatures.

We propose schemes for the detection of topological defect dark matter using pulsars and other luminous extraterrestrial systems via non-gravitational signatures. The dark matter field, which makes up a defect, may interact with standard model particles, including quarks and the photon, resulting in the alteration of their masses. When a topological defect passes through a pulsar, its mass, radius and internal structure may be altered, resulting in a pulsar `quake'. A topological defect may also function as a cosmic dielectric material with a distinctive frequency-dependent index of refraction, which would give rise to the time delay of a periodic extraterrestrial light or radio signal, and the dispersion of a light or radio source in a manner distinct to a gravitational lens. A topological defect passing through Earth may give rise to temporary non-zero electric dipole moments for an electron, proton, neutron, nuclei and atoms. The biggest advantage of such astrophysical observations over recently proposed terrestrial detection methods is the much higher probability of a defect been found in the vast volumes of outer space compared with one passing through Earth itself.

We propose methods for extracting limits on the strength of P-odd interactions of pseudoscalar and pseudovector cosmic fields with electrons, protons and neutrons. Candidates for such fields are dark matter (including axions) and dark energy, as well as several more exotic sources described by standard-model extensions. Calculations of parity nonconserving amplitudes and atomic electric dipole moments induced by these fields are performed for H, Li, Na, K, Rb, Cs, Ba+, Tl, Dy, Fr, and Ra+. From these calculations and existing measurements in Dy, Cs and Tl, we constrain the interaction strengths of the parity-violating static pseudovector cosmic field to be 7*10^(-15) GeV with an electron, and 3*10^(-8) GeV with a proton.

We propose a new probe of the dependence of the fine structure constant, alpha, on a strong gravitational field using metal lines in the spectra of white dwarf stars. Comparison of laboratory spectra with far-UV astronomical spectra from the white dwarf star G191-B2B recorded by the Hubble Space Telescope Imaging Spectrograph gives limits on the fractional variation of alpha of (Delta alpha/alpha)=(4.2 +- 1.6)x10^(-5) and (-6.1 +- 5.8)x10^(-5) from Fe V and Ni V spectra, respectively, at a dimensionless gravitational potential relative to Earth of (Delta phi) ~ 5x10^(-5). With better determinations of the laboratory wavelengths of the lines employed these results could be improved by up to two orders of magnitude.

The quark mass dependences of light element binding energies and nuclear scattering lengths are derived using chiral perturbation theory in combination with non-perturbative methods. In particular, we present new, improved values for the quark mass dependence of meson resonances that enter the nuclear force. A detailed analysis of the theoretical uncertainties arising in this determination is presented. As an application we derive from a comparison of observed and calculated primordial deuterium and helium abundances a stringent limit on the variation of the light quark mass, $\delta m_q/m_q = 0.02 \pm 0.04$. Inclusion of the neutron lifetime modification under the assumption of a variation of the Higgs vacuum expectation value that translates into changing quark, electron, and weak gauge boson masses, leads to a stronger limit, $|\delta m_q/m_q| < 0.009$.

We statistically analyse a recent sample of data points measuring the fine-structure constant alpha (relative to the terrestrial value) in quasar absorption systems. Using different statistical techniques, we find general agreement with previous authors that a dipole model is a well-justified fit to the data. We determine the significance of the dipole fit relative to that of a simple monopole fit, discuss the consistency of the interpretation, and test alternate models for potential variation of alpha against the data. Using a simple analysis we find that the monopole term (the constant offset in (delta alpha)/alpha) may be caused by non-terrestrial magnesium isotope abundances in the absorbers. Finally we test the domain-wall model against the data.

In their Letter, Kentosh and Mohageg [Phys. Rev. Lett. 108, 110801 (2012)] seek to use data from clocks aboard global positioning system (GPS) satellites to place limits on local position invariance (LPI) violations of Planck's constant, h. It is the purpose of this comment to show that discussing limits on variation of dimensional constants (such as h) is not meaningful; that even within a correct framework it is not possible to extract limits on variation of fundamental constants from a single type of clock aboard GPS satellites; and to correct an important misconception in the authors' interpretation of previous Earth-based LPI experiments.

(abridged) We present a new analysis of a large sample of quasar absorption-line spectra obtained using UVES (the Ultraviolet and Visual Echelle Spectrograph) on the VLT (Very Large Telescope) in Chile. In the VLT sample (154 absorbers), we find evidence that alpha increases with increasing cosmological distance from Earth. However, as previously shown, the Keck sample (141 absorbers) provided evidence for a smaller alpha in the distant absorption clouds. Upon combining the samples an apparent variation of alpha across the sky emerges which is well represented by an angular dipole model pointing in the direction RA=(17.3 +/- 1.0) hr, dec. = (-61 +/- 10) deg, with amplitude (0.97 +0.22/-0.20) x 10^(-5). The dipole model is required at the 4.1 sigma statistical significance level over a simple monopole model where alpha is the same across the sky (but possibly different to the current laboratory value). The data sets reveal a number of remarkable consistencies: various data cuts are consistent and there is consistency in the overlap region of the Keck and VLT samples. Assuming a dipole-only (i.e. no-monopole) model whose amplitude grows proportionally with `lookback-time distance' (r=ct, where t is the lookback time), the amplitude is (1.1 +/- 0.2) x 10^(-6) GLyr^(-1) and the model is significant at the 4.2 sigma confidence level over the null model [Delta alpha]/alpha = 0). We apply robustness checks and demonstrate that the dipole effect does not originate from a small subset of the absorbers or spectra. We present an analysis of systematic effects, and are unable to identify any single systematic effect which can emulate the observed variation in alpha.

We previously reported Keck telescope observations suggesting a smaller value of the fine structure constant, alpha, at high redshift. New Very Large Telescope (VLT) data, probing a different direction in the universe, shows an inverse evolution; alpha increases at high redshift. Although the pattern could be due to as yet undetected systematic effects, with the systematics as presently understood the combined dataset fits a spatial dipole, significant at the 4.2-sigma level, in the direction right ascension 17.5 +/- 0.9 hours, declination -58 +/- 9 degrees. The independent VLT and Keck samples give consistent dipole directions and amplitudes, as do high and low redshift samples. A search for systematics, using observations duplicated at both telescopes, reveals none so far which emulate this result.

A possibility to constrain axion-like particles from precision atomic physics is considered.

We compare the redshifts of neutral carbon and carbon monoxide in the redshifted sources in which the fine structure transition of neutral carbon, [CI], has been detected, in order to measure space-time variation of the fundamental constants. Comparison with the CO rotational lines measures gives the same combination of constants obtained from the comparison fine structure line of singly ionised carbon, [CII]. However, neutral carbon has the distinct advantage that it may be spatially coincident with the carbon monoxide, whereas [CII] could be located in the diffuse medium between molecular clouds, and so any comparison with CO could be dominated by intrinsic velocity differences. Using [CI], we obtain a mean variation of dF/F = (-3.6 +/- 8.5) x 10^-5, over z = 2.3 - 4.1, for the eight [CI] systems, which degrades to (-1.5+/- 11) x 10^-5, over z = 2.3 - 6.4 when the two [CII] systems are included. That is, zero variation over look-back times of 10.8-12.8 Gyr. However, the latest optical results indicate a spatial variation in alpha, which describes a dipole and we see the same direction in dF/F. This trend is, however, due to a single source for which the [CI] spectrum is of poor quality. This also applies to one of the two [CII] spectra previously used to find a zero variation in alpha^2/mu. Quantifying this, we find an anti-correlation between |dF/F| and the quality of the carbon detection, as measured by the spectral resolution, indicating that the typical values of >50 km/s, used to obtain a detection, are too coarse to reliably measure changes in the constants. From the fluxes of the known z > 1 CO systems, we predict that current instruments are incapable of the sensitivities required to measure changes in the constants through the comparison of CO and carbon lines. We therefore discuss in detail the use of ALMA for such an undertaking ... ABRIDGED

We identify the molecular ion NH^+ as a potential candidate for probing variations in the fine structure constant alpha and electron-to-proton mass ratio mu. NH^+ has an anomalously low-lying excited Sigma state, being only a few hundred cm^-1 above the ground Pi state. Being a light molecule, this proximity is such that rotational levels of the respective states are highly intermixed for low angular momenta. We find that several low-frequency transitions within the collective rotational spectrum experience enhanced sensitivity to alpha- and mu-variation. This is attributable to the close proximity of the Pi and Sigma states, as well as the ensuing strong spin-orbit coupling between them. Suggestions that NH^+ may exist in interstellar space and recent predictions that trapped-ion precision spectroscopy will be adaptable to molecular ions make NH^+ a promising system for future astrophysical and laboratory studies of alpha- and mu-variation.

Theories unifying gravity with other interactions suggest spatial and temporal variation of fundamental "constants" in the Universe. A change in the fine structure constant, alpha, could be detected via shifts in the frequencies of atomic transitions in quasar absorption systems. Recent studies using 140 absorption systems from the Keck telescope and 153 from the Very Large Telescope, suggest that alpha varies spatially. That is, in one direction on the sky alpha seems to have been smaller at the time of absorption, while in the opposite direction it seems to have been larger. To continue this study we need accurate laboratory measurements of atomic transition frequencies. The aim of this paper is to provide a compilation of transitions of importance to the search for alpha variation. They are E1 transitions to the ground state in several different atoms and ions, with wavelengths ranging from around 900 - 6000 A, and require an accuracy of better than 10^{-4} A. We discuss isotope shift measurements that are needed in order to resolve systematic effects in the study. The coefficients of sensitivity to alpha-variation (q) are also presented.

We review recent activity searching for variations in the fundamental constants of nature in quasar absorption spectra and in the laboratory. While research in this direction has been ongoing for many decades, the topic has recently been stimulated by astronomical evidence for spatial variation in the fine-structure constant, alpha. This result could be confirmed using different quasar data and atomic clock measurements, but there are significant challenges to obtain the required accuracy. We review existing measurements and discuss some of the most promising systems where any variations would be strongly enhanced.

The detection of a spatial variation of the fine-structure constant, alpha, based on study of quasar absorption systems has recently been reported. The physics that causes this alpha-variation should have other observable manifestations, and this motivates us to look for complementary astrophysical effects. In this paper we propose a method to test whether spatial variation of fundamental constants existed during the epoch of big bang nucleosynthesis. Using existing measurements of primordial deuterium abundance we find very weak indications that such a signature might exist, but the paucity of measurements precludes any firm conclusion. We also examine existing quasar absorption spectra data that are sensitive to variation of the electron-to-proton mass ratio, mu, and x = (alpha^2 mu g_p) for spatial variation.

The remarkable detection of a spatial variation in the fine-structure constant, alpha, from quasar absorption systems must be independently confirmed by complementary searches. In this letter, we discuss how terrestrial measurements of time-variation of the fundamental constants in the laboratory, meteorite data, and analysis of the Oklo nuclear reactor can be used to corroborate the spatial variation seen by astronomers. Furthermore, we show that spatial variation of the fundamental constants may be observable as spatial anisotropy in the cosmic microwave background, the accelerated expansion (dark energy), and large-scale structure of the Universe.

We consider a possibility of capture of a heavy charged massive particle $\chi^-$ by the nucleus leading to appearance of a bound state. A simple analytic formula allowing to calculate binding energies of the $N\chi^-$ bound state for different nuclei is derived. If the binding energy is sufficiently large the particle $\chi^-$ is stable inside the nucleus. The probabilities of the nuclear fissions for such states are calculated. It is shown that the bound states are more stable to a possible fission in comparison to the bare nucleus. This makes an observation of this hypothetical charged massive particle and the superheavy nuclei more probable.

Khatri and Wandelt reported that change in the value of alpha by 1% changes the mean brightness temperature $T_b$ decrement of the CMB due to 21 cm absorption by 5% over the redshift range z $<$ 50. A drawback of their approach is that the dimensionful parameters are used. Changing of units leads to the change of the magnitude and even sign of the effect. Similar problems may be identified in a large number of other publications which consider limits on the variation of alpha using dimentionful parameters. We propose a method to obtain consistent results and provide an estimate of the effect.

Observing inversion lines of ammonia (NH3), complemented by rotational lines of NH3 and other molecular species, provides stringent constraints on potential variations of the proton-to-electron mass ratio. While a limit of one part per million is derived for a lookback time of 7 billion years, nearby dark clouds might show a statistically significant variation of order 20-30 parts per billion, possibly being related to chameleon fields. The detection of radio-loud quasars with strong molecular absorption lines at redshifts z > 1 as well as the identification of a larger sample of nearby dark clouds with exceptionally narrow lines (<0.2 km/s) would be essential to improve present limits and to put the acquired results onto a firmer statistical basis.

Using the relativistic Hartree-Fock approximation, we calculate the rates of atomic ionization by absorption of pseudoscalar particles in the mass range from 10 to $\sim$ 50 keV. We present numerical results for atoms relevant for the direct dark matter searches (e.g. Ar, Ge, I and Xe), as well as the analytical formula which fits numerical calculations with few per cent accuracy and may be used for multi-electron atoms, molecules and condensed matter systems.

We use relativistic Hartree-Fock and configuration interaction methods to calculate the dependence of transition frequencies for singly ionized cobalt on the fine structure constant. The results are to be used in the search for variation of the fine structure constant in quasar absorption spectra.

The Keck telescope's HIRES spectrograph has previously provided evidence for a smaller fine-structure constant, alpha, compared to the current laboratory value, in a sample of 143 quasar absorption systems: da/a=(-0.57+/-0.11)x10^{-5}. This was based on a variety of metal-ion transitions which, if alpha varies, experience different relative velocity shifts. This result is yet to be robustly contradicted, or confirmed, by measurements on other telescopes and spectrographs; it remains crucial to do so. It is also important to consider new possible instrumental systematic effects which may explain the Keck/HIRES results. Griest et al. (2009, arXiv:0904.4725v1) recently identified distortions in the echelle order wavelength scales of HIRES with typical amplitudes +/-250m/s. Here we investigate the effect such distortions may have had on the Keck/HIRES varying alpha results. We demonstrate that they cause a random effect on da/a from absorber to absorber because the systems are at different redshifts, placing the relevant absorption lines at different positions in different echelle orders. The typical magnitude of the effect on da/a is ~0.4x10^{-5} per absorber which, compared to the median error on da/a in the sample, ~1.9x10^{-5}, is relatively small. Consequently, the weighted mean value changes by less than 0.05x10^{-5} if the corrections we calculate are applied. Nevertheless, we urge caution, particularly for analyses aiming to achieve high precision da/a measurements on individual systems or small samples, that a much more detailed understanding of such intra-order distortions and their dependence on observational parameters is important if they are to be avoided or modelled reliably. [Abridged]

We calculate the effect of variation in the light-current quark mass, $m_q$, on standard big bang nucleosynthesis. A change in $m_q$ at during the era of nucleosynthesis affects nuclear reaction rates, and hence primordial abundances, via changes the binding energies of light nuclei. It is found that a relative variation of $\delta m_q/m_q = 0.016 \pm 0.005$ provides better agreement between observed primordial abundances and those predicted by theory. This is largely due to resolution of the existing discrepancies for 7Li. However this method ignores possible changes in the position of resonances in nuclear reactions. The predicted 7Li abundance has a strong dependence on the cross-section of the resonant reactions 3He(d,p)4He and t(d,n)4He. We show that changes in $m_q$ at the time of BBN could shift the position of these resonances away from the Gamow window and lead to an increased production of 7Li, exacerbating the lithium problem.

Based on measurements with the Effelsberg 100-m telescope, a multi-line study of molecular species is presented toward the gravitational lens system PKS1830-211. Obtaining average radial velocities and performing Large Velocity Gradient radiative transfer calculations, the aims of this study are (1) to determine the density of the gas, (2) to constrain the temperature of the cosmic microwave background, and (3) to evaluate the proton-to-electron mass ratio at redshift 0.9. Analyzing data from six rotational HC3N transitions (this includes the J=7-6 line, which is likely detected for the first time in the interstellar medium) we obtain about 2000 cm-3 for the gas density of the south-western absorption component. Again toward the south-western source, excitation temperatures of molecular species with optically thin lines and higher rotational constants are, on average, consistent with the expected temperature of the cosmic microwave background, T_CMB = 5.14 K. However, individually, there is a surprisingly large scatter which far surpasses expected uncertainties. A comparison of CS J=1-0 and 4-3 optical depths toward the weaker north-western absorption component results in an excitation temperature of 11 K and a 1-sigma error of 3 K. For the south-eastern main component, a comparison of velocities determined from ten optically thin NH3 inversion lines with those from five optically thin rotational transitions of HC3N, observed at similar frequencies, constrains potential variations of the proton-to-electron mass ratio, with respect to its present value, to <1.4 x 10^-6 with 3-sigma confidence. Also including optically thin rotational lines from other molecular species, it is emphasized that systematic errors are smaller than 1 km/s, corresponding to an uncertainty of 10-6.

In the Letter [1] (also [2]) there is a claim that the generalised second law of thermodynamics (entropy increase) for black holes provides some limits on the rate of variation of the fundamental constants of nature (electric charge e, speed of light c, etc.). We have come to a different conclusion. The results in [1,2] are based on assumption that mass of a black hole does not change without radiation and accreation. We present arguments showing that this assumption is incorrect and give an estimate of the black hole mass variation due to alpha=e^2/\hbar c variation using entropy (and quantum energy level) conservation in an adiabatic process. No model-independent limits on the variation of the fundamental constants are derived from the second law of thermodynamics.

The Standard Model of particle physics assumes that the so-called fundamental constants are universal and unchanging. Absorption lines arising in molecular clouds along quasar sightlines offer a precise test for variations in the proton-to-electron mass ratio, mu, over cosmological time and distance scales. The inversion transitions of ammonia are particularly sensitive to mu compared to molecular rotational transitions. Comparing the available ammonia spectra observed towards the quasar B0218+357 with new, high-quality rotational spectra, we present the first detailed measurement of mu with this technique, limiting relative deviations from the laboratory value to |dmu/mu| < 1.8x10^{-6} (95% confidence level) at approximately half the Universe's current age - the strongest astrophysical constraint to date. Higher-quality ammonia observations will reduce both the statistical and systematic uncertainties in these measurements.

A calculation of the current-quark-mass-dependence of nucleon static electromagnetic properties is necessary in order to use observational data as a means to place constraints on the variation of Nature's fundamental parameters. A Poincare' covariant Faddeev equation, which describes baryons as composites of confined-quarks and -nonpointlike-diquarks, is used to calculate this dependence The results indicate that, like observables dependent on the nucleons' magnetic moments, quantities sensitive to their magnetic and charge radii, such as the energy levels and transition frequencies in Hydrogen and Deuterium, might also provide a tool with which to place limits on the allowed variation in Nature's constants.

In their Letter [Phys. Rev. Lett. 92, 121302 (2004)] (also [Astron. Astrophys. 417, 853 (2004)]), Srianand et al. analysed optical spectra of heavy-element species in 23 absorption systems along background quasar sight-lines, reporting limits on relative variations in the fine-structure constant: da/a=(-0.06+/-0.06) x 10^{-5}. Here we demonstrate basic flaws in their analysis, using the same data and absorption profile fits, which led to spurious values of da/a and significantly underestimated uncertainties. We conclude that these data and fits offer no stringent test of previous evidence for a varying alpha. In their Reply (arXiv:0711.1742) to this Comment, Srianand et al. state or argue several points regarding their original analysis and our new analysis. We discuss these points here, dismissing all of them because they are demonstrably incorrect or because they rely on a flawed application of simple statistical arguments.

We critically review the current null results on a varying fine-structure constant, alpha, derived from VLT/UVES quasar absorption spectra, focusing primarily on the many-multiplet analysis of 23 absorbers from which Chand et al. (2004) reported a weighted mean relative variation of da/a=(-0.06+/-0.06) x 10^{-5}. Our analysis of the same reduced data, using the same fits to the absorption profiles, yields very different individual da/a values with uncertainties typically larger by a factor of ~3. We attribute the discrepancies to flawed parameter estimation techniques in the original analysis and demonstrate that the original da/a values were strongly biased towards zero. Were those flaws not present, the input data and spectra should have given a weighted mean of da/a=(-0.44+/-0.16) x 10^{-5}. Although this new value does reflect the input spectra and fits (unchanged from the original work -- only our analysis is different), we do not claim that it supports previous Keck/HIRES evidence for a varying alpha: there remains significant scatter in the individual da/a values which may stem from the overly simplistic profile fits in the original work. Allowing for such additional, unknown random errors by increasing the uncertainties on da/a to match the scatter provides a more conservative weighted mean, da/a=(-0.64+/-0.36) x 10^{-5}. We highlight similar problems in other current UVES constraints on varying alpha and argue that comparison with previous Keck/HIRES results is premature.

We examine how the binding of light ($A\leq 8$) nuclei depends on possible variations of hadronic masses, including meson, nucleon, and nucleon-resonance masses. Small variations in hadronic masses may have occurred over time; the present results can help evaluate the consequences for big bang nucleosynthesis. Larger variations may be relevant to current attempts to extrapolate properties of nucleon-nucleon interactions from lattice QCD calculations. Results are presented as derivatives of the energy with respect to the different masses so they can be combined with different predictions of the hadronic mass-dependence on the underlying current-quark mass $m_q$. As an example, we employ a particular set of relations obtained from a study of hadron masses and sigma terms based on Dyson-Schwinger equations and a Poincar\'{e}-covariant Faddeev equation for confined quarks and diquarks. We find that nuclear binding decreases moderately rapidly as the quark mass increases, with the deuteron becoming unbound when the pion mass is increased by $\sim$60% (corresponding to an increase in $X_q=m_q/\Lambda_{QCD}$ of 2.5). In the other direction, the dineutron becomes bound if the pion mass is decreased by $\sim$15% (corresponding to a reduction of $X_q$ by $\sim$30%). If we interpret the disagreement between big bang nucleosynthesis calculations and measurements to be the result of variation in $X_q$, we obtain an estimate $\delta X_q/X_q=K \cdot (0.013 \pm 0.002)$ where $K \sim 1$ (the expected accuracy in $K$ is about a factor of 2). The result is dominated by $^7$Li data.

Radio-frequency E1 transitions between nearly degenerate, opposite parity levels of atomic dysprosium were monitored over an eight month period to search for a variation in the fine-structure constant. During this time period, data were taken at different points in the gravitational potential of the Sun. The data are fitted to the variation in the gravitational potential yielding a value of $(-8.7 \pm 6.6) \times 10^{-6}$ for the fit parameter $k_\alpha$. This value gives the current best laboratory limit. In addition, our value of $k_{\alpha}$ combined with other experimental constraints is used to extract the first limits on k_e and k_q. These coefficients characterize the variation of m_e/m_p and m_q/m_p in a changing gravitational potential, where m_e, m_p, and m_q are electron, proton, and quark masses. The results are $k_e = (4.9 \pm 3.9) \times 10^{-5}$ and $k_q = (6.6 \pm 5.2) \times 10^{-5}$.

Review of recent works devoted to the variation of the fundamental constants is presented including atomic clocks, quasar absorption spectra, and Oklo natural nuclear reactor data. Assuming linear variation with time we can compare different results. From the quasar absorption spectra: $\dot{\mu}/\mu=(1 \pm 3) \times 10^{-16}$ yr$^{-1}$. A combination of this result and the atomic clock results gives the best limt on variation of $\alpha$: $\dot{\alpha}/\alpha=(-0.8 \pm 0.8) \times 10^{-16}$ yr$^{-1}$. The Oklo natural reactor gives the best limit on the variation of $m_s/\Lambda_{QCD}$ where $m_s$ is the strange quark mass. Huge enhancement of the relative variation effects happens in transitions between close atomic, molecular and nuclear energy levels. We suggest several new cases where the levels are very narrow. Large enhancement of the variation effects is also possible in cold atomic and molecular collisions near Feshbach resonance. Massive bodies (stars or galaxies) can also affect physical constants. They have large scalar charge $S$ proportional to number of particles which produces a Coulomb-like scalar field $U=S/r$. This leads to a variation of the fundamental constants proportional to the gravitational potential, e.g. $\delta \alpha/ \alpha = k_\alpha \delta (GM/ r c^2)$. We compare different manifestations of this effect.The strongest limit $k_\alpha +0.17 k_e= (-3.5\pm 6) \times 10^{-7}$.

We calculate the sensitivity of the inversion spectrum of ammonia to possible time-variation of the ratio of the proton mass to the electron mass, mu=m_p/m_e. For the inversion transition (lambda= 1.25 cm^{-1}) the relative frequency shift is significantly enhanced: delta(omega)/omega=-4.46, delta(mu)/mu. This enhancement allows one to increase sensitivity to the time-variation of mu using NH_3 spectra for high redshift objects. We use published data on microwave spectra of the object B0218+357 to place the limit delta(mu)/mu =(0.6 +/- 1.9) 10^{-6} at redshift z=0.6847; this limit is several times better than the limits obtained by different methods and may be significantly improved. Assuming linear time dependence we obtain dot{mu}/mu=(-1 +/- 3) 10^{-16} yr^{-1}.

Light scalar fields very naturally appear in modern cosmological models, affecting such parameters of Standard Model as electromagnetic fine structure constant $\alpha$, dimensionless ratios of electron or quark mass to the QCD scale, $m_{e,q}/\Lambda_{QCD}$. Cosmological variations of these scalar fields should occur because of drastic changes of matter composition in Universe: the latest such event is rather recent (redshift $z\sim 0.5$), from matter to dark energy domination. In a two-brane model (we use as a pedagogical example) these modifications are due to changing distance to "the second brane", a massive companion of "our brane". Back from extra dimensions, massive bodies (stars or galaxies) can also affect physical constants. They have large scalar charge $Q_d$ proportional to number of particles which produces a Coulomb-like scalar field $\phi=Q_d/r$. This leads to a variation of the fundamental constants proportional to the gravitational potential, e.g. $\delta \alpha/ \alpha = k_\alpha \delta (GM/ r c^2)$. We compare different manifestations of this effect. The strongest limits $k_\alpha +0.17 k_e= (-3.5\pm 6) * 10^{-7}$ are obtained from the measurements of dependence of atomic frequencies on the distance from Sun (the distance varies due to the ellipticity of the Earth's orbit).

Current analyses of VLT/UVES quasar spectra disagree with the Keck/HIRES evidence for a varying fine-structure constant, alpha. To investigate this we introduce a simple method for calculating the minimum possible uncertainty on Delta(alpha)/alpha for a given quasar absorber. For many absorbers in Chand et al. (2004) and for the single-absorber constraint of Levshakov et al. (2006) the quoted uncertainties are smaller than the minimum allowed by the UVES data. Failure of this basic consistency test prevents reliable comparison of the UVES and HIRES results.

Nine quasar absorption spectra at 21-cm and UV rest-wavelengths are used to estimate possible variations in x=alpha^2 g_p mu, (alpha is the fine structure constant, g_p the proton g-factor and mu=me/mp the electron-to-proton mass ratio). We find <Delta x/x>^weighted_total(=Dxxwt)=(0.63+-0.99) 10^-5 over 0.23~<z_abs~<2.35 (2.7 to 10.5 Gyr, look-back time, t_lb). A linear best fit against t_lb, tied to Delta x/x=0 at z=0, gives (dot x)/x=(-0.6+-1.2) 10^-15 /yr. Our large sample demonstrates that intrinsic line-of-sight velocity differences between the 21-cm and UV absorption redshifts, (on average Delta_vlos~6km/s), with random sign and magnitude in each absorption system, limit our precision. Combining our Delta x/x measurement with absorption-line constraints on alpha-variation yields strong limits on the variation of mu. Our most conservative estimate, obtained by assuming no variations in alpha or g_p is Delta mu/mu(=Dmm)=Dxxwt. If we use only the four high-redshift absorbers in our sample, we obtain Dmm=(0.58+-1.95) 10^-5, which agrees (2sigma) with recent, more direct estimates from two absorption systems containing molecular hydrogen, also at high redshift, and which have hinted at a possible mu-variation, Dmm=(-2.0+-0.6) 10^-5. Our method of constraining Dmm is completely independent from the molecular hydrogen observations. If we include the low-redshift systems, our Dmm result differs significantly from the high-redshift molecular hydrogen results. We detect a dipole variation in mu across the sky, but this model is required by the data at only the 88 per cent confidence level. Clearly, much larger samples of 21-cm and molecular hydrogen absorbers are required to adequately resolve the issue of the variation of mu and x.(Abridged)

We present a review of recent works devoted to the variation of the fine structure constant alpha, strong interaction and fundamental masses. There are some hints for the variation in quasar absorption spectra, Big Bang nucleosynthesis, and Oklo natural nuclear reactor data. A very promising method to search for the variation of the fundamental constants consists in comparison of different atomic clocks. Huge enhancement of the variation effects happens in transition between accidentally degenerate atomic and molecular energy levels. A new idea is to build a ``nuclear'' clock based on the ultraviolet transition between very low excited state and ground state in Thorium nucleus. This may allow to improve sensitivity to the variation up to 10 orders of magnitude! Huge enhancement of the variation effects is also possible in cold atomic and molecular collisions near Feschbach resonance.

We calculate the dependence of the nuclear magnetic moments on the quark masses including the spin-spin interaction effects and obtain limits on the variation of the fine structure constant $\alpha$ and $(m_q/\Lambda_{QCD})$ using recent atomic clock experiments examining hyperfine transitions in H, Rb, Cs, Yb$^+$ and Hg$^+$ and the optical transition in H, Hg$^+$ and Yb$^+$ .

We address the emerging discrepancy between the Big Bang Nucleosynthesis data and standard cosmology, which asks for a bit longer evolution time. If this effect is real, one possible implication (in a framework of brane cosmology model) is that there is a ``dark radiation'' component which is negative and makes few percents of ordinary matter density. If so, all scales of this model can be fixed, provided brane-to-bulk leakage problem is solved.

A calculation of the current-quark mass dependence of hadron masses can help in using observational data to place constraints on the variation of nature's fundamental parameters. A hadron's sigma-term is a measure of this dependence. The connection between a hadron's sigma-term and the Feynman-Hellmann theorem is illustrated with an explicit calculation for the pion using a rainbow-ladder truncation of the Dyson-Schwinger equations: in the vicinity of the chiral limit sigma_pi = m_pi/2. This truncation also provides a decent estimate of sigma_rho because the two dominant self-energy corrections to the rho-meson's mass largely cancel in their contribution to sigma_rho. The truncation is less accurate for the omega, however, because there is little to compete with an omega->rho+pi self-energy contribution that magnifies the value of sigma_omega by ~25%. A Poincare' covariant Faddeev equation, which describes baryons as composites of confined-quarks and -nonpointlike-diquarks, is solved to obtain the current-quark mass dependence of the masses of the nucleon and Delta, and thereby sigma_N and sigma_Delta. This "quark-core" piece is augmented by the "pion cloud" contribution, which is positive. The analysis yields sigma_N~60MeV and sigma_Delta~50MeV.

We present an accurate ab initio method of calculating isotope shifts and relativistic shifts in atomic spectra. We test the method on neutral carbon and three carbon ions. The relativistic shift of carbon lines may allow them to be included in analyses of quasar absorption spectra that seek to measure possible variations in the fine structure constant, alpha, over the lifetime of the Universe. Carbon isotope shifts can be used to measure isotope abundances in gas clouds: isotope abundances are potentially an important source of systematic error in the alpha-variation studies. These abundances are also needed to study nuclear reactions in stars and supernovae, and test models of chemical evolution of the Universe.

Quasar absorption spectra at 21-cm and UV rest-wavelengths are used to estimate the time variation of x=alpha^2 g_p mu, where alpha is the fine structure constant, g_p the proton g factor, and m_e/m_p=mu the electron/proton mass ratio. Over a redshift range 0.24 ~< z_abs ~< 2.04, <Delta x/x>_total_weighted=(1.17+-1.01) x 10^-5. A linear fit gives d/dt{x}/x=(-1.43+-1.27) x 10^-15 yr^-1. Two previous results on varying alpha yield the strong limits Delta mu/mu=(2.31+-1.03) x 10^-5 and Delta mu/mu=(1.29+-1.01) x 10^-5. Our sample, 8x larger than any previous, provides the first direct estimate of the intrinsic 21-cm and UV velocity differences ~6 km s^-1.

We present an ab initio method of calculation of isotope shift and relativistic shift in atoms with a few valence electrons. It is based on an energy calculation involving combination of the configuration interaction method and many-body perturbation theory. This work is motivated by analyses of quasar absorption spectra that suggest that the fine structure constant, alpha, was smaller at an early epoch. Relativistic shifts are needed to measure this variation of alpha, while isotope shifts are needed to resolve systematic effects in this study. The isotope shifts can also be used to measure isotopic abundances in gas clouds in the early universe, which are needed to study nuclear reactions in stars and supernovae and test models of chemical evolution. This paper shows that isotope shift in magnesium can be calculated to very high precision using our new method.

Sensitivity of atomic transition frequencies to variation of the fine structure constant $\alpha=e^2/\hbar c$ increases proportional to $Z^2$ where $Z$ is the nuclear charge. Recently several lines of heavy ions Ge II, Sn II and Pb II have been detected in quasar absorption spectra. We have performed accurate many-body calculations of $\alpha^2$-dependence of transition frequencies ($q$-coefficients) for these atoms and found an order of magnitude increase in sensitivity in comparison with atomic transitions which were previously used to search for temporal and spatial variation of $\alpha$ in quasar absorption systems. An interesting feature in Pb II is highly non-linear dependence on $\alpha^2$ due to the level pseudo-crossings.

We calculate the dependence of atomic transition frequencies on the fine structure constant, $\alpha = e^2/\hbar c$, for some ions of Ti, Mn, Na, C, and O. The results of these calculations will be used in the search for variation of alpha in quasar absorption spectra.

We show that using the modified form of the Dirac Hamiltonian as suggested by Bekenstein, does not affect the analysis in [1-4] of QSO data pertaining to a measurement of alpha variation. We obtain the present time limit on Bekenstein's parameter, tan^{2}chi =(0.2 \pm 0.7)\times10^{-6}, from the measurement of the hydrogen 2p fine structure using value of alpha obtained from different experiments.

We show why a dynamically varying fine structure constant does not impact on the Webb group's analysis as suggested by Bekenstein. We also provide a limit on the size of a possible correction to the Dirac Hamiltonian caused by a perturbation due to a dynamically varying alpha in hydrogen.

At present several groups are analyzing astrophysical data in a search for time-variation of the fine structure constant $\alpha$. Here we discuss how to exclude systematic effects caused by the changes in the isotope abundances during the evolution of the universe. We suggest using particular combinations of the transition frequencies of O II, Al II, Al III, Si II and Mn II as anchors, which are insensitive to $\alpha$-variation and to changes in isotope abundances. These anchors can be used to determine the cosmological redshift. Then, one can use other combinations of frequencies as probes for the time-variation of $\alpha$ and another set as probes for the isotopic abundances. In this way it is possible to eliminate one source of the systematic errors in the search for $\alpha$-variation and get new information about evolution of the isotopes. On the level of accuracy that has already been reached in the search for $\alpha$-variation it is possible to see ~10% changes in isotope abundances of Si and Ni.

Several recently found properties of the event horizon of black holes are discussed. One of them is the reflection of the incoming particles on the horizon. A particle approaching the black hole can bounce on the horizon back, into the outside world, which drastically reduces the absorption cross section in the infrared region. Another, though related phenomenon takes place for particles inside the horizon. A locked inside particle has, in fact, an opportunity to escape into the outside world. Thus, the confinement inside the horizon is not absolute. The escape from within the interior region of the horizon allows the transfer of information from this region into the outside world. This result may help resolve the information paradox for black holes. Both the reflection and escape phenomena happen due to pure quantum reasons, being impossible in the classical approximation.

We use the relativistic Hartree-Fock method, many-body perturbation theory and configuration-interaction method to calculate the dependence of atomic transition frequencies on the fine structure constant, alpha. The results of these calculations will be used in the search for variation of the fine structure constant in quasar absorption spectra.

We perform calculations of the dependence of nuclear magnetic moments on quark masses and obtain limits on the variation of $(m_q/\Lambda_{QCD})$ from recent measurements of hydrogen hyperfine (21 cm) and molecular rotational transitions in quasar absorption systems, atomic clock experiments with hyperfine transitions in H, Rb, Cs, Yb$^+$, Hg$^+$ and optical transition in Hg$^+$. Experiments with Cd$^+$, deuterium/hydrogen, molecular SF$_6$ and Zeeman transitions in $^3$He/Xe are also discussed.

The absorption cross section for scalar particle impact on a Schwarzschild black hole is found. The process is dominated by two physical phenomena. One of them is the well-known greybody factor that arises from the energy-dependent potential barrier outside the horizon that filters the incoming and outgoing waves. The other is related to the reflection of particles on the horizon (Kuchiev 2003). This latter effect strongly diminishes the cross section for low energies, forcing it to vanish in the infrared limit. It is argued that this is a general property, the absorption cross section vanishes in the infrared limit for scattering of particles of arbitrary spin.

We present evidence for variations in the fine-structure constant from Keck/HIRES spectra of 143 quasar absorption systems over the redshift range 0.2 < z_abs < 4.2. This includes 15 new systems, mostly at high-z (z_abs > 1.8). Our most robust estimate is a weighted mean da/a=(-0.57 +/- 0.11)x10^{-5}. We respond to recent criticisms of the many-multiplet method used to extract these constraints. The most important potential systematic error at low-z is the possibility of very different Mg heavy isotope abundances in the absorption clouds and laboratory: {\it higher} abundances of {25,26}Mg in the absorbers may explain the low-z results. Approximately equal mixes of {24}Mg and {25,26}Mg are required. Observations of Galactic stars generally show {\it lower} {25,26}Mg isotope fractions at the low metallicities typifying the absorbers. Higher values can be achieved with an enhanced population of intermediate mass stars at high redshift, a possibility at odds with observed absorption system element abundances. At present, all observational evidence is consistent with the varying-alpha results. Another promising method to search for variation of fundamental constants involves comparing different atomic clocks. Here we calculate the dependence of nuclear magnetic moments on quark masses and obtain limits on the variation of alpha and m_q/Lambda_QCD from recent atomic clock experiments with hyperfine transitions in H, Rb, Cs, Hg+ and an optical transition in Hg+.

We use Big Bang Nucleosynthesis calculations and light element abundance data to constrain the relative variation of the deuteron binding energy since the universe was a few minutes old, $\delta Q = Q(BBN)-Q(present)$. Two approaches are used, first treating the baryon to photon ratio, $\eta$, as a free parameter, but with the additional freedom of varying $\delta Q$, and second using the WMAP value of $\eta$ and solving only for $\delta Q$. Including varying $Q$ yields a better fit to the observational data than imposing the present day value, rectifying the discrepancy between the $^4He$ abundance and the deuterium and $^7Li$ abundances, {\it and} yields good agreement with the independently determined $\eta_{WMAP}$. The minimal deviation consistent with the data is significant at about the 4-$\sigma$ level; $\delta Q/Q= -0.019 \pm 0.005$. If the primordial $^4$He abundance lies towards the low end of values in the literature, this deviation is even larger and more statistically significant. Taking the light element abundance data at face-value, our result may be interpreted as variation of the dimensionless ratio $X=m_s/\Lambda_{QCD}$ of the strange quark mass and strong scale: $\delta X/X=(1.1 \pm 0.3) \times 10^{-3}$. These results provide a strong motivation for a more thorough exploration of the potential systematic errors in the light element abundance data.

[Abridged] We previously presented evidence for a varying fine-structure constant, alpha, in two independent samples of Keck/HIRES QSO spectra. Here we present a detailed many-multiplet analysis of a third Keck/HIRES sample containing 78 absorption systems. We also re-analyse the previous samples, providing a total of 128 absorption systems over the redshift range 0.2<z_abs<3.7. All three samples separately yield consistent, significant values of da/a. The analyses of low- and high-z systems rely on different ions/transitions with very different dependencies on alpha, yet they also give consistent results. We identify additional random errors in 22 high-z systems characterized by transitions with a large dynamic range in apparent optical depth. Increasing the statistical errors on da/a for these systems gives our fiducial result, a weighted mean da/a=(-0.543+/-0.116)x10^-5, representing 4.7-sigma evidence for a smaller weighted mean alpha in the absorption clouds. Assuming that da/a=0 at z_abs=0, the data marginally prefer a linear increase in alpha with time: dota/a=(6.40+/-1.35)x10^-16 yr^-1. The two-point correlation function for alpha is consistent with zero over 0.2-13 Gpc comoving scales and the angular distribution of da/a shows no significant dipolar anisotropy. We therefore have no evidence for spatial variations in da/a. We extend our previous searches for possible systematic errors, identifying atmospheric dispersion and isotopic structure effects as potentially the most significant. However, overall, known systematic errors do not explain the results. Future many-multiplet analyses of QSO spectra from different telescopes and spectrographs will provide a now crucial check on our Keck/HIRES results.

This work presents a method for the ab initio calculation of isotope shift in atoms and ions with one valence electron above closed shells. As a zero approximation we use relativistic Hartree-Fock and then calculate correlation corrections. The main motivation for developing the method comes from the need to analyse whether different isotope abundances in early universe can contribute to the observed anomalies in quasar absorption spectra. The current best explanation for these anomalies is the assumption that the fine structure constant, alpha, was smaller at early epoch. We test the isotope shift method by comparing the calculated and experimental isotope shift for the alkali and alkali-like atoms Na, MgII, K, CaII and BaII. The agreement is found to be good. We then calculate the isotope shift for some astronomically relevant transitions in SiII and SiIV, MgII, ZnII and GeII.

We follow our previous paper on possible cosmological variation of weak scale (quark masses) and strong scale, inspired by data on cosmological variation of the electromagnetic fine structure constant from distant quasar (QSO) absorption spectra. In this work we identify the {\em strange quark mass} $m_s$ as the most important quantity, and the {\em sigma meson mass} as the ingredient of the nuclear forces most sensitive to it. As a result, we claim significantly stronger limits on ratio of weak/strong scale ($W=m_s/\Lambda_{QCD}$) variation following from our previous discussion of primordial Big-Bang Nucleosynthesis ($|\delta W/W|<0.006$) and Oklo natural nuclear reactor ($|\delta W/W|<1.2 \cdot 10^{-10}$; there is also a non-zero solution $\delta W/W=(-0.56 \pm 0.05) \cdot 10^{-9}$) .

We have previously presented 5.7-sigma evidence for a smaller alpha at redshifts 0.2 < z_abs < 3.7 from a sample of 128 Keck/HIRES quasar absorption systems: da/a = (-0.57 +/- 0.10) x 10^{-5}. A non-zero da/a manifests itself as a distinct pattern of shifts in the measured absorption line wavelengths. The statistical error is now small: we do detect small line shifts in the HIRES data. Whether these shifts are due to systematic errors or due to real variation in alpha is now the central question. Here we summarize the two potentially most important systematic effects: atmospheric dispersion and isotopic abundance evolution. Previously, these have been difficult to quantify/model but here we find that neither of them can explain our results. Furthermore, the HIRES spectra themselves contain no evidence for these effects. Independent measurements of da/a with a different telescope and spectrograph are now crucial if we are to rule out or confirm the present evidence for a variable alpha.

We report preliminary results from a third sample of quasar absorption line spectra from the Keck telescope which has been studied to search for any possible variation of the fine structure constant, alpha. This third sample, which is larger than the sum of the two previously published samples, shows the same effect, and also gives, as do the previous two samples, a significant result. The combined sample yields a highly significant effect, da/a = (alpha_z - alpha_0)/alpha_0 = -0.57 +/- 0.10 x 10^{-5}, averaged over the redshift range 0.2 < z < 3.7. We include a brief discussion of small-scale kinematic structure in quasar absorbing clouds. However, kinematics are unlikely to impact significantly on the averaged non-zero da/a above, and we have so far been unable to identify any systematic effect which can explain it. New measurements of quasar spectra obtained using independent instrumentation and telescopes are required to properly check the Keck results.

We present a short review of the current quasar (QSO) absorption line constraints on possible variation of the fine structure constant, alpha = e^2/(hbar*c). Particular attention is paid to recent optical Keck/HIRES spectra of 49 absorption systems which indicate a smaller alpha in the past (Murphy et al. 2001, Webb et al. 2001). Here we present new preliminary results from 128 absorption systems: da/a = (-0.57 +/- 0.10) x 10^{-5} over the redshift range 0.2 < z < 3.7, in agreement with the previous results. Known potential systematic errors cannot explain these results. We compare them with strong `local' constraints and discuss other (radio and millimetre-wave) QSO absorption line constraints on variations in alpha^2 * g_p and alpha^2 * g_p * m_e/m_p (g_p is the proton g-factor and m_e/m_p is the electron/proton mass ratio). Finally, we discuss future efforts to rule out or confirm the current 5.7 sigma optical detection.

We discuss limits on variation of $(m_q/\Lambda_{QCD})$. The results are obtained by studying $n-\alpha$-interaction during Big Bang, Oklo natural nuclear reactor data and limits on variation of the proton $g$-factor from quasar absorpion spectra.

A recent paper [DDL] claims that the increase of the proton electric charge e leads to a (forbidden) decrease of black hole entropy, therefore, possible evidence for variation of $\alpha=e^2/\hbar c$ [alpha] must be interpreted as a decrease of the speed of light. We argue that purely theoretical consideration of black holes possibly cannot give any model-independent limitations on variation of fundamental constants.

In a recent paper [1] Lineweaver and Davis performed a statistical analysis to claim that the rapidity of biogenesis on Earth indicates high probability of biogenesis on terrestrial- type planets. We argue that the rapid appearance of life on Earth hardly tells us anything about the probability of life to appear on other planet. The conclusion should be different. The rapid initial biogenesis is consistent with a large number ($N \sim 10$) of crucial steps in evolution from simplest life forms to humans.

We present a catalogue of the 322 damped Lyman alpha absorbers taken from the literature. All damped Lyman alpha absorbers are included, with no selection on redshift or quasar magnitude. Of these, 123 are candidates and await confirmation using high resolution spectroscopy. For all 322 objects we catalogue the radio properties of the background quasars, where known. Around 60 quasars have radio flux densities above 0.1 Jy and approximately half of these have optical magnitudes brighter than V = 18. This compilation should prove useful in several areas of extragalactic/cosmological research.

Recent data on cosmological variation of the electromagnetic fine structure constant from distant quasar (QSO) absorption spectra have inspired a more general discussion of possible variation of other constants. We discuss variation of strong scale and quark masses. We derive the limits on their relative change from (i) primordial Big-Bang Nucleosynthesis (BBN); (ii) Oklo natural nuclear reactor, (iii) quasar absorption spectra, and (iv) laboratory measurements of hyperfine intervals.

We performed accurate calculation of $\alpha$-dependence ($\alpha=e^2/hc$) of the transition frequencies for ions, which are used in a search for the variation of the fine structure constant $\alpha$ in space-time. We use Dirac-Hartree-Fock method as a zero approximation and then the many-body perturbation theory and configuration interaction methods to improve the results. An important problem of level pseudocrossing (as functions of $\alpha$) is considered. Near the crossing point the derivative of frequencies over $\alpha$ varies strongly (including change of the sign). This makes it very sensitive to the position of the crossing point. We proposed a semiempirical solution of the problem which allows to obtain accurate results.

An experimental search for variation in the fundamental coupling constants is strongly motivated by modern high-energy physics theories. Comparison of quasar absorption line spectra with laboratory spectra provides a sensitive probe for variability of the fine structure constant, alpha, over cosmological time-scales. We have previously developed and applied a new method providing an order of magnitude gain in precision over previous optical astrophysical constraints. Here we extend that work by including new quasar spectra of damped Lyman-alpha absorption systems. We also re-analyse our previous lower redshift data and confirm our initial results. The constraints on alpha come from simultaneous fitting of absorption lines of subsets of the following species: Mg I, Mg II, Al II, Al III, Si II, Cr II, Fe II, Ni II and Zn II. We present a detailed description of our methods and results based on an analysis of 49 quasar absorption systems (towards 28 QSOs) covering the redshift range 0.5 < z < 3.5. There is statistical evidence for a smaller alpha at earlier epochs: da/a = (-0.72 +/- 0.18) * 10^{-5}. The new and original samples are independent but separately yield consistent and significant non-zero values of da/a. We summarise the results of a thorough investigation of systematic effects published in a companion paper. The value we quote above is the raw value, not corrected for any of these systematic effects. The only significant systematic effects so far identified, if removed from our data, would lead to a more significant deviation of da/a from zero.

We describe the results of a search for time variability of the fine structure constant, alpha, using absorption systems in the spectra of distant quasars. Three large optical datasets and two 21cm/mm absorption systems provide four independent samples, spanning 23% to 87% of the age of the universe. Each sample yields a smaller alpha in the past and the optical sample shows a 4-sigma deviation: da/a = -0.72 +/- 0.18 x 10^{-5} over the redshift range 0.5 < z < 3.5. We find no systematic effects which can explain our results. The only potentially significant systematic effects push da/a towards positive values, i.e. our results would become more significant were we to correct for them.

Comparison of quasar absorption line spectra with laboratory spectra provides a precise probe for variability of the fine structure constant, alpha, over cosmological time-scales. We constrain variation in alpha in 21 Keck/HIRES Si IV absorption systems using the alkali doublet (AD) method in which changes in alpha are related to changes in the doublet spacing. The precision obtained with the AD method has been increased by a factor of 3: da/a = (-0.5 +/- 1.3) * 10^{-5}. We also analyse potential systematic errors in this result. Finally, we compare the AD method with the many-multiplet method which has achieved an order of magnitude greater precision and we discuss the future of the AD method.

Comparison of quasar absorption spectra with laboratory spectra allow us to probe possible variations in the fundamental constants over cosmological time-scales. In a companion paper we present an analysis of Keck/HIRES spectra and report possible evidence suggesting that the fine structure constant, alpha, may have been smaller in the past: da/a = (-0.72 +/- 0.18) * 10^{-5} over the redshift range 0.5 < z < 3.5. In this paper we describe a comprehensive investigation into possible systematic effects. Most of these do not significantly influence our results. When we correct for those which do produce a significant systematic effect in the data, the deviation of da/a from zero becomes more significant. We are lead increasingly to the interpretation that alpha was slightly smaller in the past.

QSO absorption spectra provide an extremely useful probe of possible cosmological variation in various physical constants. Comparison of H I 21cm absorption with corresponding molecular (rotational) absorption spectra allows us to constrain variation in y=alpha^2*g_p where alpha is the fine structure constant and g_p is the proton g-factor. We analyse spectra of two QSOs, PKS 1413+135 and TXS 0218+357, and derive values of dy/y at absorption redshifts of z=0.2467 and 0.6847 by simultaneous fitting of the H I 21cm and molecular lines. We find dy/y=(-0.20 +/- 0.44)*10^{-5} and dy/y=(-0.16 +/- 0.54)*10^{-5} respectively, indicating an insignificantly smaller y in the past. We compare our results with other recent constraints from the same two QSOs (Drinkwater et al. 1998; Carilli et al. 2000) and with our recent optical constraints which indicated a smaller alpha at higher redshifts.

Theories unifying gravity and other interactions suggest the possibility of spatial and temporal variation of physical ``constants'' in the Universe. Detection of high redshift absorption systems intersecting the sight lines towards distant quasars provide a powerful tool for measuring these variations. In the present paper we demonstrate that high sensitivity to variation of the fine structure constant alpha can be obtained by comparing cosmic and laboratory spectra of the Ni II ion. Relativistic effects in Ni II reveal many interesting features. The Ni II spectrum exhibits avoided level crossing phenomenon under variation of alpha and the intervals between the levels have strong nonlinear dependencies on relativistic corrections. The values of the transition frequency shifts, due to the change of alpha, vary significantly from state to state including change of the sign. This enhances the sensitivity to the variation of alpha and reduces possible systematic errors. The calculations of alpha-dependence of the nickel ion spectral lines that are detectable in quasar absorption spectra have been performed using a relativistic configuration interaction method.

It is believed that there may have been a large number of black holes formed in the very early universe. These would have quantised masses. A charged ``elementary black hole'' (with the minimum possible mass) can capture electrons, protons and other charged particles to form a ``black hole atom''. We find the spectrum of such an object with a view to laboratory and astronomical observation of them, and estimate the lifetime of the bound states. There is no limit to the charge of the black hole, which gives us the possibility of observing Z>137 bound states and transitions at the lower continuum. Negatively charged black holes can capture protons. For Z>1, the orbiting protons will coalesce to form a nucleus (after beta-decay of some protons to neutrons), with a stability curve different to that of free nuclei. In this system there is also the distinct possibility of single quark capture. This leads to the formation of a coloured black hole that plays the role of an extremely heavy quark interacting strongly with the other two quarks. Finally we consider atoms formed with much larger black holes.

A method offering an order of magnitude sensitivity gain is described for using quasar spectra to investigate possible time or space variation in the fine structure constant, alpha. Applying the technique to a sample of 30 absorption systems, spanning redshifts 0.5 < z< 1.6, obtained with the Keck I telescope, we derive limits on variations in alpha over a wide range of epochs. For the whole sample Delta(alpha)/alpha = -1.1 +/- 0.4 x 10^{-5}. This deviation is dominated by measurements at z > 1, where Delta(alpha)/alpha = -1.9 +/- 0.5 x 10^{-5}. For z < 1, Delta(alpha)/alpha = -0.2 +/- 0.4 x 10^{-5}, consistent with other known constraints. Whilst these results are consistent with a time-varying alpha, further work is required to explore possible systematic errors in the data, although careful searches have so far not revealed any.

Theories unifying gravity and other interactions suggest the possibility of spatial and temporal variation of physical ``constants'' in the Universe. Detection of high-redshift absorption systems intersecting the sight lines towards distant quasars provide a powerful tool for measuring these variations. We have previously demonstrated that high sensitivity to the variation of the fine structure constant $\alpha$ can be obtained by comparing spectra of heavy and light atoms (or molecules). Here we describe new calculations for a range of atoms and ions, most of which are commonly detected in quasar spectra: Fe II, Mg II, Mg I, C II, C IV, N~V, O I, Al III, Si II, Si IV, Ca I, Ca II, Cr II, Mn II, Zn II, Ge II (see the results in Table 3). The combination of Fe II and Mg II, for which accurate laboratory frequencies exist, have already been used to constrain $\alpha$ variations. To use other atoms and ions, accurate laboratory values of frequencies of the strong E1-transitions from the ground states are required. We wish to draw the attention of atomic experimentalists to this important problem. We also discuss a mechanism which can lead to a greatly enhanced sensitivity for placing constraints on variation on fundamental constants. Calculations have been performed for Hg II, Yb II, Ca I and Sr II where there are optical transitions with the very small natural widths, and for hyperfine transition in Cs I and Hg II.

Detection of high-redshift absorption in the optical spectra of quasars have provided a powerful tool to measure spatial and temporal variations of physical ``constants'' in the Universe. It is demonstrated that high sensitivity to the variation of the fine structure constant alpha can be obtained from a comparison of the spectra of heavy and light atoms (or molecules). We have performed calculations for the pair FeII and MgII for which accurate quasar and laboratory spectra are available. A possibility of $10^5$ times enhanced effects of the fundamental constants variation suitable for laboratory measurements is also discussed.

Recent detections of high-redshift absorption by both atomic hydrogen and molecular gas in the radio spectra of quasars have provided a powerful tool to measure possible temporal and spatial variations of physical `constants' in the Universe. We have compared the frequency of high-redshift Hydrogen 21 cm absorption with that of associated molecular absorption in two quasars to place new (1 sigma) upper limits on any variation in y = g_p.alpha^2 (where alpha is the fine structure constant and g_p is the proton g-factor) of |Delta(y)/y| < 5.10^{-6} at redshifts z=0.25 and z=0.68. These quasars are separated by a comoving distance of 3000 Mpc (H_0=75 km/s/Mpc and q_0=0). We also derive limits on the time rates of change of |(d.g_p/d.t)/g_p| < 10^{-15}/y and |(d.alpha/d.t)/\alpha| < 5.10^{-16}/y between the present epoch and z=0.68. These limits are more than an order of magnitude smaller than previous results derived from high-redshift measurements.

We have compared the frequency of high-redshift Hydrogen HI 21 cm absorption with that of associated molecular absorption in two quasars to place new (1 sigma) upper limits on any variation in y = g_p.alpha^2 (where alpha is the fine structure constant and g_p is the proton g-factor) of |Delta(y)/y| < 5.10^{-6} at redshifts z=0.25 and z=0.68.