Search

Motivated by the recently discovered high-$T_c$ bilayer nickelate superconductor La$_3$Ni$_2$O$_7$, we comprehensively research a bilayer $2\times2\times2$ cluster for different electronic densities $n$ by using the Lanczos method. We also employ the random-phase approximation to quantify the first magnetic instability with increasing Hubbard coupling strength, also varying $n$. Based on the spin structure factor $S(q)$, we have obtained a rich magnetic phase diagram in the plane defined by $n$ and $U/W$, at fixed Hund coupling. We have observed numerous states, such as A-AFM, Stripes, G-AFM, and C-AFM. For half-filling $n=2$ (two electrons per Ni site, corresponding to $N$ = 16 electrons), the canonical superexchange interaction leads to a robust G-AFM state $(\pi,\pi,\pi)$ with antiferromagnetic couplings in plane and between layers. By increasing or decreasing electronic densities, ferromagnetic tendencies emerge from the ``half-empty'' and ``half-full'' mechanisms, leading to many other interesting magnetic tendencies. In addition, the spin-spin correlations become weaker both in the hole or electron doping regions compared with half-filling. At $n = 1.5$ (or $N=12$), density corresponding to La$_3$Ni$_2$O$_7$, we obtained the ``Stripe 2'' ground state (antiferromagnetic coupling in one in-plane direction, ferromagnetic coupling in the other, and antiferromagnetic coupling along the $z$-axis) in the $2\times2\times2$ cluster. In addition, we obtained a much stronger AFM coupling along the $z$-axis than the magnetic coupling in the $xy$ plane. The random-phase approximation calculations with varying $n$ give very similar results as Lanczos. Meanwhile, a state with $q/\pi = (0.6, 0.6, 1)$ close to the E-phase wavevector is found in our RPA calculations by slightly reducing the filling to $n=1.25$, possibly responsible for the E-phase SDW recently observed in experiments., Comment: 12 pages, 9 figures

Quantum computers hold the promise of solving certain problems that lie beyond the reach of conventional computers. Establishing this capability, especially for impactful and meaningful problems, remains a central challenge. One such problem is the simulation of nonequilibrium dynamics of a magnetic spin system quenched through a quantum phase transition. State-of-the-art classical simulations demand resources that grow exponentially with system size. Here we show that superconducting quantum annealing processors can rapidly generate samples in close agreement with solutions of the Schr\"odinger equation. We demonstrate area-law scaling of entanglement in the model quench in two-, three- and infinite-dimensional spin glasses, supporting the observed stretched-exponential scaling of effort for classical approaches. We assess approximate methods based on tensor networks and neural networks and conclude that no known approach can achieve the same accuracy as the quantum annealer within a reasonable timeframe. Thus quantum annealers can answer questions of practical importance that classical computers cannot.

Quantum Error Mitigation (QEM) presents a promising near-term approach to reduce error when estimating expectation values in quantum computing. Here, we introduce QEM techniques tailored for quantum annealing, using Zero-Noise Extrapolation (ZNE). We implement ZNE through zero-temperature extrapolation as well as energy-time rescaling. We conduct experimental investigations into the quantum critical dynamics of a transverse-field Ising spin chain, demonstrating the successful mitigation of thermal noise through both of these techniques. Moreover, we show that energy-time rescaling effectively mitigates control errors in the coherent regime where the effect of thermal noise is minimal. Our ZNE results agree with exact calculations of the coherent evolution over a range of annealing times that exceeds the coherent annealing range by almost an order of magnitude., Comment: 10 pages, 5 figures

Motivated by the recent experimental realization of minimal Kitaev chains using quantum dots, we investigate the Majorana zero modes (MZM) in $Y$-shape Kitaev wires. We solve the associated Kitaev models analytically at the sweet spot ($t_h=\Delta$) and derive the exact form of MZM wave-functions in this geometry. The novelty of our result is the observation of multi-site MZMs located near the junction center on the nearby edge sites of each leg. This result is important for potential braiding of Majoranas and the performance of $Y$-junctions made from arrays of quantum dots. Furthermore, we study the stability of local (single-site) and multi-site MZMs modes in the presence of Coulomb repulsion, using density matrix renormalization group theory. Our local density-of-states calculation shows that these multi-site MZMs are as equally topologically protected as the single-site MZMs when in the presence of Coulomb repulsion or when away from the sweet-spot., Comment: 8 pages and 5 figures

Ab-initio and continuum model studies predicted that the $\Gamma$ valley transition metal dichalcogenide (TMD) homobilayers could simulate the conventional multi-orbital Hubbard model on the moir\'e Honeycomb lattice. Here, we perform the Wannierization starting from the continuum model and show that a more general Kanamori-Moir\'e-Hubbard model emerges, beyond the extensively studied standard multi-orbital Hubbard model, which can be used to investigate the many-body physics in the $\Gamma$ valley TMD homobilayers. Using the unrestricted Hartree-Fock and Lanczos techniques, we study these half-filled multi-orbital moir\'e bands. By constructing the phase diagrams we predict the presence of an antiferromagnetic state and in addition we found unexepected and dominant states, such as a S = 1 ferromagnetic insulator and a charge density wave state. Our theoretical predictions made using this model can be tested in future experiments on the $\Gamma$ valley TMD homobilayers.

Lattice Wigner crystal states stabilized by long-range Coulomb interactions have recently been realized in two-dimensional moir\'e materials. We employ large-scale unrestricted Hartree-Fock techniques to unveil the magnetic phase diagrams of honeycomb lattice Wigner crystals. For the three lattice filling factors with the largest charge gaps, n = 2/3, 1/2, 1/3, the magnetic phase diagrams contain multiple phases, including ones with non-collinear and non-coplanar spin arrangements. We discuss magnetization evolution with external magnetic field, which has potential as an experimental signature of exotic spin states. Our theoretical results could potentially be validated in moir\'e materials formed from group VI transition metal dichalcogenide twisted homobilayers.

We studied a multi-orbital Hubbard model at half-filling for two and three orbitals per site on a two-site cluster via full exact diagonalization, in a wide range for the onsite repulsion $U$, from weak to strong coupling, and multiple ratios of the Hund coupling $J_H$ to $U$. The hopping matrix elements among the orbitals were also varied extensively. At intermediate and large $U$, we mapped the results into a Heisenberg model. For two orbitals per site, the mapping is into a $S=1$ Heisenberg model where by symmetry both nearest-neighbor $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})$ and $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{2}$ are allowed, with respective couplings $J_1$ and $J_2$. For the case of three orbitals per site, the mappping is into a $S=3/2$ Heisenberg model with $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})$, $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{2}$, and $(\mathbf{S}_{i}\cdot\mathbf{S}_{j})^{3}$ terms, and respective couplings $J_1$, $J_2$, and $J_3$. The strength of these coupling constants in the Heisenberg models depend on the $U$, $J_H$, and hopping amplitudes of the underlying Hubbard model. Our study allows to establish bounds on how large the ratios $J_2/J_1$ and $J_3/J_1$ can be. We show that those ratios are severely limited and, as a crude guidance, we conclude that $J_2/J_1$ is less than 0.4 and $J_3/J_1$ is less than 0.2, establishing bounds on effective models for strongly correlated Hubbard systems., Comment: 8 pages, 7 figures

The competition between spin-orbit coupling $\lambda$ and electron-electron interaction $U$ leads to a plethora of novel states of matter, extensively studied in the context of $t_{2g}^4$ and $t_{2g}^5$ materials, such as ruthenates and iridates. Excitonic magnets -- the antiferromagnetic state of bounded electron-hole pairs -- is a prominent example of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground-state of spin-orbit-coupled $t_{2g}^4$ Hubbard models. Here, we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the $\lambda=0$ state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wavevector $q\to 0$. In the excitonic condensation regime at large $U$, we also have found a novel high-energy $\pi$-mode composed solely of orbital excitations. These unique fingerprints of the $t_{2g}^4$ excitonic magnet are important in the analysis of neutron and RIXS experiments.

We report the results of a Hartree-Fock study applied to interacting electrons moving in two different bipartite lattices: the dice and the Lieb lattices, at half-filling. Both lattices develop ferrimagnetic order in the phase diagram $U$-$\lambda$, where $U$ is the Hubbard onsite repulsion and $\lambda$ the Rashba spin-orbit coupling strength. Our main result is the observation of an unexpected multitude of topological phases for both lattices. All these phases are ferrimagnetic, but they differ among themselves in their set of six Chern numbers (six numbers because the unit cells have three atoms). The Chern numbers $|C|$ observed in our study range from 0 to 3, showing that large Chern numbers can be obtained by the effect of electronic correlations, adding to the recently discussed methodologies to increase $|C|$ based on extending the hopping range in tight-binding models, using sudden quenches, or photonic crystals, all without including electronic interactions., Comment: 12 pages, 12 figures

Recent DFT calculations for Ba2CoO4 (BCO) and neutron scattering experiments for SrRuO3 (SRO) have shown that oxygen develops a magnetic polarization. Moreover, DFT calculations for these compounds also unveiled unexpected nodes in the spin density, both along Co-O and Ru-O. For BCO, the overall antiferromagnetic state in its triangular lattice contains unusual zigzag spin patterns. Here, using simple model calculations supplemented by DFT we explain and extend these results. We predict that ligands that in principle should be spinless, such as O$^{2-}$, will develop a net polarization when they act as electronic bridges between transition metal (TM) spins ferromagnetically ordered, regardless of the number of intermediate ligand atoms. The reason is the hybridization between atoms and mobility of the electrons with spins opposite to those of the closest TM atoms. Moreover, for bonds with TMs antiferromagnetically ordered, counterintuitively our calculations show that oxygens should also have a net magnetization for the super-super-exchange cases TM-O-O-TM while for only one oxygen, as in Cu-O-Cu, the O-polarization should cancel. Our simple model also allows us to explain the presence of nodes based on the antibonding character of the dominant singly occupied molecular orbitals along the TM-O bonds. Finally, the zigzag pattern order becomes the ground state mainly due to the influence of the Hubbard $U$, that creates the moments, in combination with a robust easy-axis anisotropy that suppresses the competing 120$^{\circ}$ degree antiferromagnetic order of a triangular lattice. Our predictions are generic and should be applicable to any other compound with characteristics similar to those of BCO and SRO., Comment: 15 pages, 17 figures

In strongly correlated electronic systems, several novel physical properties are induced by the orbital degree of freedom. In particular, orbital degeneracy near the Fermi level leads to spontaneous symmetry breaking, such as the nematic state in FeSe and the orbital ordering in several perovskite systems. Here, the novel layered perovskite material CsVF$_4$, with a $3d^2$ electronic configuration, was systematically studied using density functional theory and a multiorbital Hubbard model within the Hatree-Fock approximation. Our results show that CsVF$_4$ should be magnetic, with a G-type antiferromagnetic arrangement in the $ab$ plane and weak antiferromagnetic exchange along the $c$-axis, in agreement with experimental results. Driven by the Jahn-Teller distortion in the VF$_6$ octahedra that shorten the $c$-axis, the system displays an interesting electron occupancy $d_{xy}^1(d_{xz}d_{yz})^1$ corresponding to the lower nondegenerate $d_{xy}$ orbital being half-filled and the other two degenerate $d_{yz}$ and $d_{xz}$ orbitals sharing one electron per site. We show that this degeneracy is broken and a novel $d_{yz}$/$d_{xz}$ staggered orbital pattern is here predicted by both the first-principles and Hubbard model calculations. This orbital ordering is driven by the electronic instability associated with degeneracy removal to lower the energy., Comment: 10 pages, 12 figures

Motivated by recent experimental progress in transition metal oxides with the K$_2$NiF$_4$ structure, we investigate the magnetic and orbital ordering in $\alpha$-Sr$_2$CrO$_4$. Using first principles calculations, first we derive a three-orbital Hubbard model, which reproduces the {\it ab initio} band structure near the Fermi level. The unique reverse splitting of $t_{2g}$ orbitals in $\alpha$-Sr$_2$CrO$_4$, with the $3d^2$ electronic configuration for the Cr$^{4+}$ oxidation state, opens up the possibility of orbital ordering in this material. Using real-space Hartree-Fock for multi-orbital systems, we constructed the ground-state phase diagram for the two-dimensional compound $\alpha$-Sr$_2$CrO$_4$. We found stable ferromagnetic, antiferromagnetic, antiferro-orbital, and staggered orbital stripe ordering in robust regions of the phase diagram. Furthermore, using the density matrix renormalization group method for two-leg ladders with the realistic hopping parameters of $\alpha$-Sr$_2$CrO$_4$, we explore magnetic and orbital ordering for experimentally relevant interaction parameters. Again, we find a clear signature of antiferromagnetic spin ordering along with antiferro-orbital ordering at moderate to large Hubbard interaction strength. We also explore the orbital-resolved density of states with Lanczos, predicting insulating behavior for the compound $\alpha$-Sr$_2$CrO$_4$, in agreement with experiments. Finally, an intuitive understanding of the results is provided based on a hierarchy between orbitals, with $d_{xy}$ driving the spin order, while electronic repulsion and the effective one dimensionality of the movement within the $d_{xz}$ and $d_{yz}$ orbitals driving the orbital order., Comment: 13 pages, 12 figures

The magnetic and electronic phase diagram of a model for the quasi-one-dimensional alkali metal iron selenide compound Na$_2$FeSe$_2$ is presented. The novelty of this material is that the valence of iron is Fe$^{2+}$ contrary to most other iron-chain compounds with valence Fe$^{3+}$. Using first-principles techniques, we developed a three-orbital tight-binding model that reproduces the {\it ab initio} band structure near the Fermi level. Including Hubbard and Hund couplings and studying the model via the density matrix renormalization group and Lanczos methods, we constructed the ground state phase diagram. A robust region where the block state $\uparrow \uparrow \downarrow \downarrow \uparrow \uparrow \downarrow \downarrow$ is stabilized was unveiled. The analog state in iron ladders, employing 2$\times$2 ferromagnetic blocks, is by now well-established, but in chains a block magnetic order has not been observed yet in real materials. The phase diagram also contains a large region of canonical staggered spin order $\uparrow \downarrow \uparrow \downarrow \uparrow \downarrow \uparrow$ at very large Hubbard repulsion. At the block to staggered transition region, a novel phase is stabilized with a mixture of both states: an inhomogeneous orbital-selective charge density wave with the exotic spin configuration $\uparrow \uparrow \downarrow \uparrow \downarrow \downarrow \uparrow \downarrow$. Our predictions for Na$_2$FeSe$_2$ may guide crystal growers and neutron scattering experimentalists towards the realization of block states in one-dimensional iron-selenide chain materials., Comment: 10 pages, 5 figures

The condensation of spin-orbit-induced excitons in $(t_{2g})^4$ electronic systems is attracting considerable attention. In the large Hubbard U limit, antiferromagnetism was proposed to emerge from the Bose-Einstein Condensation (BEC) of triplons ($J_{\textrm{eff}} = 1$). In this publication, we show that even for the weak and intermediate U regimes, the spin-orbit exciton condensation is possible leading also to staggered magnetic order. The canonical electron-hole excitations (excitons) transform into local triplon excitations at large U , and this BEC strong coupling regime is smoothly connected to the intermediate U excitonic insulator region. We solved the degenerate three-orbital Hubbard model with spin-orbit coupling ($\lambda$) in one-dimensional geometry using the Density Matrix Renormalization Group, while in two-dimensional square clusters we use the Hartree-Fock approximation (HFA). Employing these techniques, we provide the full $\lambda$ vs U phase diagrams for both one- and two- dimensional lattices. Our main result is that at the intermediate Hubbard U region of our focus, increasing $\lambda$ at fixed U the system transitions from an incommensurate spin-density-wave metal to a Bardeen-Cooper-Schrieffer (BCS) excitonic insulator, with coherence length r coh of O(a) and O(10a) in 1d and 2d, respectively, with a the lattice spacing. Further increasing $\lambda$, the system eventually crosses over to the BEC limit (with r coh << a).

We introduce a variational state for one-dimensional two-orbital Hubbard models that intuitively explains the recent computational discovery of pairing in these systems when hole doped. Our Ansatz is an optimized linear superposition of Affleck-Kennedy-Lieb-Tasaki valence bond states, rendering the combination a valence bond liquid dubbed Orbital Resonant Valence Bond. We show that the undoped (one electron/orbital) quantum state of two sites coupled into a global spin singlet is exactly written employing only spin-1/2 singlets linking orbitals at nearest-neighbor sites. Generalizing to longer chains defines our variational state visualized geometrically expressing our chain as a two-leg ladder, with one orbital per leg. As in Anderson's resonating valence-bond state, our undoped variational state contains preformed singlet pairs that via doping become mobile leading to superconductivity. Doped real materials with one-dimensional substructures, two near-degenerate orbitals, and intermediate Hubbard U/W strengths -- W the carrier's bandwidth -- could realize spin-singlet pairing if on-site anisotropies are small. If these anisotropies are robust, spin-triplet pairing emerges.

We study ribbons of the dice two-dimensional lattice (that we call ``dice ladders'') known to have nontrivial topological properties, such as Chern numbers 2 [Wang and Y. Ran, Phys. Rev. B {\bf 84}, 241103 (2011)]. Our main results are two folded: (1) Analyzing the tight-binding model in the presence of Rashba spin-orbit coupling and an external magnetic field, we observed that dice ladders qualitatively display properties similar to their two-dimensional counterpart all the way to the limit of only two legs in the short direction. This includes flat bands near the Fermi level, edge currents and edge charge localization near zero energy when open boundary conditions are used, two chiral edge modes, and a nonzero Hall conductance. (2) We studied the effect of Hubbard correlation $U$ in the two-leg dice ladder using Lanczos and density matrix renormalization group techniques. We show that increasing $U$ the flat bands split without the need of introducing external fields. Moreover, robust ferrimagnetic order develops. Overall, our work establishes dice ladders as a promising playground to study the combined effect of topology and correlation effects, one of the frontiers in Quantum Materials.

Theoretical studies recently predicted the condensation of spin-orbit excitons at momentum $q$=$\pi$ in $t_{2g}^4$ spin-orbit coupled three-orbital Hubbard models at electronic density $n=4$. In parallel, experiments involving iridates with non-integer valence states for the Ir ions are starting to attract considerable attention. In this publication, using the density matrix renormalization group technique we present evidence for the existence of a novel excitonic condensate at $n=3.5$ in a one-dimensional Hubbard model with a degenerate $t_{2g}$ sector, when in the presence of spin-orbit coupling. At intermediate Hubbard $U$ and spin-orbit $\lambda$ couplings, we found an excitonic condensate at the unexpected momentum $q$=$\pi/2$ involving $j_{\textrm{eff}}=3/2,m=\pm1/2$ and $j_{\textrm{eff}}=1/2,m=\pm1/2$ bands in the triplet channel, coexisting with an also unexpected block magnetic order. We also present the entire $\lambda$ vs $U$ phase diagram, at a fixed and robust Hund coupling. Interestingly, this new `block excitonic phase' is present even at large values of $\lambda$, unlike the $n=4$ excitonic phase discussed before. Our computational study helps to understand and predict the possible magnetic phases of materials with $d^{3.5}$ valence and robust spin-orbit coupling.

Motivated by recent proposals of correlation induced insensitivity of d-wave superconductors to impurities, we develop a simple pairing theory for these systems for up to a moderate strength of disorder. Our description implements the key ideas of Anderson, originally proposed for disordered s-wave superconductors, but in addition takes care of the inherent strong electronic repulsion in these compounds, as well as disorder induced inhomogeneities. We first obtain the self-consistent one-particle states, that capture the effects of disorder exactly, and strong correlations using Gutzwiller approximation. These `normal states', representing the interplay of strong correlations and disorder, when coupled through pairing attractions following the path of Bardeen-Cooper-Schrieffer (BCS), produce results nearly identical to those from a more sophisticated Gutzwiller augmented Bogoliubov-de Gennes analysis., Comment: 10 pages, 7 figures

Using the Density Matrix Renormalization Group technique we study the effect of spin-orbit coupling on a three-orbital Hubbard model in the $(t_{2g})^{4}$ sector and in one dimension. Fixing the Hund coupling to a robust value compatible with some multiorbital materials, we present the phase diagram varying the Hubbard $U$ and spin-orbit coupling $\lambda$, at zero temperature. Our results are shown to be qualitatively similar to those recently reported using the Dynamical Mean Field Theory in higher dimensions, providing a robust basis to approximate many-body techniques. Among many results, we observe an interesting transition from an orbital-selective Mott phase to an excitonic insulator with increasing $\lambda$ at intermediate $U$. In the strong $U$ coupling limit, we find a non-magnetic insulator with an effective angular momentum $\langle(\textbf{J}^{eff})^{2}\rangle \ne 0$ near the excitonic phase, smoothly connected to the $\langle(\textbf{J}^{eff})^{2}\rangle = 0$ regime. We also provide a list of quasi-one dimensional materials where the physics discussed in this publication could be realized.

We study the effect of spin-orbit interaction on one-dimensional U(1)-invariant frustrated magnets with dominant critical nematic fluctuations. The spin-orbit coupling explicitly breaks the U(1) symmetry of arbitrary global spin rotations about the high-symmetry axis down to $Z_2$ (invariance under a $\pi$-rotation). Given that the nematic order parameter is invariant under a $\pi$-rotation, it is relevant to ask if other discrete symmetries can be spontaneously broken. Here we demonstrate that the spin-orbit coupling induces a bond density wave that spontaneously breaks the translational symmetry and opens a gap in the excitation spectrum.

We employ a recently developed computational many-body technique to study for the first time the half-filled Anderson-Hubbard model at finite temperature and arbitrary correlation ($U$) and disorder ($V$) strengths. Interestingly, the narrow zero temperature metallic range induced by disorder from the Mott insulator expands with increasing temperature in a manner resembling a quantum critical point. Our study of the resistivity temperature scaling $T^{\alpha}$ for this metal reveals non Fermi liquid characteristics. Moreover, a continuous dependence of $\alpha$ on $U$ and $V$ from linear to nearly quadratic was observed. We argue that these exotic results arise from a systematic change with $U$ and $V$ of the "effective" disorder, a combination of quenched disorder and intrinsic localized spins.

A recently introduced one-dimensional three-orbital Hubbard model displays orbital-selective Mott phases with exotic spin arrangements such as spin block states [J. Rinc\'on {\em et al.}, Phys. Rev. Lett. \textbf{112}, 106405 (2014)]. In this publication we show that the Constrained-Path Quantum Monte Carlo (CPQMC) technique can accurately reproduce the phase diagram of this multiorbital one-dimensional model, paving the way to future CPQMC studies in systems with more challenging geometries, such as ladders and planes. The success of this approach relies on using the Hartree-Fock technique to prepare the trial states needed in CPQMC. We also study a simplified version of the model where the pair-hopping term is neglected and the Hund coupling is restricted to its Ising component. The corresponding phase diagrams are shown to be only mildly affected by the absence of these technically difficult-to-implement terms. This is confirmed by additional Density Matrix Renormalization Group and Determinant Quantum Monte Carlo calculations carried out for the same simplified model, with the latter displaying only mild Fermion sign problems. We conclude that these methods are able to capture quantitatively the rich physics of the several orbital-selective Mott phases (OSMP) displayed by this model, thus enabling computational studies of the OSMP regime in higher dimensions, beyond static or dynamic mean field approximations., Comment: 13 pages, 14 figures

Temporal variation of secondary cosmic rays (SCR) flux was measured during the total lunar eclipse on December 10, 2011 and the subsequent full moon on January 8, 2012. The measurements were done at Department of Physics, University of Mumbai, Mumbai (Geomagnetic latitude: 10.6 N), India using NaI (Tl) scintillation detector by keeping energy threshold of 200 KeV. The SCR flux showed approximately 8.1% enhancement during the lunar eclipse as compared to the average of pre- and post-eclipse periods. Weather parameters (temperature and relative humidity) were continuously monitored and their correlations with temporal variation in SCR flux were examined. The influences of geomagnetic field, interplanetary parameters and tidal effect on SCR flux were considered. Qualitative analysis of SCR flux variation indicates that the known factors affecting SCR flux fail to explain observed enhancement during the eclipse. This enhancement during lunar eclipse and widely reported decrease during solar eclipses may unravel hitherto unnoticed factors modulating SCR flux., Comment: 12 pages, 8 figures