Optically levitated particles have great potential to form the basis of novel quantum- enhanced sensors. These systems are very well suited for inertial sensing, as the particles are isolated from the environment when they are levitated at low pressures. However, there are many challenges in the experimental realization that may affect the performance of these systems. For example, optical aberrations in the wavefront of the trapping laser which arise from optical elements or misalignment have a great impact on the trapping potential. The detrimental effect of optical aberrations has not been thoroughly studied, and usually they are iteratively corrected, giving some conflicting results depending on the figures of merit that are used. In this work, we present a thorough study of the effects of structuring the wavefront of the trapping beams. We observe that clean beams, i.e. highly focused beams with unaberrated wavefronts, may be used to optimize the longitudinal frequencies, at the cost of the transversal ones. Our work is based in a combination of experimental studies using a complete basis of orthogonal polynomials (Zernike polynomials) to control the wavefront and a set of numerical calculations, which allow us to compare the impact of structured wavefronts on the quality of traps for optically levitated particles in vacuum. This will have direct applications in quantum sensing and fundamental studies of quantum mechanics, as it allows the reduction of optical backaction and thermal decoherence of the particles.

The calculation of quantum-geometric properties of Bloch electrons -- Berry curvature, quantum metric, orbital magnetic moment and effective mass -- was implemented in a pseudopotential plane-wave code. The starting point was the first derivative of the periodic part of the wavefunction psi_k with respect to the wavevector k. This was evaluated with perturbation theory by solving a Sternheimer equation. Comparison of effective masses obtained from perturbation theory for silicon and gallium arsenide with carefully-converged numerical second derivatives of band energies confirms the high precision of the method. Calculations of quantum-geometric quantities for gapped graphene were performed by adding a bespoke symmetry-breaking potential to first-principles graphene. As the two bands near the opened gap are reasonably isolated, the results could be compared with those obtained from an analytical two-band model, allowing to assess the strengths and limitations of such widely-used models. The final application was trigonal tellurium, where quantum-geometric quantities flip sign with chirality.

We present a numerical implementation, based on Wannier interpolation, of a Kubo-Greenwood formalism for computing the spatially dispersive optical conductivity in crystals at first order in the wave vector of light. This approach is more efficient than direct $\textit{ab initio}$ methods because, with less computational cost, it allows for a much finer sampling of reciprocal space, resulting in better resolved spectra. Moreover, Wannier interpolation avoids errors arising from truncation of the sums over conduction bands when evaluating the spatially dispersive optical matrix elements. We validate our method by computing the optical activity spectrum of selected crystals, both polar (GaN) and chiral (trigonal Te, trigonal Se, and $\alpha$-quartz), and comparing with existing literature.

We present a systematic study of bound relations between different electronic properties of magnetic crystals: electron density, effective mass, orbital magnetization, localization length, Chern invariant, and electric susceptibility. All relations are satisfied for a group of low-lying bands, while some remain valid for upper bands. New results include a lower bound on the electric susceptibility of Chern insulators, and an upper bound on the sum-rule part of the orbital magnetization. In addition, bounds involving the Chern invariant are generalized from two dimensions (Chern number) to three (Chern vector). Bound relations are established for metals as well as insulators, and are illustrated for model systems. The manner in which they approach saturation in a model Chern insulator with tunable flat bands is analyzed in terms of the optical absorption spectrum.

We present a systematic study of bound relations between different electronic properties of magnetic crystals: electron density, effective mass, orbital magnetization, localization length, Chern invariant, and electric susceptibility. All relations are satisfied for a group of low-lying bands, while some remain valid for upper bands. New results include a lower bound on the electric susceptibility of Chern insulators, and an upper bound on the sum-rule part of the orbital magnetization. In addition, bounds involving the Chern invariant are generalized from two dimensions (Chern number) to three (Chern vector). Bound relations are established for metals as well as insulators, and are illustrated for model systems. The manner in which they approach saturation in a model Chern insulator with tunable flat bands is analyzed in terms of the optical absorption spectrum.

We discuss a Kohn-Luttinger-like mechanism for superconductivity in Bernal bilayer graphene and rhombohedral trilayer graphene. Working within the continuum model description, we find that the screened long-range Coulomb interaction alone gives rise to superconductivity with critical temperatures that agree with experiments. We observe that the order parameter changes sign between valleys, which implies that both materials are valley-singlet, spin-triplet superconductors. Adding Ising spin-orbit coupling leads to a significant enhancement in the critical temperature, also in line with experiment, and the superconducting order parameter shows locking between the spin and valley degrees of freedom.

The current density ${\bf j}^{\rm{\bf B}}$ induced in a clean metal by a slowly-varying magnetic field ${\bf B}$ is formulated as the low-frequency limit of natural optical activity, or natural gyrotropy. Working with a multiband Pauli Hamiltonian, we obtain from the Kubo formula a simple expression for $\alpha^{\rm gme}_{ij}=j^{\rm{\bf B}}_i/B_j$ in terms of the intrinsic magnetic moment (orbital plus spin) of the Bloch electrons on the Fermi surface. An alternate semiclassical derivation provides an intuitive picture of the effect, and takes into account the influence of scattering processes in dirty metals. This "gyrotropic magnetic effect" is fundamentally different from the chiral magnetic effect driven by the chiral anomaly and governed by the Berry curvature on the Fermi surface, and the two effects are compared for a minimal model of a Weyl semimetal. Like the Berry curvature, the intrinsic magnetic moment should be regarded as a basic ingredient in the Fermi-liquid description of transport in broken-symmetry metals.

Once only accessible in nearby galaxies, we can now study individual stars across much of the observable universe aided by galaxy-cluster gravitational lenses. When a star, compact object, or multiple such objects in the foreground galaxy-cluster lens become aligned, they can magnify a background individual star, and the timescale of a magnification peak can limit its size to tens of AU. The number and frequency of microlensing events therefore opens a window into the population of stars and compact objects, as well as high-redshift stars. To assemble the first statistical sample of stars in order to constrain the initial mass function (IMF) of massive stars at redshift z=0.7-1.5, the abundance of primordial black holes in galaxy-cluster dark matter, and the IMF of the stars making up the intracluster light, we are carrying out a 192-orbit program with the Hubble Space Telescope called "Flashlights," which is now two-thirds complete owing to scheduling challenges. We use the ultrawide F200LP and F350LP long-pass WFC3 UVIS filters and conduct two 16-orbit visits separated by one year. Having an identical roll angle during both visits, while difficult to schedule, yields extremely clean subtraction. Here we report the discovery of more than a dozen bright microlensing events, including multiple examples in the famous "Dragon Arc" discovered in the 1980s, as well as the "Spocks" and "Warhol" arcs that have hosted already known supergiants. The ultradeep observer-frame ultraviolet-through-optical imaging is sensitive to hot stars, which will complement deep James Webb Space Telescope infrared imaging. We are also acquiring Large Binocular Telescope LUCI and Keck-I MOSFIRE near-infrared spectra of the highly magnified arcs to constrain their recent star-formation histories.

Over the last two decades, following the early developments on maximally localized Wannier functions, an ecosystem of electronic-structure simulation techniques and software packages leveraging the Wannier representation has flourished. This environment includes codes to obtain Wannier functions and interfaces with first-principles simulation software, as well as an increasing number of related post-processing packages. Wannier functions can be obtained for isolated or extended systems (both crystalline and disordered), and can be used to understand chemical bonding, to characterize electric polarization, magnetization, and topology, or as an optimal basis set, providing very accurate interpolations in reciprocal space or large-scale Hamiltonians in real space. In this review, we summarize the current landscape of techniques, materials properties and simulation codes based on Wannier functions that have been made accessible to the research community, and that are now well integrated into what we term a \emph{Wannier function software ecosystem}. First, we introduce the theory and practicalities of Wannier functions, starting from their broad domains of applicability to advanced minimization methods using alternative approaches beyond maximal localization. Then we define the concept of a Wannier ecosystem and its interactions and interoperability with many quantum simulations engines and post-processing packages. We focus on some of the key properties and capabilities that are empowered by such ecosystem\textemdash from band interpolations and large-scale simulations to electronic transport, Berryology, topology, electron-phonon couplings, dynamical mean-field theory, embedding, and Koopmans functionals\textemdash concluding with the current status of interoperability and automation. [...]

Intercellular exchange networks are essential for the adaptive capabilities of populations of cells. While diffusional exchanges have traditionally been difficult to map, recent advances in nanotechnology enable precise probing of exchange fluxes with the medium at single-cell resolution. Here we introduce a tiling-based method to reconstruct the dynamic unfolding of exchange networks from flux data, subsequently applying it to an experimental mammalian co-culture system where lactate exchanges affect the acidification of the environment. We observe that the network, which initially exhibits a dense matrix of exchanges, progressively breaks up into small disconnected clusters of cells. To explain this behaviour, we develop a two-parameter Maximum-Entropy multicellular metabolic model that incorporates diffusion-driven exchanges through a set of global constraints that couple cellular behaviors. The model predicts a transition from a densely interconnected network to a sparse, motif-dominated state as glucose and oxygen consumption levels shift. We characterize such a crossover both numerically, revealing a power-law decay in the cluster-size distribution at the critical transition, and analytically, by computing the critical line through a mean-field approximation based on percolation theory. By comparing empirical data with theoretical predictions, we find that populations evolve towards the sparse phase by remaining near the crossover point between these two regimes. These findings offer new insights into the collective organization driving the adaptive dynamics of cell populations.

We investigate how crystal symmetry tailors the non-Hermitian electro-optic effect driven by the Berry curvature dipole. Specifically, we demonstrate the critical influence of the material's point group symmetry and external electric biases in shaping this effect, leading to current-induced optical gain and non-reciprocal optical responses. Through a symmetry-based analysis of the crystallographic point groups, we identify how different symmetries affect the electro-optic response, enabling the engineering of polarization-dependent optical gain without the need for gyrotropic effects. In particular, we demonstrate that the non-Hermitian electro-optic response in a broad class of crystals is characterized by linear dichroic gain. In this type of response, the eigenpolarizations that activate the gain or dissipation are linearly polarized. Depending on the specific symmetry point group, it is possible to achieve gain (or dissipation) for all eigenpolarizations or to observe polarization-dependent gain and dissipation. Weyl semimetals emerge as promising candidates for realizing significant non-Hermitian electro-optic effects and linear dichroic gain. We further examine practical applications by studying the reflectance of biased materials in setups involving mirrors, demonstrating how optical gain and attenuation can be controlled via symmetry and bias configurations.

The insulating state of matter can be probed by means of a ground state geometrical marker, which is closely related to the modern theory of polarization (based on a Berry phase). In the present work we show that this marker can be applied to determine the metal-insulator transition in disordered systems. In particular, for non-interacting systems the geometrical marker can be obtained from the configurational average of the norm-squared one-body density matrix, which can be calculated within open as well as periodic boundary conditions. This is in sharp contrast to a classification based on the static conductivity, which is only sensible within periodic boundary conditions. We exemplify the method by considering a simple lattice model, known to have a metal-insulator transition as a function of the disorder strength and demonstrate that the transition point can be obtained accurately from the one-body density matrix. The approach has a general {\it ab-initio} formulation and can be applied to realistic disordered materials by standard electronic structure methods.

We present a systematic study of bound relations between different electronic properties of magnetic crystals: electron density, effective mass, orbital magnetization, localization length, Chern invariant, and electric susceptibility. All relations are satisfied for a group of low-lying bands, while some remain valid for upper bands. New results include a lower bound on the electric susceptibility of Chern insulators, and an upper bound on the sum-rule part of the orbital magnetization. In addition, bounds involving the Chern invariant are generalized from two dimensions (Chern number) to three (Chern vector). Bound relations are established for metals as well as insulators, and are illustrated for model systems. The manner in which they approach saturation in a model Chern insulator with tunable flat bands is analyzed in terms of the optical absorption spectrum.

Here we consider the anharmonic oscillator that is a dynamical system given by $y_{xx}+\delta y^{n}=0$. We demonstrate that to this equation corresponds a new example of a superintegrable two-dimensional metric with a linear and a transcendental first integrals. Moreover, we show that for particular values of $n$ the transcendental first integral degenerates into a polynomial one, which provides an example of a superintegrable metric with additional polynomial first integral of an arbitrary even degree. We also discuss a general procedure of how to construct a superintegrable metric with one linear first integral from an autonomous nonlinear oscillator that is cubic with respect to the first derivative. We classify all cubic oscillators that can be used in this construction. Furthermore, we study the Liénard equations that are equivalent to the anharmonic oscillator with respect to the point transformations. We show that there are nontrivial examples of the Liénard equations that belong to this equivalence class, like the generalized Duffing oscillator or the generalized Duffing--Van der Pol oscillator.

We derive a single-cell level understanding of metabolism in an isogenic cyanobacterial population by integrating secondary ion mass spectrometry (SIMS) derived multi-isotope uptake measurements of Synechocystis sp. PCC6803 with a statistical inference protocol based on Liebig's law of the minimum, the maximum entropy principle, and constraint-based modeling. We find the population is structured in two metabolically distinct clusters: cells optimizing carbon yield while excessively turning over nitrogen, and cells which act reciprocally, optimizing nitrogen yield and excessively turning over carbon. This partition enables partial heterotrophy within the population via metabolic exchange, likely in the form of organic acids. Exchange increases the feasible metabolic space, and mixotrophic cells achieve the fastest growth rates. Metabolic flux analysis at the single-cell level reveals heterogeneity in carbon fixation rates, Rubisco specificity, and nitrogen assimilation. Our results provide a necessary foundation for understanding how population level phenotypes arise from the collective contributions of distinct individuals.

In crystals with broken time-reversal symmetry, zone-center phonons can acquire a finite angular momentum via velocity-dependent forces on the nuclei. Despite having the same order of magnitude as the electron spin angular momentum, the phonon angular momentum can be hard to detect because the frequency splitting is small. Here, by developing a theory of lattice magneto-optical effects in reflection and transmission, we show that infrared magnetic circular dichroism is a sensitive probe of zone-center phonon chirality. We evaluate the infrared magneto-optical Faraday, Kerr, and circular-dichroism spectra of CrI$_3$ from time-dependent density-functional theory in the adiabatic local-density approximation. We find sizeable circular dichroism from the infrared-active E$_u$ mode at $\approx 214$ cm$^{-1}$, even though the calculated splitting is only 0.22 cm$^{-1}$.

The shift current is part of the second-order optical response of materials with a close connection to topology. Here we report a sign inversion in the band-edge shift photoconductivity of the Haldane model when the system undergoes a topological phase transition. This result is obtained following two complementary schemes. On one hand, we derive an analytical expression for the band-edge shift current in a two-band tight-binding model showing that the sign reversal is driven by the mass term. On the other hand, we perform a numerical evaluation on a continuum version of the Haldane model. This approach allows us to include off-diagonal matrix elements of the position operator, which are discarded in tight-binding models but can contribute significantly to the shift current. Explicit evaluation of the shift current shows that while the model predictions remain accurate in the deep tight-binding regime, significant deviations arise for shallow potential landscapes. Notably, the sign reversal across the topological phase transition is observed in all regimes, implying it is a robust effect that could be observable in a wide range of topological insulators.

The Beyond Ultra-deep Frontier Fields and Legacy Observations (BUFFALO) is a 101 orbit + 101 parallel Cycle 25 Hubble Space Telescope Treasury program taking data from 2018-2020. BUFFALO will expand existing coverage of the Hubble Frontier Fields (HFF) in WFC3/IR F105W, F125W, and F160W and ACS/WFC F606W and F814W around each of the six HFF clusters and flanking fields. This additional area has not been observed by HST but is already covered by deep multi-wavelength datasets, including Spitzer and Chandra. As with the original HFF program, BUFFALO is designed to take advantage of gravitational lensing from massive clusters to simultaneously find high-redshift galaxies which would otherwise lie below HST detection limits and model foreground clusters to study properties of dark matter and galaxy assembly. The expanded area will provide a first opportunity to study both cosmic variance at high redshift and galaxy assembly in the outskirts of the large HFF clusters. Five additional orbits are reserved for transient followup. BUFFALO data including mosaics, value-added catalogs and cluster mass distribution models will be released via MAST on a regular basis, as the observations and analysis are completed for the six individual clusters.

We study the role of hopping matrix elements of the position operator $\mathbf{\hat{r}}$ in tight-binding calculations of linear and nonlinear optical properties of solids. Our analysis relies on a Wannier-interpolation scheme based on \textit{ab initio} calculations, which automatically includes matrix elements of $\mathbf{\hat{r}}$ between different Wannier orbitals. A common approximation, both in empirical tight-binding and in Wannier-interpolation calculations, is to discard those matrix elements, in which case the optical response only depends on the on-site energies, Hamiltonian hoppings, and orbital centers. We find that interatomic $\mathbf{\hat{r}}$-hopping terms make a sizeable contribution to the shift photocurrent in monolayer BC$_2$N, a covalent acentric crystal. If a minimal basis of $p_z$ orbitals on the carbon atoms is used to model the band-edge response, even the dielectric function becomes strongly dependent on those terms.