author:
Astronomy Department, University of Cape Town
E-mail: ; Michelle E. Cluver
IPAC, California Institute of Technology
E-mail: ; Tom H. Jarrett
IPAC, California Institute of Technology
E-mail: ; Patrick A. Woudt
Astronomy Department, University of Cape Town
E-mail:
title: Properties of the extremely HI-massive galaxy HIZOA J0836-43

# Introduction

The number of galaxies with HI-masses above a characteristic value, $M^*_{\rm{HI}} = 6\cdot10^9$`<!-- -->`{=html} declines exponentially (Zwaan et al. 2005). The most massive galaxies are very scarce, their number density is ill-constrained, and little is known about their properties.

Hierarchical structure formation predicts these galaxies to be forming only recently ($z < 1$; Mo et al. 1998) although models overpredict their current population ($z = 0 - 1$; Renzini 2007). Local massive and gas-rich galaxies are thus expected to be in an early stage of formation. But the very few known ones - like Malin 1 (Impey & Bothun 1989) - are rather quiescent, providing few clues about their past or future evolution.

A detailed investigation of the recently discovered, fairly nearby ($v_{hel}$=10689) galaxy HIZOA J0843-36 (Donley et al. 2006), one of the most massive HI-galaxies known to date, will therefore provide a valuable local probe to be put into context with the forthcoming higher redshift pathfinder surveys.

HIZOA J0843-36 was discovered with the Parkes radio telescope in the deep systematic HI survey of the Zone of Avoidance (ZOA; e.g. Kraan-Korteweg et al. 2005). Follow-up ATCA and AAT radio observations confirmed it to be amongst the most massive spiral galaxies (Donley et al. 2006). It contains 7.5$\times10^{10}$`<!-- -->`{=html} of  gas, has a total dynamical mass of 1.4$\times10^{12}$`<!-- -->`{=html} and a 20-cm derived star formation rate of $\sim$`<!-- -->`{=html}35/yr. It also has a prominent bulge in the near-infrared (bulge-to-disk ratio in the $K_{s}$ band of $\sim 0.80$) which lies central to its enormous, rapidly rotating HI disk. It is amongst the most HI-massive spiral galaxies known to date and, at a distance of 148 Mpc, also the nearest. Contrary to other HI-rich galaxies, like Malin 1, it is, however, not a giant low surface brightness galaxy.

This leads to the question: What is happening in this unusual galaxy? Is the galaxy undergoing a starburst, or is it an AGN? Can its environment provide clues about its origin/formation? How did it grow to its current size? Was it through merger/accretion processes?

# Near and Mid-Infrared observations of HIZOA J0836-43

A detailed investigation of HIZOA J0843-3 is rendered difficult due to its location behind the Vela supernova remnant in the ZOA. Lying behind a dust layer of $A_B = 10\fm0$ it is optically invisible. To minimize the effects of foreground extinction, we launched an observational campaign in the infrared of the galaxy and its local environment. The program comprises imaging of an area of $28\arcmin \times 28\arcmin$, centered on the galaxy, obtained by using *Spitzer* with IRAC (3.6, 4.5, 5.8, 8.0$\micron$) and MIPS (24 .0, 70, 160$\micron$), as well as IRS high and low resolution spectroscopy of the galaxy itself. We also performed a deep imaging survey ($J H K_s$) of $1\fdg5 \times 1\fdg5$ with the Infrared Survey Facility (IRSF) at the 1.4 m Japanese telescope in Sutherland (SAAO), see Cluver (2009) for details.

## Infrared properties of HIZOA J0836-43

The near-infrared light profiles of the galaxy appear typical for an early spiral or lenticular with its prominent bulge, while the radio observations suggest recent, active star formation. The new data find a NIR – MIR Spectral Energy Distribution that is consistent with an Sc morphology (from GRASIL modeling). However, the total IR luminosity of $L_{\rm{TIR}}$=$1.2\times{10}^{11}$`<!-- -->`{=html} and star-formation rate in the far infrared (FIR) of $\sim 21$`<!-- -->`{=html}/yr is much higher than that of typical disk spirals and is more consistent with the properties of local Luminous Infrared Galaxies (LIRGs). The galaxy also has substantial cold molecular gas ($M_{\rm{H_{2}}}$=$1.3\times{10}^{10}$`<!-- -->`{=html}) as estimated from the FIR, but little warm dust. The galaxy possesses a prominent bulge of evolved stars and a total stellar mass of 4.4$\times10^{10}$`<!-- -->`{=html} (Cluver et al. 2009).

In summary, we have a starbursting galaxy (the spectrum shows no signs of AGN activity) with a prominent bulge comprised of evolved stars. But its specific star formation rate ($\sim$`<!-- -->`{=html}0.5 $\rm{Gyr}^{-1}$ ) implies an early stage of stellar mass building. With its huge reservoir reservoir of gas it could double its stellar mass in about 2 Gyr. Based on the stellar surface brightness profile and emission from tracers of star formation, the galaxy seems to be undergoing inside-out star formation similar to processes observed at higher redshift ($z\sim 0.7$) where gas fractions of disks were higher compared to local galaxies (e.g. Bell et al. 2005, Pérez-González et al. 2005). Despite its extreme  mass, both the near- and mid-infrared luminosity appear consistent with expectations from /luminosity relations in these bands. This suggests relatively "normal" evolution in HIZOA J0836-43, despite lying at the extreme high end of the relation.

### Environment of HIZOA J0836-43

The question that remains unanswered is: did HIZOA J0843-36 grow to its current size through cannibalism in a dense environment or accretion along filaments? Early NIR observations found two possible companions, further supported by the subsequent Spitzer observations.

The deep NIR IRSF survey of 2.24$\Box^\circ$ around the massive galaxy revealed 404 galaxies, most previously unknown. Of the 44 2MASX object in the survey area, only 16 coincided with our galaxies. The deeper and higher resolution IRSF images found the 28 other 2MASX objects to either be blended (unresolved) star images or Galactic objects. According to photometric redshifts (Jarrett et al. 2009, in prep.), only a handful of the IRSF galaxies lie at the distance of the massive galaxy.

The absolute magnitude distribution of the galaxies that are – according to their photometric redshifts – located within a volume of radius $r = 10$ Mpc around HIZOA J0836-43 are compared to four regions of varying galaxy density at similar redshifts. This includes a field sample, a supercluster filament, a fork apex, and a cluster (Hercules). The comparison reveals that the population of galaxies in the surroundings of the massive galaxy is constituted of an overdensity of low stellar mass galaxies ($M^o_{K{_s}} \ga -23\fm0$) with only a handful of brighter galaxies. Even fewer than in the comparison *field* sample.

Hence, the massive galaxy HIZOA J0836-43 lives in a low $L^*$ density environment. That may have allowed its formation and survival, enabling it to evolve into the unusual LIRG starburst galaxy we observe today. It also seems consistent with the observation that HI-massive galaxies are preferentially found in low density regions.

Inspection of the global larger-scale structures in which HIZOA J0836-43 is embedded finds this galaxy to lie in a slightly enhanced density of fainter galaxies within a (newly identified) filamentary structure encircling a void (Cluver 2009). The conditions appear favorable for the galaxy to have grown this large HI disk, both through accretion of gas due to minor merging, as well as infall of gas along the filament.

# Next Steps

This is the only known HI-massive galaxy undergoing an infrared-luminous starburst. It appears to be in the process of building a disk, fueled by its large gas reservoir. We have obtained millimeter observations with MOPRA to complement the infrared study of the galaxy's star formation phase and have detected the very broad CO (1-0) line of HIZOA J0836-43. We find a total molecular gas mass (H$_2$ + Helium, where the He fraction is 1.38) of 4.5$\times{10}^{9}$`<!-- -->`{=html} corresponding to a H$_2$-to-HI gas-mass fraction of $\sim$`<!-- -->`{=html}6%. We, however, require deeper observations to obtain the sensitivity to detect the low velocity CO gas and map the extent of molecular gas in the disk. This will provide clues regarding the galaxy's past and present evolution.

A future deep HI survey with MeerKAT (or ASKAP) will shed further light on the environment of the galaxy with its association of low stellar mass, but possibly HI-rich galaxies. Previous simulations of MeerKAT performance in the ZOA in the Great Attractor region (Kraan-Korteweg et al. 2009) have proven its effectiveness in answering these questions on reasonable timescales.

Such an HI-survey will furthermore give substance to the findings by Kraan-Korteweg (in prep.) of the existence of a massive overdensity lurking in the ZOA in Vela at $v_{hel} \sim 20000$`<!-- -->`{=html}. This overdensity is a complete surprise. Some suggestions of its existence are seen in the reconstructed density maps based on the MASX Redshift Survey (Erdogdu et al. 2006). If confirmed, this would have major implications for the the dynamics in the nearby Universe. It may well help explain the large-scale flows towards that direction found by Hudson et al. (2009).

**Acknowledgment:** Financial support from the National Research Foundation and the University of Cape Town is gratefully acknowledged.


abstract: Collective coherent light scattering by polarizable particles creates surprisingly strong, long range inter-particle forces originating from interference of the light scattered by different particles. While for monochromatic laser beams this interaction decays with the inverse distance, we show here that in general the effective interaction range and geometry can be controlled by the illumination bandwidth and geometry. As generic example we study the modifications inter-particle forces within a 1D chain of atoms trapped in the field of a confined optical nanofiber mode. For two particles we find short range attraction as well as optical binding at multiple distances. The range of stable distances shrinks with increasing light bandwidth and for a very large bandwidth field as e.g. blackbody radiation we find a strongly attractive potential up to a critical distance beyond which the force gets repulsive. Including multiple scattering can even lead to the appearance of a stable configuration at a large distance. Such broadband scattering forces should be observable contributions in ultra-cold atom interferometers or atomic clocks setups. They could be studied in detail in 1D geometries with ultra-cold atoms trapped along or within an optical nanofiber. Broadband radiation force interactions might also contribute in astrophysical scenarios as illuminated cold dust clouds.
address: Institute for Theoretical Physics, University of Innsbruck, Technikerstraße 25, A-6020 Innsbruck, Austria
author: D. Holzmann and H. Ritsch
title: Tailored long range forces on polarizable particles by collective scattering of broadband radiation

# Introduction

Light scattering from point like particles is connected to momentum and energy exchange between the particles and the field. The interference of the fields scattered from different particles in an extended ensemble leads to important modifications of the total force on the ensemble  and introduces long range inter-particle light-forces . These forces, which originate from coherent scattering stay sizable even if the light fields are far detuned from any optical resonance . The full coupled nonlinear interaction in a cloud thus leads to a very rich and complex dynamics  including trapping, optical binding  and selfordering . Interestingly, many of the key physical effects can already be found and studied in effective 1D geometries . One particularly interesting example are atoms coupled to 1D optical micro structures  as e. g. an optical nanofiber, where even a single atom can significantly modify light propagation .

In a milestone experiment, Rauschenbeutel and coworkers recently managed to trap cold atoms alongside a tapered optical nanofiber , a result which has been repeated and improved in several new setups . As a key consequence of the strong atom-fiber coupling, optical dipole interaction and light forces between two atoms are now mediated over basically the whole fiber length. Corresponding calculations of the collective particle dynamics exhibit enhanced inter-particle forces supporting the formation of periodical self-ordered regular arrays , where light and motion is strongly coupled and correlated .

Even for large ensembles, illumination with monochromatic laser light generates a translation invariant 1D geometry. To lowest order this implements equal strength infinite range coupling of each particle to all the others . At sufficiently low kinetic temperature this induces crystallization of the particles and light fields with characteristic phonon excitations . It has already been shown, that adding a second laser frequency enhances particle-particle interaction at controllable distances . Note that using circularly polarized light, asymmetric chiral scattering and interaction can be implemented in such systems as well to generate very exotic chiral spin models . Here for simplicity we will assume transverse linear polarization though to keep inversion symmetry.

In this work we generalize this coupled atom-field dynamics to incoherent light with a large bandwidth and no stable phase relation between the different frequency components. Hence the incoming radiation field is space and time translation invariant, but now possesses a finite correlation length. Surprisingly there appear substantial and partly even enhanced inter particle forces, for which we can control the effective interaction length via the spectral shape of the incoming radiation. Note that such broadband finite range forces were recently also predicted and observed in a seminal nano particle experiment with incoherent 3D illumination in a solution . As the corresponding full 3D model is rather complex to solve and understand, we will restrict ourselves to a quasi 1D setup as mentioned above and depicted in allowing an at least partially analytic treatment. Nevertheless we expect the model to exhibit the essential underlying multiple scattering physics. Note that instead of shaping the incoming light spectrum, one can alternatively envisage to tailor the waveguide dispersion relation by nano-structuring to manipulate propagation lengths .

This work is organized as follows. We first generalize the basic definitions and dynamical equations of the proven multiple scattering model for light forces in 1D systems  towards fields with finite bandwidth, orthogonally impinging on small particles in a 1D trap along an optical waveguide structure cf. . For two particles we analytically calculate the resulting inter-particle forces as function of bandwidth and separation. These results are then numerically extended to higher particle numbers studying the self-consistent coupled particle field dynamics for small ensembles and varying bandwidth. As a final case we will study temperature dependent forces in thermal radiation fields.

# Multiple scattering approach for a point particle chain in broadband light

Let us assume an incoming broadband radiation field with all the wavelength components larger than the particle size and detuned from any optical resonance. The particles move in a 1D trap close enough to an optical fiber so that they can coherently scatter light into and out of the propagating fiber modes. For simplicity we assume a nanofiber, which supports only a single transverse mode propagating in each direction and a single peaked frequency distribution for the incident light field with center frequency $\omega_0$ and width $\gamma_0$. Its maximum intensity is $I_0$ so that: $$I(\omega)=\frac{I_0}{\pi}\frac{\gamma_0}{\gamma_0^2+(\omega-\omega_0)^2},
\label{i0}$$ which can be considered as a generic model of a broadband laser field or fluorescence of a large atomic ensemble. The scattered field spectrum by a single point like scatterer will than largely be proportional to the incoming distribution. For very broad light fields as generated from an LED source or ultimately using blackbody radiation with a very broad Planck spectrum, one would of course get corrections from the frequency dependent scattering amplitude $\eta$. We will only address this limit a bit at the end of this work.

## Optical binding forces between two particles

To get some first insight we start with the simplest nontrivial example of two particles at a distance $d=x_2-x_1$ scattering light into the fiber. Inside the fiber the fields scattered from the two particles will interfere and a part of the field scattered by the first particle will be backscattered from the second one and vice versa. As shown in previous work, this situation can be modeled by a beam splitter approach via a $3\times 3$ coupling matrix . Each beam splitter is parametrized by a complex polarizability $\zeta$ with dispersive part $\zeta_r$, absorptive part $\zeta_i$ and an effective scattering amplitude $\eta_0$, which for simplicity we assume to be frequency independent. In this case the multiple scattering can be summed to all orders and still allows for analytic calculations of the total self-consistent fields and forces by integrating over the whole spectrum. Note that as only the product of the spectral polarizability and the illumination spectrum appears in the calculations any frequency dependence can be absorbed into a more complex effective pump distribution.

For two particles and small $\zeta=\zeta_r+i\zeta_i$  as well as small relative bandwidth $\gamma_0\ll\omega_0$, integration over the whole frequency range  yields: $$\begin{aligned}
\fl
\eqalign{F_1 &= \frac{I_{\eta_0}}{\pi c} \int_{-\infty}^{\infty }\frac{\gamma_0}{\gamma_0 ^2+(\omega -\omega_0)^2}\cos \left(\frac{\omega d}{c}\right) \left(1-2\zeta_i\cos\left(\frac{\omega d}{c}\right)-2\zeta_r  \sin \left(\frac{\omega d}{c}\right)\right) d\omega\\
&=\frac{I_{\eta_0}}{c}\left( e^{-\frac{\gamma_0 d}{c}}\cos\left(\frac{\omega_0 d}{c}\right)-\zeta_i\left(1+e^{-2\frac{\gamma_0 d}{c}} \cos\left(\frac{2\omega_0 d}{c}\right)\right)-\zeta_r e^{-2\frac{\gamma_0 d}{c}}\sin\left(\frac{2\omega_0 d}{c}\right) \right).}
\label{F2eq}
\end{aligned}$$

In contrast to the infinite range for a monochromatic field  we find an additional exponential decay of the inter-particle force determined solely by the illumination bandwidth. This behavior is exhibited in figure , where the forces on the first particle for different bandwidths are compared. We still find stable configurations (optical binding) at several different distances, which can be calculated from the zero-points of the force in . Note that as a reference value we also have included the case of an unrealistically broad distribution of width $\gamma_0=\omega_0$ to get some qualitative insight for the infinite bandwidth limit.

For a stable point we need zero average force and that the derivative of the force on the first beam-splitter with respect to distance is positive and negative for the second one. Surprisingly for non-absorbing particles with vanishing imaginary part (negligible absorption) of $\zeta$, $\zeta_i=0$, these stable distances $d_0$ do not depend on the bandwidth and are equal as for a monochromatic case :

$$d_0=\left(\frac{3}{4}+n\right)\lambda_0,~n\in\mathbb{N}.$$

When we also take absorption into account by introducing an imaginary part of $\zeta$, this adds extra outward radiation pressure to the force and shifts the stable distances to larger values. Eventually for too large $\zeta_i$ no stable configuration can be found as shown in figure  and . Note that scattering of light to non propagating field modes simply appears as effective loss in the imaginary part $\zeta_i$.

Once knowing the forces we can also simulate the coupled particle field dynamics for non-equilibrium conditions. Adding some effective friction to particle motion we simply integrate the equations of motion and let the coupled particle field system evolve towards stable configurations. This is shown in figures  and , where such trajectories for different illumination bandwidths $\gamma_0$ are compared. Indeed the simulations confirm the independence of the stable distance on the line width of the radiation.

Note that in all cases the particles act like resonators creating a field intensity maximum in between them, while they themselves do not sit at either a field maximum or minimum as one might expect from simple light shift considerations. As a clear signature of multiple scattering we also get a strong spectral filtering response reminiscent of a Fabry Perot cavity from our particle pair. While some frequencies are strongly confined, others are dominantly transmitted depending on the particle distance. As a consequence the spectral distribution of the intensity emitted into the fiber outside the particle pair is strongly modified from the input. It also substantially changes with varying polarizability $\zeta$ and the line width $\gamma_0$.

This is depicted in some typical examples in figure  and . Unexpectedly the intensities are not spectrally symmetric, which is a consequence of the strong frequency-dependent scattering of the beam-splitters. In this case some frequencies are filtered by particles, while others can pass them. Which frequencies can pass does not only depend on $\zeta$ and $\gamma_0$ but also at which of the stable distances the particles order. Thus, one can directly get information on the distance of the particles by measuring the frequency distribution of the outgoing intensity.

## Forces and selfordering for several beam splitters

In principle determining the field distributions, the forces and stationary states for a larger number of beam splitters is straightforward by multiplying the corresponding scattering and propagation matrices. In practice, useful analytic results can only be obtained by ignoring backscattering and setting $\zeta=0$ as multiple scattering leads to very complex expressions. In contrast to the monochromatic case the particles in steady state here are not equidistantly distributed, which leads to a much more complex transfer matrix. Nevertheless, for $\zeta=0$ the transfer matrix still can be calculated analytically for $N$ particles and integrated over the frequency distribution to yield the following expression for the total force: $$\begin{aligned}
\fl
\eqalign{F_{j,N}&=\frac{I_{\eta_0}}{c\pi}\int_{-\infty}^{\infty} d\omega \frac{\gamma_0}{\gamma_0^2+(\omega-\omega_0)^2}\left(\sum _{l=j}^{N-1} \cos \left(\frac{\omega\sum _{i=j}^l d_i}{c}\right)-\sum _{l=1}^{j-1} \cos \left(\frac{\omega\sum _{i=l}^{j-1} d_i}{c}\right)\right)\\
&=\frac{I_{\eta_0}}{c}\left(\sum _{l=j}^{N-1}e^{-\sum _{i=j}^l\frac{\gamma_0d_i}{c}} \cos \left(\frac{\omega\sum _{i=j}^l d_i}{c}\right)-\sum _{l=1}^{j-1} e^{-\sum _{i=l}^{j-1}\frac{\gamma_0d_i}{c}}\cos \left(\frac{\omega\sum _{i=l}^{j-1} d_i}{c}\right)\right).}
\label{FomegaN}
\end{aligned}$$

As the zeros of this function yielding the stable configurations are determined by a transcendent equation they can not be calculated analytically. Nevertheless a numerical solution requires very little effort. So in figure  we plot the closed configuration of stationary distances for six particles as a function of the bandwidth $\gamma_0$. Interestingly the particles are only equidistantly ordered in the two extreme limits of zero bandwidth $\gamma_0=0$ and a very large bandwidth, where interactions are short distance. For finite $\gamma_0$ the distances between the inner particles are larger than the distances between the outer ones.

For large bandwidth $\gamma_0$ the particles only significantly interact with their direct neighbors yielding equidistant order at a distance $d\approx\frac{3}{4}~\lambda_0$. This is the same distance as for two particles as a consequence of the short range of the force. Note that by controlling the input bandwidth we thus can switch from infinite range to nearest neighbor interaction.
This behavior is also shown in figure , where the dependence of the stable distance as function of the number of particles is compared for different bandwidth $\gamma_0$. The distance between the particles for larger $\gamma_0$ is closer to the two-particle-distance $\frac{3}{4}~\lambda_0$. Additionally they again confirm that the outer particles are closer than the inner ones and that the distances for the outer ones become more and more equal for large particle numbers.

Again the couple atom-field dynamics shows intriguing physics. Figure  and show two examples for the time evolution of four particles in a broadband field. Again the distance between the outer particles is larger than between the inner particles. Interestingly the particles tend to create high field intensities between them, with the outer two providing trap sites for the inner pair and a total field maximum generated at the center. This effect, more prominent for smaller bandwidth, reminds of a self assembled cavity QED system.

# Interparticle forces induced by blackbody radiation

As a final and extreme limit of broadband radiation let us consider the particles to be illuminated by a thermal radiation field (BBR) with a Planck spectrum. At temperature $T$ with Boltzmann constant $k_B$, Planck constant $\hbar$ and light velocity $c$ the spectral distribution reads: $$I_0(\omega)=\frac{\hbar}{2\pi^2 c^2}\frac{\omega^3}{e^{\frac{\hbar\omega}{k_B T}}-1}.
\label{bbint}$$ Of course the assumption of constant linear polarizability is only approximately valid throughout the whole Planck spectrum. Nevertheless it has been shown recently, that for atoms with a first excited state in the UV region as e.g. atomic hydrogen, this approximation is surprisingly good up to temperatures of $T\approx 6000~K$ . Introducing a more realistic frequency dependent polarizability and resonances will introduce some quantitative changes, but one can expect qualitative agreement as long as the majority of the incident power is sufficiently detuned from resonances. This is generally the case for *low temperature* radiation with a peak below the visible.

It has been shown before that the emitted BBR field from a hot source itself induces a surprising attractive force pulling atoms towards the hot object . Here we start with two temperature less particles illuminated from the side by a BBR field. Again we integrate the fiber mediated scattering force between them over the BBR frequency range in , which still can be done analytically to obtain $$\label{bbr2}
F_1=-F_2=\frac{\hbar\omega_T^4}{2\pi^2 c^3}\left(\frac{\cosh\left(\frac{2d}{r_T}\right)+2}{\sinh^4\left(\frac{d}{r_T}\right)}-3\left(\frac{r_T}{d}\right)^4\right),$$

where we defined a thermal length $r_T=\frac{\hbar c}{\pi k_B T}$ and a thermal frequency $\omega_T=\frac{\pi k_B T}{\hbar}$. For particles in an environment with room temperature the thermal length gives $r_{300K}\approx \frac{7.3\cdot 10^{-4}}{T}~m\cdot K\approx 2.4\cdot 10^{-6}~m$, which is much larger then a typical optical wavelength. Of course the long wavelength components will not be significantly guided by the fiber, but due to their low energy and momentum content, these parts are not relevant for the resulting force anyway.

From  we can expect a distance and temperature dependent sign change of this force including a zero force distance. In fact as depicted in figure , the particles will attract each other if they are close but very surprisingly this changes to repulsion at distances larger than the thermal radius $r_T$. We thus find no stable equilibrium distance. The strongest force appears at short distances and will induce a capture range below which two particles will collapse together. This will large determine the effective scattering size of two particles. As said, the zero force equilibrium point ($d_0=1.37~r_T$) is not stable and a small perturbation causes collapse or separation. This behavior is plotted in figure  and , where in the first case the particles are attracted by the high intensity peak forming between them, while they are repelled in the second case.

While the quantitative accuracy of this prediction can be doubted such a qualitative behavior should be generally true, as 3D calculations yield similar attractive behavior . This could imply important changes in the dynamics of ultra cold particle ensembles in thermal radiation fields with even astrophysical consequences.

Interestingly, the physics changes when we include backscattering and absorption of the radiation by the particles via introducing a finite $\zeta$-parameter. Figure  and  show that for real $\zeta>0$ we get long range attraction and can even find a second zero point of the force at a large distance. As this effect is very tiny it might be more of academic than practical interest, but this separation is even stable against perturbations. Including absorption via imaginary parts of $\zeta$ has the opposite effect on the particles as expected. The force is modified in a way that in the limit $\zeta_r\ll 1$ no stable configurations can be found due to long range radiation pressure contributions.

Similar to the previous case the particle also filter the light and sculpt an oscillating transmission pattern into the outgoing light. Figure  shows the outgoing intensities for a blackbody radiation field, which keep the gross shape of BBR but again the particles tend to filter some frequencies.

## Larger ensembles in BBR fields

BBR induced forces on larger ensembles can be calculated in the same way as before. As expected this leads to the rather complex analytic expression in terms of simple sums for the forces that govern the dynamics of the cloud, but allows simple numerical evaluation: $$\begin{aligned}
\fl
\eqalign{F_\Omega&=\frac{\hbar}{2\pi^2c^3}\int_0^\infty d\omega \frac{\omega^3}{e^{\frac{\hbar\omega}{k_B T}}-1}\left(\sum _{l=j}^{N-1} \cos \left(k \sum _{i=j}^l d_i\right)-\sum _{l=1}^{j-1} \cos \left(k\sum _{i=l}^{j-1} d_i\right)\right)\\
&=\frac{\hbar\omega_T^4}{2\pi^2c^3}\left(\sum_{l=j}^{N-1}\left(\frac{\cosh \left(\frac{2\sum_{i=j}^l d_i}{r_T}\right)+2}{\sinh^4\left(\frac{\sum_{i=j}^l d_i}{r_T}\right)}-\frac{3 r_T^4}{\left(\sum_{i=j}^l d_i\right)^4}\right)\right.\\
&\left.-\sum_{l=1}^{j-1}\left(\frac{\cosh \left(\frac{2 \sum_{i=l}^{j-1}d_i}{r_T}\right)+2}{\sinh^4\left(\frac{\sum_{i=l}^{j-1}d_i}{r_T}\right)}-\frac{3r_T^4}{\left(\sum_{i=l}^{j-1}d_i\right)^4}\right)\right).}
\label{Fnbbeq}
\end{aligned}$$ As it is rather hopeless to analytically find the zeros of this function, in order to get some first insight in figure  we simply plot the force on the outermost particle as function of the distance for an equidistant particle array and weak backscattering amplitude. Note that the magnitude of the near field attraction strongly increases with the particle number, while its range decreases at the same time. Hence this would trigger a fast collapse of a particle cloud once a certain small distance (high density) has been reached. At larger distances the force on the first particle does not depend on the total particle number and is repulsive. This effect is also visible in figure .

Figure , which shows the force on the first five of ten particles clearly shows that the force on the inner particles is smaller than on the outer ones. Only for some smaller distances this is not true, which is also shown in figure , while we choose a larger distance in figure . The simulations confirm that the particles do not form stable configurations but tend either to collapse together or repulse each other towards infinity. This effect can also be observed for four particles in figure  and . For this case large intensities between the particles push them apart. Here definitely more extensive and detailed simulations would be needed to fully understand the dynamics, which, however, is beyond the scope of this work.

# Conclusions

We have shown that fiber mediated broadband light scattering from linear arrays of polarizable particles induces a much more rich and complex inter particle forces than simple repulsive radiation pressure. In fact one can tailor the interaction from infinite range to nearest neighbor coupling by a suitable choice of the bandwidth. This originates from the frequency and distance dependent interference of the various frequency components in the field, which appears for dispersive as well as for absorptive scattering. For a wide range of bandwidths, one finds multi-stable regular ordering of the particles and long range forces as for optical binding in monochromatic fields. The range of inter particle interactions can be well tuned via the bandwidth of the incoming field. While we have only studied single peaked frequency distributions to limit interaction ranges, multiple peaks can be expected to allow for a much more general design of controlling interaction strength with distance. The physical effect of these generally quite weak interactions can be significantly enhanced and studied in a controlled way by coupling via optical nanostructures as nanofibers. Some preliminary calculations, nevertheless, show that qualitatively analogous behavior survives in 2D and 3D situations.

Surprisingly even for very broadband BBR fields one finds multi particle collective effects and a complex interplay of attraction and repulsion at distances around the thermal radius, which will lead to a complex nonlinear dynamics and response of larger ensembles. While this certainly is hard to observe in standard setups  it still could have significant influence on longer time scales. Alternatively it could be measurable in very force sensitive setups such as atom interferometers where it should induce extra density dependent shifts and dephasing.

We thank Juan José Sáenz, Holger Muller and Tobias Griesser for stimulating discussions. We acknowledge support via the Austrian Science Fund grant SFB F4013 .

# References


abstract: I discuss some observational properties of aspherical nebulae around massive stars, and conclusions inferred for how they may have formed. Whether or not these ideas are applicable to the shaping of planetary nebulae is uncertain, but the observed similarities between some PNe and bipolar nebulae around massive stars is compelling. In the well-observed case of Eta Carinae, several lines of observational evidence point to a scenario where the shape of its bipolar nebula resulted from an intrinsically bipolar explosive ejection event rather than an interacting winds scenario occurring after ejection from teh star. A similar conclusion has been inferred for some planetary nebulae. I also briefly mention bipolar nebulae around some other massive stars, such as the progenitor of SN 1987A and related blue supergiants.
author: Nathan Smith
title: Eta Carinae and Nebulae Around Massive Stars: Similarities to Planetary Nebulae?

# Introduction: Massive Stars and PNe

Although this is a meeting on aspherical planetary nebulae (PNe), which are the descendants of intermediate/low-mass stars, I'd briefly like to shift gears and discuss massive stars. In the nebulae around massive stars, like PNe, we see a wide variety of non-spherical geometries with a common theme of bipolar shapes in the ejecta. For most of the more stunning examples of massive star nebulae, one can usually find a PN with nearly identical appearance, at least superficially. Some of the more familiar comparisons are $\eta$ Car to Hb 5, Mz 3, or even the Red Rectangle (Soker 2007), as well as the similar multiple rings seen around SN 1987A, the luminous blue variable (LBV) star HD 168625 (Smith 2007), the Red Square (Tuthill & Lloyd 2007), and He 2-104 (Corradi et al. 2001), or the peculiar double rings around RY Scuti (Smith et al. 2002) and Abel 14, to name a few. Not surprisingly, discussions of the shaping mechanisms for massive stars and PNe share common themes: binaries/mergers vs. interacting winds vs. rotating ejections. For massive stars, though, magnetic fields still seem to be mostly taboo for the time-being.

Also, like lower mass stars during the AGB phase, it seems to be the case that massive stars shed most of their mass in a brief post-MS evolutionary phase, either as a RSG or an LBV. This was not always thought to be the case: for very massive stars, the relative importance of LBV eruptions vs. steady winds has been appreciated fairly recently because of the revised lower mass-loss rates estimated for O stars on the main sequence, and because of the very high masses of LBV nebulae (see Smith & Owocki 2006).

Despite vastly different amounts of mass and energy, the compelling similarities between massive star nebulae and PNe make it worthwhile to ask if conclusions gleaned from massive stars can inform the shaping mechanisms of PNe, and vice versa. In the interest of being provocative, then, I'll mention some results for Eta Carinae and a few other massive stars that have been studied in detail, which challenge some familiar ideas developed from the study of PNe. But first, some general comments on winds.

# Interacting Winds Scenarios?

A fast wind sweeping into a slower and denser wind is a natural avenue to pursue for shaping nebulae, and such models have had varying degrees of success in reproducing PNe shapes (there are dozens of potential references to cite here, including many in these proceedings). A similar process may occur in *some* massive stars if they pass through a very slow-wind phase as a RSG and then evolve through a faster wind phase as a BSG/LBV or a Wolf-Rayet (WR) star. This can and does produce a wind-blown bubble around the WR star in some cases.

A key point, though, is that in the case of massive stars we have some problems if we want interacting winds to account for most bipolar nebulae. First, the slower nebulae around RSGs and yellow hypergiants are generally not axisymmetric. They are often asymmetric or chaotic, but they almost never have clear signs of organized axisymetry (see, e.g., VY CMa \[Smith et al. 2001\]; NML Cyg \[Schuster et al. 2006\]; IRC+10420 \[Humphreys et al. 2002\]). Second, the resulting wind-blown bubbles around WR stars are *NOT* bipolar![^1] This lack of axisymmetry occurs despite the apparent fact that massive stars have high binary fractions. What does that mean? If binary mergers and jets blown by binaries are dominant shaping mechanisms, shouldn't we see signs of axisymmetry at all stages? Why is it the case that the only bipolar/pinched-waist nebulae around massive stars are those seen around blue supergiants such as LBVs and B\[e\] supergiants?

A critical point, I think, is that *in these blue supergiants, their escape speeds, observed ejecta/wind speeds, and surface rotation speeds are all comparable*. They are around 100–200 km s$^{-1}$, as opposed to 10–20 km s$^{-1}$ for RSGs and 1000-2000 km s$^{-1}$ for WR stars. I suspect that this is an important clue that for many nebulae around massive stars, intrinsically aspherical ejection from the surface of a rotating star is a prime agent in shaping their nebulae (see Smith & Townsend 2007). In the context of interacting winds, then, I suspect that a very interesting avenue to pursue is an aspherical fast wind interacting with a slow spherical wind or thin shell. In fact, due to the sporadic nature of episodic mass loss from massive stars, thin shells rather than steady winds is probably where most of the circumstellar mass resides.

Now, that discussion of interacting winds was for massive stars that go through a slow-wind RSG phase followed by a fast wind phase in their evolution...but that only occurs up to initial masses of about 40 M$_{\odot}$. Stars with higher initial mass (like Eta Car and most LBVs) never pass through a RSG phase. So for these stars, *ejection as an LBV is the slowest the wind speeds ever get...but it is precicesly those LBV nebulae that are observed to be bipolar*. How can this be? This means that they can't be shaped by interacting winds, because a slow dense wind blowing into a faster rarefied wind doesn't produce much interaction. Instead, the slow dense wind that follows the fast wind needs to be shaped on its own; this is discussed and amplified below. Despite the lack of interaction between the fast and slow wind, they produce shapes very similar to some PNe — a fact worth considering.

# Eta Carinae

Eta Carinae is a key object for trying to understand the shaping of bipolar nebulae, partly because is it bright and so well-observed, and partly because we have caught it so soon (only 160 yr) after its violent mass ejection, before its shape has been corrupted by interaction with the ISM. Despite its status as the most luminous and most massive star known, there is considerable overlap with some topics in PNe research, as I will highlight here.

## Energy and Momentum

Studies of the mass, kinematics, and detailed structure have led to the following basic results (summarized from Smith 2006; Smith et al. 2003):

1\. The nebula follows a Hubble-like expansion law, with the same age for the equatorial and polar ejecta.

2\. The walls of the nebula are very thin, indicating that the duration of mass ejection was less than 10% of the time elapsed since ejection.

3\. Essentially all the mass is in the thin molecular shell, formed from material ejected by the star in the outburst, not in swept-up material.

4\. The large mass, momentum, and kinetic energy came from a single exlposive event, and could not have been driven by radiation pressure alone or by the stellar wind that has blown after the eruption.

All these clues point to a single violent bipolar explosion that ejected the nebula seen today. Interestingly, all these same basic conclusions were inferred by Alcolea et al. (2007) from a similar detailed study of the PN M 1-92.

In the case of Eta Carinae, though, the difficulties for an interacting winds scenario are compounded further. The mass as a function of latitude has been measured in the bipolar lobes around Eta Car, showing that most of the mass comes from high latitudes near the pole (Smith 2006). This rules-out the familiar type of interacting-wind scenario where a spherical wind plows into a disk or torus (e.g., Frank et al. 1995), because in that scenario, the pinched waist is essentially the result of mass loading at low latitudes (note that the bipolar nebula formed in a merger model is a variation of this).[^2] Similarly, a different type of interacting winds scenario where a fast aspherical wind[^3] plows into a slower wind doesn't work either (Frank et al. 1998; Gonzalez et al. 2004). This is because we can observe the stellar wind that has been blowing after the 19th century eruption, potentially inflating and shaping the nebula. However, it is about 1000 times too weak to shape the polar lobes (like a light breeze blowing on a brick wall), and besides, the post-outburst wind speed is almost the same as that of the nebula, so the winds are not interacting anyway! There seems to be little way to escape the conclusion that the bipolar shape of the Homunculus nebula around Eta Car resulted from an intrinsically bipolar ejection by the star itself, and not from any sort of interacting winds scenario. A possible avenue to pursue is discussed after the next section.

## Double-Shell Structure

I'd like to diverge for a moment to talk about the detailed ionization structure in the walls of the nebula around Eta Car, as opposed to its overal bipolar shape, where additional similarities to some PNe can be seen. Eta Car's nebula has a distinct double-shell structure, with a thin outer shell composed of molecular gas and cool dust, and a thicker inner shell of partially-ionized atomic gas and warmer dust (Smith 2006; Smith et al. 2003). In high-resolution spectra and images of H$_2$ 2.122 $\mu$m and \[Fe ii\] 1.644 $\mu$m, this structure is almost identical to that seen in some PNe, most notably in M 2-9 (Hora & Latter 1994; Smith et al. 2005). These near-IR H$_2$ and \[Fe ii\] emission lines are usually taken as signposts for shock excitation (Shull & Hollenbach 1978) because they are seen in supernova remnants, and the double-shell structure is reminiscent of a forward/reverse shock structure that one might expect for interacting winds (e.g., Chevalier 1982).

However, they can also arise from dense atomic and molecular gas that is heated radiatively in a dense PDR (e.g., Sternberg & Dalgarno 1989). Using CLOUDY simulations, Smith & Ferland (2007) demonstrated that the observed IR emission tracers, ionization structure, and the observed dust temperatures can arise naturally from radiative heating if the two shells contain roughly the amount of mass inferred from studies of the dust (Smith et al. 2003). In fact, in the case of Eta Car, radiative heating dominates the energy budget compared to shock heating. This is comforting, because as noted earlier, the post-eruption wind speed is very similar to that of the nebula ejected in the eruption, so there is little reason to expect a strong shock anyway. So, in Eta Car, it seems clear that the observed ionization structure arises from radiative excitation, not shocks. If this is not true for M 2-9, then the almost identical ionization structure is quite a coincidence, especially since the spectra of the central objects are so similar as well (Balick 1989).

## How to Get Bipolar Lobes and a Disk

Observations of the bipolar nebula around Eta Car seem to dictate that it did not arise as a result of an interacting winds scenario, but instead, from an intrinsically bipolar wind or explosion. In other words, gas was launched from the surface of the star imprinted with the basic bipolar shape seen today.

The present-day, post-eruption wind of Eta Car is also bipolar in shape with a speed comparable to that of the nebula (Smith et al. 2003), although it is much weaker than the mass-loss rate during the 19th century eruption that made the nebula. Its almost as if the present-day bipolar wind density was simply "cranked-up" by a factor of 1000 during the outburst, maintaining the same basic speed and shape (e.g., Dwarkadas & Owocki 2002). If the bipolar nebula was created by some other external mechanism (such as jets blown by accretion onto a companion; see Soker, these proceedings) then it is a remarkable coincidence that the wind shape and speed so closely match the present-day properties of the primary star's wind. On the other hand, if the primary star ejected that material, then it is not such a coincidence at all.

In the case of an intrinsically bipolar ejection, what determines the resulting shape of a nebula is not the *density* as a function of latitude, but the *speed* as a function of latitude. When mass is driven off the surface of a star, it is a general property that this material leaves at nearly the star's surface escape speed. This is why RSGs with large radii have slow winds, and compact WR stars have very fast winds. Now, if a star is rotating fast enough to significantly modify the effective gravity at the equator (i.e. a non-negligible fraction of the critical rotation velocity), then the star's escape speed will vary with latitude, being faster at the poles and slower at the equator. Because of this simple effect, the *default* shape we should expect for material driven from the surface of a rotating star is a bipolar nebula. Admittedly, for this effect to shape the wind, the rotation speeds should be comparable to the wind speed, but this is indeed the case for blue supergiants, as noted earlier. The degree to which the waist is pinched depends on how close the star is to critical rotation. In some cases, such models can also make an equatorial disk like that seen in Eta Car.

Smith & Townsend (2007) described this type of model in detail, with test particles launched from the surface of a rotating star, following simple ballistic trajectories thereafter. They showed that it could account for the shape and speed of the polar lobes of Eta Car, as well as the basic properties of its peculiar equatorial disk. Again, this is the simplest, default shape one should expect for ejection from a rotating object. Smith & Townsend (2007) noted that it may have application to some PNe as well. Please see that paper for further details. Matt & Balick (2004) present a somewhat different intrinsic shaping model for the present-day stellar wind and disk involving MHD effects; this type of mechanism might also be relevant if the magnetic field was strong enough during the outburst.

# SN 1987A, Rings, and Mergers

The triple ring system observed around SN 1987A inspired a great deal of theoretical work on interacting winds as a potential explanation for its equatorial ring and bipolar ejecta, as well as the formation of bipolar nebulae in general. The basic favored picture is that the star had a blue loop, with a fast BSG wind pushing into a slower RSG wind. The bipolar shape and equatorial ring could arise if that RSG wind had denser material near the equator (Blondin & Lundqvist 1993; Martin & Arnett 1995), but in order for that to happen, the RSG needed to have an extra source of angular momentum, such as a binary merger event (Collins et al. 1999). This has evolved into a complex model that gives an impressive fit to the observed structure of the nebula (Morris & Podsiadlowski 2007) seen in *HST* images.

While this view is the result of considerble effort and thought, I wish to note a few flies in the ointment, which suggest that the merger model for SN 1987A might not be the final word, and may need to be revisited.

1\. A merger model followed by a transition from a RSG to BSG requires that these two events be synchronized with the supernova event itself, requiring that the best observed supernova in history happens to be a rare event.

2\. After the RSG swallowed a companion star and then contracted to become a BSG, it should have been rotating at its critical breakup velocity. Even though pre-explosion spectra (Walborn et al. 1989) do not have sufficient resolution to measure line profiles, Sk–69$^{\circ}$`<!-- -->`{=html}202 showed no evidence of rapid rotation (e.g., like a B\[e\] star spectrum). Instead, Sk–69$^{\circ}$`<!-- -->`{=html}202 had the spectrum of an entirely normal B3 supergiant.

3\. Particularly troublesome is that this merger and RSG/BSG transition would need to occur twice. From an analysis of light echoes for up to 16 yr after the supernova, Sugerman et al. (2005) have identified a much larger bipolar nebula with the *same axis orientation* as the more famous inner triple ring nebula. If a merger and RSG/BSG transition are to blame for the bipolarity in the triple-ring nebula, then what caused it in the older one?

Now, these points may seem silly at first, they are hard to reconcile with the merger model. The last one, in particular, could even be considered to be a strong rebuke. Given that more luminous blue supergiants can eject intrinsically bipolar nebulae without resorting to RSG/BSG transitions or mergers, could SN 1987A's nebula be the result of a massive star ejection instead, like an LBV?

In addition to these problems with the specific case of SN 1987A itself, we also need to take into account the growing number of observed nebulae around massive stars with rings similar to SN 1987A. Is there any evidence that they also formed from mergers? The example most people are familiar with is Sher 25, in the massive cluster NGC 3603 in our Galaxy (see Brandner et al. 1997). It has an equatorial ring with the same physical radius as that of SN 1987A, and it also has bipolar lobes. Yet, studies of the central star in Sher 25 show that its abundances indicate that *it has not gone through a RSG phase* (Smartt et al. 2002). The newly discovered SBW1 in Carina also has a 0.2 pc radius identical to the equatorial ring around 87A; its central star has about the same luminosity as that of the progenitor of SN 1987A, but its ring has Solar N abundances, so it has also *not been through a RSG phase* (Smith et al. 2007). Finally, there's the triple-ring nebula around the LBV star HD 168625 (Smith 2007), which appears to be almost identical to that of SN 1987A. While this star could indeed be a post-RSG because it is mildly N-rich, it is an LBV – stars well known for their unstable bipolar mass ejections. Furthermore, Smith (2007) argued that if the triple ring nebula of HD 168625 is coeval (polar rings have the same age as the equator, as in SN1987A), then they could not have been ejected in a RSG phase because they would be far too fast.

So while the progenitor of SN 1987A may very well have passed through the RSG phase needed for the merger model, two of its twins did not, and a third is an LBV, where the shell was likely created in an LBV ejection. If a binary merger in the RSG phase really is required to make the ring around SN 1987A, how and why did at least three other objects make *nearly identical* nebulae in a different way? At the very least, this implies that a merger model is not the only viable option. Keep in mind that their central stars still exist and are easy to observe, yet there is currently no evidence that they are post-merger products. In fact, they appear as fairly normal blue supergiants; like the progenitor of SN 1987A, these other rings do not appear to be unusually rapid rotators as one would expect for a post-merger product (their inclinations are known from the equatorial ring nebulae). If further studies reveal that any of them survive today as close binaries (binaries that have not yet merged, so that the observed rings are obviously not from a merger event), it will critically wound the merger hypothesis for SN 1987A.

**Acknowledgements.** I would like to thank the conference organizers for the invitiation to speak and for generous partial financial support that made my trip possible.

[^1]: Some WR nebulae are a little egg-shaped, but we don't see any with pinched waists. Also, I'm not including the "pinwheels" and related phenomena (see, e.g., Tuthill et al. 1999) around dust-producing WC stars, which represent a very different type of "interacting winds".

[^2]: By the way, note as well that the disk seen in *HST* images of Eta Car is *not* the agent responsible for pinching the waist of the bipolar nebula, because it is the same age or younger.

[^3]: The present-day wind appear to be bipolar on size scales smaller than the binary separation (see Smith et al. 2003; van Boekel et al. 2003).


abstract:
author: Peter Flach; Kacper Sokol
bibliography: paper.bib
title:  Simply Logical – Intelligent Reasoning by Example
  (Fully Interactive Online Edition)

# Evolution of The Book

"Simply Logical – Intelligent Reasoning by Example" by Peter Flach was first published by John Wiley in 1994 . It could be purchased as book-only or with a 3.5 inch diskette containing the SWI-Prolog  programmes printed in the book (for various operating systems). In 2007 the copyright reverted back to the author at which point the book and programmes were made freely available online; the print version is no longer distributed through John Wiley publishers. In 2015, as a pilot, we ported most of the original book into an online, interactive website using SWI-Prolog's SWISH platform . Since then, we launched the *Simply Logical* open source organisation[^1] committed to maintaining a suite of freely available interactive online educational resources[^2] about Artificial Intelligence and Logic Programming with Prolog. With the advent of new educational technologies we were inspired to rebuild the book from the ground up using the Jupyter Book platform  enhanced with a collection of bespoke plugins that implement, among other things, interactive SWI-Prolog code blocks that can be executed directly in a web browser. This new version is more modular, easier to maintain, and can be split into custom teaching modules, in addition to being modern-looking, visually appealing, and compatible with a range of (mobile) devices of varying screen sizes.

# Artificial Intelligence and Logic Programming with Prolog

The book discusses methods to implement intelligent reasoning by means of Prolog programmes. It is written from the shared viewpoints of *Computational Logic* – which strives to automate various kinds of reasoning – and *Artificial Intelligence* – which seeks to implement aspects of intelligent behaviour as computation. The combination of these two perspectives distinguishes this book from among its peers and offers a unique learning experience. The readers working in Artificial Intelligence will find a detailed treatment of how the power of logic can be harnessed to solve some of the (practical) problems they may be facing. To support this learning objective the book provides a variety of interactive code examples in the domains of natural language interpretation, abductive and inductive reasoning, and reasoning by default. Those acquainted with Logic Programming, on the other hand, will be interested in the practical side of many topics that in other educational materials mostly receive theoretical treatment – many advanced programmes presented and explained in this book are not, in a didactic form, available elsewhere. The readers unfamiliar with either field will benefit from a comprehensive learning resource that collects diverse subjects and presents them consistently across accessible and engaging modules. This approach offers the best of both worlds: it introduces the theory of Logic Programming that does not intimidate novices to the field with superfluous mathematical machinery, nonetheless it offers strong theoretical foundations for the practical aspects of programming for Artificial Intelligence.

The book consists of three parts. Part I deals with Logic and Logic Programming and covers:

-   key concepts in Logic Programming, such as programme clauses, query answering, proof trees, and recursive data structures (Chapter 1);

-   a rigorous discussion of resolution theorem proving in clausal logic, stepping through propositional clausal logic, relational clausal logic (without functors), full clausal logic, and definite clause logic, on the way dealing with concepts such as Herbrand models and resolution refutations, as well as meta-theoretical notions like soundness and completeness (Chapter 2); and

-   practical aspects of Prolog programming, including SLD-trees, cuts, arithmetic expressions, second-order predicates (e.g., `setof`), various programming techniques (e.g., accumulators and difference lists) as well as a general programming methodology and meta-interpreters (Chapter 3).

Part II shifts the focus from the Logic Programming perspective to the Artificial Intelligence viewpoint, predominantly dealing with graphs and search. Chapter 4 discusses graphs found naturally in Prolog: in the form of trees represented by terms (e.g., parse trees) and as search spaces spanned by predicates (e.g., SLD-trees). Next, Chapter 5 overviews depth-first search, iterative deepening, and breadth-first search in the context of Logic Programming; moreover, it develops a breadth-first Prolog meta-interpreter and an (inefficient) interpreter for full clausal logic, in addition to illustrating forward chaining with a programme that generates Herbrand models of a set of clauses. Chapter 6 discusses best-first search and its optimality, leading to the A$^\star$ algorithm and a brief discussion of non-exhaustive heuristic search strategies (beam search and hill-climbing).

Finally, Part III encompasses advanced topics. Chapter 7 covers natural language parsing and interpretation, in particular definite clause grammars. It also looks at semantics and natural language generation, culminating in a simple question-answering agent. Chapter 8 discusses reasoning with incomplete information in view of the Closed World Assumption and Predicate Completion. Specifically, it introduces abductive reasoning as well as default reasoning by means of negation as failure and defeasible default rules. Then, Chapter 9 deals with inductively inferring a logic programme from examples, introducing and implementing concepts such as generality between clauses and anti-unification, which are fundamental to Inductive Logic Programming.

Appendices offer answers to selected exercises, an overview of built-in predicates in Prolog, a library of utility predicates underpinning many programmes introduced throughout the book, and listings of two larger programmes for transforming a Predicate Logic formula to clausal logic and performing Predicate Completion.

As suggested by the title, this book presents intelligent reasoning techniques by example, therefore every method is accompanied by a Prolog implementation that, in the online edition, can be executed directly in a web browser. These code listings serve two didactic purposes: they engender understanding through hands-on experience, but more importantly the declarative reading of each implementation constitutes an integral part of the explanation of the underlying technique. The more elaborate programmes are carefully designed to explicitly communicate the steps taken to achieve the goal, explaining key issues along the way. The book therefore does not just focus on "How is this done in Prolog?", but rather on "How should I solve this problem, were I to start from scratch?" In other words, it embodies the "teaching by showing, learning by doing" philosophy. Despite the strong focus on practical examples, a substantial portion of the book is devoted to the theoretical underpinnings of clausal logic and Logic Programming, however their presentation is primarily motivated by their instrumental role in understanding and implementing each particular technique. [^3]

The book has been the foundation of multiple courses taught to undergraduate and master's students at the University of Bristol and elsewhere, e.g., King's College London, Queen Mary University of London, University of Cardiff, University College Cork, Tufts University, St. Joseph's University, University of Freiburg, Uppsala University, Free University Brussels, and University of Sevilla[^4]. Specifically, the *Artificial Intelligence with Logic Programming* undergraduate course[^5] offered at the University of Bristol follows the book quite closely; a collection of slides whose structure mirrors the presentation of material in the book is freely available as a supplementary educational resource[^6]. A spiritual successor of this unit is *Computational Logic for Artificial Intelligence* now taught as part of the curriculum delivered by the Interactive Artificial Intelligence Centre for Doctoral Training[^7] (IAI CDT). This unit takes a more modern approach by discussing the interconnections between logic programming (with SWI-Prolog) and recent developments in data-driven machine learning[^8].

# The Online Edition

This online book offers an intuitive and comprehensive introduction to *Artificial Intelligence and Logic Programming* in the form of an accessible educational resource with a modern look and feel. The learning experience is enhanced with *exercises* and *practical coding examples*; the latter are delivered through interactive code boxes that can be executed directly in a web browser, thus fulfilling the original vision for the paperback edition that was curtailed by the lack of appropriate technologies at the time of its publication. It is written with three kinds of readers in mind: *Artificial Intelligence* researchers or practitioners; *Logic Programming* researchers or practitioners; and students (advanced undergraduate or graduate level) in both fields. The book can be used as a teaching aid, but it is also suitable for self-study. It can be adopted as is or the content can be rearranged into a bespoke learning resource, which is made possible by the modular design of the book source that is split into a collection of files based on their type:

-   text is written in Markdown – each sub-section is placed in a separate file;

-   code is released as SWI-Prolog scripts; and

-   figures are available as Scalable Vector Graphics (SVG).

All of these are published on GitHub[^9], which allows to easily reuse these materials, incorporate them into bespoke courses, or adapt them into alternative educational resources such as practical training sessions.

Having used the book for undergraduate and postgraduate teaching for many years, we decided against making major structural changes or adding substantial new material. We have, of course, taken the opportunity to include small improvements and corrections, such as using anonymous variables wherever appropriate, and slight punctuation improvements in clauses and queries. Other enhancements include additional examples collected during years of teaching, formatted in clearly identifiable boxes.

From a technical perspective, the development of the online edition required us to implement a collection of Jupyter Book plugins (that also work with Sphinx) spanning functionality specific to SWI-Prolog[^10] , later extended to be compatible with the cplint[^11]  and ProbLog[^12]  programming languages. Specifically, we released `sphinx-prolog`[^13], which allows to embed interactive SWI-Prolog and cplint code boxes – see Figure  – by including their source directly in a Markdown file via a custom Prolog code listing or load the programme from an external code file. This plugin is based on SWISH[^14] – SWI-Prolog for Sharing – which is an online, interactive SWI-Prolog coding environment akin to Jupyter Notebooks . The `sphinx-problog`[^15] extension, on the other hand, allows to embed interactive ProbLog code boxes, the implementation of which is based on the online execution environment underpinning the code examples published as part of the ProbLog website. Additionally, building these plugins prompted us to explore alternative technologies for composing (interactive) training resources, which in turn has inspired a prototype of a new publishing workflow in which multiple artefacts such as online documents, slides, and computational notebooks can be generated from a unified collection of source materials .

[^1]: <https://simply-logical.space/>

[^2]: <https://github.com/simply-logical/>

[^3]: The book's preface offers a more comprehensive explanation of its instructional design and content: <https://book.simply-logical.space/src/simply-logical.html#author-s-preface>.

[^4]: <http://people.cs.bris.ac.uk/~flach/SL/Courses.html>

[^5]: <https://coms30106.github.io/>

[^6]: <https://github.com/simply-logical/slides/>

[^7]: <https://www.bristol.ac.uk/cdt/interactive-ai/>

[^8]: <https://github.com/simply-logical/ComputationalLogic/>

[^9]: <https://github.com/simply-logical/simply-logical/>

[^10]: <https://www.swi-prolog.org/>

[^11]: <https://cplint.ml.unife.it/>

[^12]: <https://dtai.cs.kuleuven.be/problog/>

[^13]: <https://github.com/simply-logical/sphinx-prolog/>

[^14]: <https://swish.swi-prolog.org/>

[^15]: <https://github.com/simply-logical/sphinx-problog/>