# Topological transitions for lattice bosons in a magnetic field

See allHide authors and affiliations

Edited by Allan H. MacDonald, University of Texas at Austin, Austin, TX, and approved September 12, 2011 (received for review July 7, 2011)

## Abstract

The Hall response provides an important characterization of strongly correlated phases of matter. We study the Hall conductivity of interacting bosons on a lattice subjected to a magnetic field. We show that for any density or interaction strength, the Hall conductivity is characterized by an integer. We find that the phase diagram is intersected by topological transitions between different values of this integer. These transitions lead to surprising effects, including sign reversal of the Hall conductivity and extensive regions in the phase diagram where it acquires a negative sign, which implies that flux flow is reversed in these regions—vortices there flow upstream. Our findings have immediate applications to a wide range of phenomena in condensed matter physics, which are effectively described in terms of lattice bosons.

The Hall response is a key theoretical and experimental tool for characterizing emergent charge carriers (1) in strongly correlated systems, ranging from high temperature superconductors (2–4) to the quantum Hall effect (5). In this paper, we study the Hall conductivity of strongly correlated bosons on a lattice. We find that the entire phase diagram of such systems can be characterized using an integer *p*, and inevitably contains topological transitions between regions with different values of *p*. These observations allow us to calculate the Hall conductivity throughout the whole phase diagram, and we show that they lead to surprising consequences, such as sign reversals of the Hall conductivity. The model we study describes a wide range of systems in condensed matter physics, to which our results have immediate implications. Examples are cold atoms on optical lattices (6, 7), Josephson junction arrays (8), granular superconductors (9, 10), and perhaps even high temperature superconductors such as the underdoped cuprates (11–13).

In the absence of disorder and at weak magnetic fields the Hall conductivity of bosonic systems is dominated by the flow of superfluid vortices. For a continuum (Galilean invariant) superfluid, vortex flow gives a Hall conductivity that is proportional to the ratio of the particle density and the applied magnetic field. We find that on the lattice, vortex dynamics is strongly modified. As a result the Hall conductivity is characterized, in addition to the particle density, by the integer *p*. We show how emergent particle-hole symmetry points in the ground-state phase diagram necessarily lead to a nontrivial behavior of this integer, and we discuss the topological transitions between different *p* sectors. As we show, these transitions are attributed to degeneracies in the many-body spectrum, which serve as sources for the Berry curvature.

Specifically, we focus on the conventional Bose–Hubbard model (14) in two dimensions. We restrict our study to a dissipationless system, at zero temperature and without disorder. Within the phase diagram of this model we find large parameter regions corresponding to a negative Hall conductivity, *σ*_{xy} < 0, and reversed vortex motion where vortices flow upstream (see Fig. 1). We discuss methods to directly test these results in cold atom systems where the neutral atoms are subjected to synthetic magnetic fields introduced through rotation or phase imprinting (15, 16).

## Hall Conductivity and Vortex Motion

We begin by giving a semiclassical description of vortex dynamics in bosonic systems. A vortex moving with respect to a current experiences a force arising from the interaction of the velocity field of the vortex with the external current. This hydrodynamical force is called the Magnus force, and it acts perpendicularly to the current, as depicted in Fig. 2*A*. Similarly, a superfluid vortex (of unit vorticity) in two dimensions experiences a force [1]where *n*_{s} is the number density of superfluid bosons, and **v**_{s} is their velocity. The unit vector is a normal to the plane.

The force **F**_{M} in Eq. **1** arises from the dynamical phase (time integral of the energy) in a Lagrangian describing the superfluid. Such a Lagrangian necessarily also contains a term corresponding to the Berry phase picked up by the vortex motion. The Berry phase acquired by a vortex moving around a loop of area *S* is given by 2*παS*, where *α* is a proportionality factor that depends on the microscopic details of the Hamiltonian. Therefore, an equation of motion for the vortex leading to dissipationless flow is linear in the vortex velocity and given by (17) [2]

Eq. **2** can also be understood from the perspective of momentum balance. A moving vortex imprints a phase discontinuity on the superfluid wave function. The Josephson relation Δ*μ* = *ℏ*∂_{t}Δ*φ* connects the resulting chemical potential to the time derivative of the relative phase difference (see Fig. 2*B*). The chemical potential drop will be balanced by a flow of particles, which results in momentum transfer from the particles to the moving vortex, perpendicular to the vortex velocity **v**_{v} and proportional to its magnitude. The proportionality factor *α* relates the change in the system’s momentum to the vortex velocity.

The Hall conductivity can be related to the drift velocity of a vortex. From Eqs. **1** and **2** we get [3]In a system with a low density of vortices we can neglect the effects of vortex–vortex interactions. Considering strictly dissipationless flow, we obtain from the Josephson relation a semiclassical expression for the Hall conductivity [4]where *n*_{v} is the density of vortices and *q* is the boson charge.

In systems with Galilean invariance, in a reference frame moving at the vortex velocity, there should be no forces acting on the vortex, requiring **v**_{v} = **v**_{s} and therefore setting *α* = *n*_{s} (18, 19). This relation is modified in the presence of a lattice, as we discuss now.

We consider the standard model for interacting bosons on a lattice (14) [5]where creates a boson on site **r**, *t* is the hopping amplitude, *U* is the on-site repulsion; is the phase factor due to applied gauge field **A**. We work in units where *ℏ* = *c* = 1, and likewise we set the lattice constant *a* = 1.

We first note that vortices live on the center of the plaquettes of the lattice. We define 2*πα* to be the Berry phase acquired by moving a vortex around a lattice site (i.e., around a plaquette of the dual lattice).* In order to compute *α* we explicitly construct operators and that translate a vortex by one lattice constant in the *x* and *y* directions. We show that they obey the commutation relation (see *SI Text* for a full derivation) [6]Here, is the total boson number operator, and *N* is the number of sites. We denote the particle filling by *n*_{b} = *N*_{b}/*N*. From Eq. **6**, we obtain [7]The integer *p* arises from the 2*π* ambiguity in Eq. **6**.

Eq. **7** can also be understood in terms of momentum balance (20). Following Paramekanti and Vishwanath (21), we note that Eq. **6** implies that when a vortex is transported by Δ*y* sites along , the momentum of the system changes by Δ*P*_{x} = 2*πn*_{b}Δ*y*. At the same time we can integrate Eq. **2** to obtain Δ*P*_{x} = 2*πα*Δ*y*. Combining the two results and taking into account that momentum is only conserved up to a reciprocal lattice vector 2*πp* leads us to Eq. **7**.

The above results, together with Eq. **4**, imply a similar relation for the Hall conductivity, [8]Although Eqs. **4** and **8** are a semiclassical derivation of the Hall conductivity, below we derive an *exact* relation between α and the Hall conductivity for a system containing *one vortex*,^{†} *N*_{v} = 1, [9]

In the remainder of this paper we investigate the relations **7**–**9** throughout the phase diagram of the Bose–Hubbard model. First, we study these relations in the Gross–Pitaevskii and Mott transition limits. In the following sections we then study the transition between different *p* sectors in the hard-core boson limit. Finally, we complete the phase diagram using numerical calculations.

## Low-Energy Limits

We start by discussing low-energy limits of Eq. **5** where a diverging length scale enables the derivation of a continuum low-energy theory. In these limits α and *σ*_{xy} can be deduced directly. We review the derivation of the low-energy theories for weak (*U* ≪ *t*) and strong (*U*≫*t*) interactions. In both cases we start by rewriting Eq. **5** as a coherent state path integral with the following action for the complex valued field *ψ*_{i}: [10]where c.c. is the complex conjugate.

### The Gross–Pitaevskii Limit.

In the weakly interacting limit, the Gross–Pitaevskii (GP) healing length is much larger than the lattice spacing *a*, which enables a straightforward gradient expansion of Eq. **10**. To lowest order in gradients we obtain the continuum action [11]

Using the above expression we can now derive the coefficient α in the Gross–Pitaevskii limit. When written in terms of , the term *ψ*^{∗}∂_{τ}*ψ* in Eq. **11** leads to a purely imaginary , where *m* is a field that counts the winding of the phase *ϑ* (17). Consider the action of a field configuration associated with taking a vortex around a closed loop of area *S*. A little reflection shows that the regions outside the loop do not change the value of *m*, whereas those inside contribute unity per particle. Hence, gives rise to a Berry phase of 2*πin*_{b}*S* (18). This observation fixes [12]

### Around the Mott Insulator.

At strong interactions and integer filling, , the Hamiltonian (Eq. **5**) stabilizes a localized Mott insulating phase with vanishing superfluid fraction (14). In the insulating phase, all sites are occupied by exactly *n*_{b} bosons. Both the addition or removal of a particle is protected by a finite gap. This gap closes at the boundary of the Mott lobes in the phase diagram of Fig. 1, whereby at the lower (upper) boundary of the Mott lobe the hole (particle) gap vanishes. Hence, the tip of the Mott lobe represents a multicritical point where the particle and hole gap close simultaneously (22). In the following, we focus on this multicritical point.

When both the particle and hole gap vanish, an enhanced symmetry in the low-energy sector emerges. Instead of going through the standard procedure of deriving the low-energy theory from microscopic considerations (23, 24), we motivate the effective action via its symmetry properties. We expect the following particle-hole symmetry (PHS) to hold and **A** → -**A** (25). To leading order in powers of we find [13]The gradient expansion leading to an effective continuum theory is controlled by the diverging correlation length close to the second-order phase transition into the Mott insulating state.

A direct consequence of PHS in the *continuum theory* (Eq. **13**) is *σ*_{xy}(**A**) = *σ*_{xy}(-**A**). Together with the Onsager relation *σ*_{xy}(**A**) = -*σ*_{xy}(-**A**) we obtain *σ*_{xy} = 0. This result can also be understood in terms of vortex motion. As opposed to the Gross–Pitaevskii action, the continuum theory (Eq. **13**) is real. Hence it does not give rise to any Berry phase when a vortex is moved around a closed loop, and we conclude [14]

Starting from the PHS points, we expect to find lines with *α* = 0 in the *μ*/*U* - *t*/*U* phase diagram of the Bose–Hubbard model. From Eq. **7**, on the other hand, we know that at a fixed density, α can only change by an integer, which leads to the conclusion that the lines with *α* = 0 are bound to lines of integer fillings in the phase diagram (see Fig. 1).

## Hard-Core Bosons Limit

We now consider the limit *t*/*U* → 0 and *μ*/*U* → *m*, where *m* is an integer. In Fig. 1, these limits lie in-between two Mott lobes. The two states with *m* and *m* + 1 bosons per site are degenerate single-site states of the Hamiltonian [**5**]. States with different fillings are separated by a gap of order *U* and do not appear in the low-energy theory.

We use a Schrieffer–Wolf transformation (26) to project the Hamiltonian [**5**] onto the subspace with only *m* and *m* + 1 bosons per site. The resulting Hamiltonian corresponds to hard-core bosons (HCB) and can be written using spin-half operators; is the on-site number operator, and () raises (lowers) the occupation by one. At zeroth order in *t*/*U*, the HCB Hamiltonian is given by [15]where h.c. represents the hermitian conjugate.

The HCB Hamiltonian Eq. **15** has an emergent charge conjugation symmetry. One defines the unitary transformation [16]*C* transforms particles into holes (i.e., ), and [17]At half-filling for the hard-core bosons, the Hamiltonian [**15**] is independent of *μ* and hence Eq. **17** implies invariance under **A** → -**A**. Hence, the Onsager relation *σ*_{xy}(**A**) = -*σ*_{xy}(-**A**) implies that for half-integer fillings () [18]

Note that the situation at the HCB limits and at the tip of the Mott lobes are qualitatively different. The Hall conductivity at the tip of the Mott lobe vanishes due to a zero crossing of α when *p* = -*n*_{b}. In the HCB limit, the integer *p* jumps exactly at . In other words, in the first case it reflects a particle-hole-symmetry between *n*_{b} - 1 and *n*_{b} + 1, whereas in the latter, the symmetry connects *n*_{b} - 1 and *n*_{b} at . The symmetry at the HCB limits has a remarkable consequence for *σ*_{xy} in the full phase diagram of the model, as we shall now show.

## Away from the Hard-Core Boson Limit

We now consider the effect of a finite but small value of *t*/*U*. Second-order processes in which a virtual excitation with an on-site occupation of *m* - 1 or *m* + 2 bosons are created lead to corrections to the HCB Hamiltonian [**15**] of order *t*^{2}/*U*. Taking into account all the different processes, we obtain up to irrelevant renormalizations of the parameters in [19]where denote sites **r** and **r**^{′} which are nearest neighbors of site **r**^{′′}, and *ϵ*_{m} = (*m* + 1)(*m* + 2)*t*^{2}/*U*.

The new terms in Eq. **19** break the charge conjugation symmetry. Therefore, for finite *t*/*U*, the Hall conductivity at exactly half integer filling does not vanish; below we calculate it in the limit of small *t*/*U*.

We consider the model Eq. **19** on a torus of size *N* = *L*_{x}*L*_{y}, with *N* even. The gauge field **A** describes a uniform flux penetrating the surface of the torus. We take the total flux to be one flux quantum, which induces one vortex into the system. An important gauge invariant quantity described by the gauge field are the two Wilson line functions (27) [20]We define Θ_{x} = Φ_{x}(*y* = 0) and Θ_{y} = Φ_{y}(*x* = 0). Changing the values of Θ_{x} and Θ_{y} corresponds to threading Aharonov–Bohm (AB) fluxes through the two holes of the torus (27).

The Hall conductivity at zero temperature for a general many-body Hamiltonian can be calculated by integrating the Berry curvature (28) [21] is given by [22]Here Ψ_{0}(**Θ**) is the many-body ground-state wave function that depends on the AB fluxes through the holes of the torus.

Remarkably, the Hall conductivity in the presence of one vortex can be calculated analytically at half-filling. The key ingredients are degeneracies in the spectrum that occur for *t*/*U* = 0 (27, 29) and serve as point (monopole) sources for the Berry curvature (30).

To understand the degeneracies, we consider an effective Hamiltonian for the vortex hopping between dual lattice sites. As shown in refs. 27 and 29, this Hamiltonian is given by [23]where creates a vortex on a dual lattice site, and *t*_{V} ≈ *t*. The dual gauge field’s flux is given by the boson density, ∇ × **A**^{D} = Φ^{D} = 2*πn*_{b}. The potential *U*(**R** - **R**_{V}) for the vortex position arises due to the fact that the Wilson lines [**20**] break translational symmetry on the torus. In fact, as shown in refs. 27 and 29, for one flux quantum penetrating the surface of the torus, all translational symmetries are absent, and the potential *U*(**r** - **R**_{V}) acquires its minimum at a point **R**_{V} for which the Wilson lines both take on the value *π*.

If the point **R**_{V} lies on a site of the direct lattice, the eigenstates of (in a symmetric gauge) can be written as *ψ*(**R** - **R**_{V}) = *f*(|**R** - **R**_{V}|)*e*^{imφ(R-RV)}. Here *φ*(**R** - **R**_{V}) denotes the angle between **R** - **R**_{V} and the *x* axis, and *m* = 0, ± 1, ± 2. At half-filling, the average dual flux per plaquette is and the ground state is doubly degenerate with *m* = 0,1. The two states (*m* = 0) and (*m* = 1) represent states with clockwise and counterclockwise vortex currents, respectively, as depicted in Fig. 3*B*. Note that this twofold degeneracy occurs for *N* distinct values of **Θ**.

To calculate *σ*_{xy}, we need to analyze the spectrum around the *N* degeneracy points. Around these points, the Hamiltonian restricted to the and basis is of the form *H*_{V} = **h**·** σ**, where

**are the pauli matrices. To find**

*σ***h**, we first notice that if is a degeneracy point, tuning away from it by

**Θ**=

**Θ**

^{0}+ Δ

**Θ**, moves

**R**

_{V}as (27) [24]where

*α*,

*β*=

*x*,

*y*are not summed. Thus, tuning away from

**Θ**

^{0}breaks the degeneracy between the two ground states and , as it shifts the minimum of the potential

*U*(

**R**-

**R**

_{V}). Second, the terms ∝

*ϵ*

_{m}in [

**19**] lift the degeneracy even at

**Θ**=

**Θ**

^{0}; the assisted hopping terms through

**R**

_{V}(blue arrow in Fig. 3

*B*) favor over .

Together, these two effects give rise to the following low-energy Hamiltonian for each degeneracy point (see *SI Text* for details), [25]where the energy scales appearing above are , and .

We use Eqs. **21**–**25** to calculate the Hall conductivity for *one vortex*. Consider the Hamiltonian of Eq. **19** where we let the parameter *ϵ*_{m} take on both negative and positive values. For 0 < *ϵ*_{m}/*t* ≪ 1, the many-body ground-state Ψ_{+}(**Θ**,*ϵ*_{m}) is nondegenerate, and likewise Ψ_{-}(**Θ**,*ϵ*_{m}) for -1 ≪ *ϵ*_{m}/*t* < 0. For *ϵ*_{m} = 0, the two states become degenerate at a set of *N* values of **Θ** space. The Berry connection of the ground-state manifold at *ϵ*_{m} = 0 is given by [26]and must satisfy due to particle-hole symmetry at *ϵ*_{m} = 0. As a result, [27]

Next, we consider the space of Θ_{x},Θ_{y},*ϵ*_{m}, which has the topology of a thick torus, as depicted in Fig. 4. We are interested in the integral of the Berry curvature on the surfaces with *ϵ*_{m} > 0 and *ϵ*_{m} < 0, which yield *σ*_{xy}(*ϵ*_{m} > 0) and *σ*_{xy}(*ϵ*_{m} < 0), respectively. At the same time, an analog of Gauss’s law for implies that the integral of on these surfaces counts the number of sources for (30, 31). These are just the degeneracy points discussed above, which are all described by Eq. **25** and therefore correspond to sources with charge , which leads to [28]Combining Eqs. **27** and **28** gives [29]

Before concluding this section, we note that Eq. **24** leads to an exact relation between α and the Hall conductivity of one vortex. From Eq. **24**, the Berry phase for moving a vortex around a plaquette is given by [30]where the contour defines a square of size 2*π*/*L*_{x} × 2*π*/*L*_{y} in flux space, and is the surface it bounds. Therefore, for one vortex, [31]

Consider again the phase diagram of the Bose–Hubbard model. From the above discussion, we conclude that the transitions lines between two integer values of *p* emanate from the HCB points and move to higher densities with increasing *t*/*U*. These lines all correspond to changes of the integer *p* by unity. The PHS lines emanating from the neighboring Mott lobe tips *terminate* at these topological transition lines (see Fig. 1). Together, these two types of lines define regions with negative α and Hall conductivity.

## Evolution of the Transition Lines

We numerically calculate the Chern number (21) for one vortex, to obtain the behavior of the integer *p* in the full parameter regime of the Bose–Hubbard model. We use a Lanczos algorithm (32) to find the ground-state wave function Ψ_{0}(**Θ**) for different AB fluxes. Using a standard procedure (33) to numerically integrate the Berry curvature [**22**] we obtain the Hall conductivity for different values of *t*/*U* and *n*_{b}.

In Fig. 5, we show the results obtained for a 3 × 3 cluster (see *SI Text*). We indicate which integer *p* describes the Hall conductivity in Fig. 5*B*. Regions in the phase diagram where *σ*_{xy} < 0 are marked with yellow hatches. Fig. 5*A* shows a trace of *σ*_{xy} for different particle numbers at *t*/*U* = 0.2. As expected from the calculation at half-filling, the transition lines between integer values of *p* move to higher densities for *t*/*U* > 0. Remarkably, the transition lines intersect the integer density line at increasing values of *t*/*U* for higher densities. As a direct consequence, the area of negative Hall conductivity *increases* for higher densities, which is in contrast to the *decreasing* extent of the Mott insulating phases indicated by the yellow bars in Fig. 5*B*. This surprising behavior of the Bose–Hubbard model is explained below.

In order to see at which values of *t*/*U* a sign change of α should be expected at integer fillings, we consider the healing length *ξ*_{GP}, which sets the size of a vortex. For , the size of a vortex is much smaller then the lattice spacing *a* (34), and the Bose–Hubbard model maps onto the quantum rotor model (10, 35), which has an emergent PHS at integer filling (17, 36). We have seen that PHS implies *α* = 0. The dependence of *ξ*_{GP} on the *mean-field* interaction *Un*_{b} therefore explains the growing extent of the negative Hall conductivity. This behavior has to be contrasted to the bosonic enhancement factors in the hopping terms which facilitate the melting of the Mott insulator and lead to smaller Mott lobes at increasing densities.

Finally, in Fig. 1 we present the numerical results as a function of *t*/*U* and *μ*/*U*. To translate from the results at fixed density *n*_{b} to a fixed chemical potential *μ*, we use a standard mean-field approach (37).

## Experimental Verification

Our results can be directly verified experimentally using two complimentary approaches. First, α, and in particular its sign, can be determined by directly measuring the vortex velocity. Second, the transverse response of collective modes of a superfluid can be used to measure its Hall conductivity. In the following, we discuss methods to implement these approaches in a system of cold atoms.

In the presence of a superfluid current, the vortex velocity is related to α via Eq. **3**. In a cold-atoms setup, the direction and speed of the vortex flow can be measured with in situ imaging techniques (38), which image the density profile near the vortex core. However, in the strongly interacting regime, the vortex core is on the order of the lattice spacing *a*. Hence, to make the vortex visible in the density profile, the system parameters have to be ramped into the weakly interacting regime before imaging.

The sign change of the Hall conductivity can also be detected by studying collective modes of an atomic superfluid in a harmonic trap. We illustrate the effect of the Hall conductivity on the dynamics of the “sloshing” mode (39). This mode is excited by suddenly displacing the minimum of the harmonic trap with respect to the atomic cloud. In the absence of a Hall response, this displacement causes the center of the atomic cloud to oscillate in the trap along the displacement axis, as shown in Fig. 6. A nonzero Hall conductivity yields a transverse force, which deflects the motion of the atomic cloud (see Fig. 6). The direction in which the center of the cloud is deflected depends on the sign of *σ*_{xy}. Similar effects can be studied in collective modes with higher angular momentum, such as the quadrupole mode (39). We note that density inhomogeneities due to the trapping potential will lead to a distortion of the atomic cloud during its oscillation in the trap. However, the effect of a sign change in *σ*_{xy} should be visible even in the presence of these distortions.

## Discussion and Outlook

In this paper, we focused on vortex dynamics for the Bose–Hubbard model. We mapped the sectors corresponding to different integers *p*, which characterizes the Hall conductivity and vortex motion throughout the phase diagram. We found that, close to the Mott insulating phases, the sign of *σ*_{xy} is reversed and vortices flow against the applied current.

Our results are obtained neglecting vortex–vortex interactions or disorder. We note, however, that there are an infinite number of particle-hole symmetric points in the zero temperature phase diagram: at the tip of every Mott lobe, and in between two adjacent Mott lobes. We saw that the latter necessarily slice the full phase diagram into an infinite number of different *p* sectors. This underlying structure cannot be removed by the inclusion of vortex–vortex interactions or disorder. However, the transition lines are expected to change their exact location and to be smoothed out by these effects, as well as by finite temperature.

Incidentally, reversal of the Hall conductivity has been repeatedly measured in many strongly correlated electronic materials, including high-temperature superconductors (e.g., refs. 2–4). These experiments are beyond the direct applicability of our model. An extension to treat these materials is an interesting future direction. As discussed above, a clean verification of our predictions is possible in systems of cold atoms.

## Acknowledgments

We thank Assa Auerbach, Ehud Altman, Joseph Avron, Hans-Peter Büchler, Olexi Motrunich, Ian Spielman, and Ady Stern for fruitful discussions. Special thanks to Daniel Podolsky for his enlightening comments. N.H.L. acknowledges support by The Gordon and Betty Moore Foundation through Caltech’s Center for the Physics of Information, National Science Foundation Grant PHY-0803371, and the Israel Rothschild Foundation. S.D.H. acknowledges support by the Swiss Society of Friends of the Weizmann Institute of Science. This research was supported in part by the National Science Foundation under Grant PHY05-51164.

## Footnotes

- ↵
^{1}To whom correspondence should be addressed. E-mail: sebastian.huber{at}weizmann.ac.il.

Author contributions: S.D.H. and N.H.L. designed research, performed research, and wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1110813108/-/DCSupplemental.

↵

^{*}The Berry phase acquired by moving a vortex around a dual plaquette is related to*α*in the semiclassical theory (see Eq.**2**), by dividing by the area of a plaquette*a*^{2}.↵

^{†}Note that the exact relation in Eq.**9**is related to Eq.**8**by setting*n*_{v}= 1/*N*.

## References

## Citation Manager Formats

## Article Classifications

- Physical Sciences
- Physics