abstract: # **Summary**
  .
  The strong CP problem is that $SU(3)$ gauge field instantons naturally induce a CP violating term in the QCD Lagrangian which is constrained by experiment to be very small for no obvious reason. We show that this problem disappears if one assumes the existence of at least one black hole somewhere in the universe. The argument is reminiscent of Dirac's argument for the quantization of charge, in which the existence of one monople anywhere in the universe suffices to require the quantization of electric charge everywhere.
author: John Swain
date: revised version June 21, 2010
title: Black Holes and the Strong CP Problem

# Introduction

In quantum chromodynamics (QCD) - the generally accepted theory of quarks and gluons - there was a prediction that there should be a light pseudoscalar particle associated with the conserved current associated with global chiral rotations of the quarks. No such meson was observed, and this was called the "U(1) problem". It was realized that the quantum effects spoiled the conservation of the quark axial current, making its divergence proportional to $TrF_{\mu\nu}F^{\ast\mu\nu}$ where $F_{\mu\nu}$ is the $SU(3)$ field strength, $F^{\ast\mu\nu}$ its Hodge dual, and the trace is taken over $SU(3)$ indices.

This divergence corresponds to a CP-violating Lagrangian density of the form

$$L_{CP-violating} = \frac{\theta g^2}{32\pi^2}TrF_{\mu\nu}F^{\ast\mu\nu}
\label{eqn:CP}$$

where $g$ is the $SU(3)$ gauge coupling constant, and $\theta$ is a free parameter. Overall, this expression is proportional to the Pontryagin density, which on integration over spacetime is an integer topological charge representing the number $n$ of times that $S^3$ (considered as physical space ${\mathbb{R}}^3$ plus a point at infinity) nontrivially "winds around" $SU(3)$. The physical gauge-invariant vacuum is constructed as a superposition of states of winding number $n$, each weighted by $e^{in\theta}$ with the sum running from $n=-\infty$ to $n=\infty$ in order to preserve gauge invariance under the "large" gauge transformations which are not continuously connected to the identity. $\theta$ is not determined by the theory, and can, in principle, take any value between $0$ and $2\pi$.

When the weak interactions and quark masses are included, $\theta$ is shifted by an amount $arg(det(m))$ where $m$ is the quark mass matrix, but the basic form of the expression remains the same and unless the shifted $\theta$ is zero (or $\pi$, but this subtlety will not concern us here) this term leads to a (CP-violating) electric dipole moment $|d_n|$ for the neutron. The present upper bound $|d_n|<2.9\times 10^{-26}e\cdot$cm, where $e$ is magnitude of the electron charge , implies $\theta <10^{-9}$. The puzzle of why $\theta$ is so small is the "strong CP problem".

A wide variety of solutions have been proposed, generally involving new physics. Many postulate particles called axions associated with an additional $U(1)$ symmetry which can be used to rotate $\theta$ to zero. These have not been observed and are in general quite constrained by astrophysical considerations. Other ideas include adding dimensions to the usual 3+1 that we know , or making them some fractional value a little less than four . Two-dimensional fundamental objects (2-branes) have been considered, as have microscopic wormholes , hypothetical new interactions , new (non-axion) particles , supersymmetry , and magnetic monopoles. It has also been claimed that certain choices of regularization techniques could solve the problem. It has even been suggested that a staggering $10^{32}$ standard model copies could do the job . It has also been argued that the strong CP problem might naturally not appear at all if one simply reformulated QCD in terms of holonomies (gauge invariant traces of Wilson loops). This list of ideas and references is not meant to be complete, but rather to show that the problem has driven theorists to a wide range of quite exotic scenarios in the search for an explanation. Despite all this creativity, the strong CP problem is still generally considered unsolved.

The point of this paper is that the problem could be resolved without unobserved exotica, and without spoiling the solution of the $U(1)$ problem, if the spacetime integral of the Pontryagin density were somehow to be zero – something I now argue will happen if one allows for the existence of even one black hole.

In elementary particle physics one usually ignores gravity, and works with quantum field theory in flat and topologically trivial spacetime. While quantum field theory in a general curved spacetime is highly nontrivial, the question asked here only requires a little topology.

First let us recall where the instantons come from that lead to the strong CP problem. We look for $SU(3)$ gauge field configurations $A_\mu$ which go to the identity (up to a gauge transformation) at spatial infinity, with all the directions at infinity identified. These turn out to fall into topologically distinct classes labelled by elements of $\pi_3(SU(3))$.

For completeness, and to make clear the origin of the $\pi_3(SU(3))$, let us repeat the argument in more detail. Pure gauge configurations are of the form $A_\mu=iU^\dag(x)\partial_\mu U(x)$ where $U(x)$ takes its values in $SU(3)$ and $x=(t,\vec{x})$. Using the gauge freedom to set $A_0=0$, which essentially means we consider time-independent fields, only partially fixes the gauge. If we require $U(\vec{x})=1$ as the spatial $\vec{x}$ goes to infinity in all directions this is the same as looking for maps from spatial ${\mathbb{R}}^3$ compactified at infinity (that is, $S^3$) into $SU(3)$. Instantons and correspond to homotopy classes $[S^3,SU(3)]$ of these maps. By definition, $[S^3,SU(3)]=\pi_3(SU(3))$ and since $\pi_3(SU(3))=Z$ we have homotopically distinct maps labelled by the integers, which turn out to be the very topological charges that come from the integration of the Pontryagin density in equation .

The key observation of this paper is that if we have black holes present and repeat the argument, we should replace $[S^3,SU(3)]$ by $[M,SU(3)]$ where $M$ is a manifold (now with boundary) created from the $S^3$ described above with a 3-ball bounded by a 2-sphere excised for each black hole present – effectively we are removing a set of distinct points (and balls around them) from space.

Physically, we require that the gauge fields go to the identity (up to gauge invariance) on the surfaces of black holes (as well as at infinity), in a similar spirit to reference in which this condition is invoked to argue for spacetime foam as a universal regulator. Note that we don't need to assume spacetime foam or wormholes (as have been used to argue for solutions to the strong CP problem before) and the black holes in question need not be virtual or microscopic – any astrophysical black holes (or, indeed just one) would do. In particular, $\pi_1(M)$ is assumed to be trivial, as is usual in considerations of the strong CP problem (and for which there is no experimental evidence to the contrary) . In many ways, this is meant to be a very conservative solution to the strong CP problem invoking essentially no new physics beyond what is generally known.

Now let us consider the homotopy classes of maps $\left[M,SU(3)\right]$ from $M$ to $SU(3)$. $M$ is clearly simply connected ($\pi_1(M)=0$), as every closed loop can be continuously shrunk to a point. If we consider possibly topologically nontrivial maps from $M$ to $SU(3)$ then the usual Postnikov construction tells one that one has to consider $\pi_2(SU(3))$, but this is zero, and one is left with nothing to worry about except $\left[M,K(\pi_1(SU(3)),1)\right]$ with $K(\pi_1(SU(3)),1)$ being the relevant Eilenberg-MacLane space.

By definition that means that $\left[M,K(\pi_1(SU(3)),1)\right]=H^1(M,\pi_1(SU(3)))$. Since $M$ is simply connected, one immediately sees that this is zero, and thus all maps from $M$ to $SU(3)$ are homotopically trivial (continuously deformable to the identity). We could argue directly that it is also zero due to the fact that $SU(3)$ is simply connected and $\pi_1(SU(3))=0$.

If one wants to argue that the true gauge group should be $SU(3)$ with its $Z_3$ center divided out, making the first homotopy group nonzero, then the first argument given in the above paragraph still makes the case.

The integral in equation now vanishes as the corresponding topological charge is zero, and the strong CP problem would seem to be solved. Clearly, analogous arguments hold for any finite-dimensional Lie group $G$ in place of $SU(3)$ since $\pi_2(G)$ is always zero in this case and the same reasoning applies. Some care is needed if multiply connected $M$ is considered since one does not want to induce a $\theta$-like term for the $U(1)$ of electromagnetism. Such an electromagnetic $\theta$ term is absent in standard analyses since $\pi_3(U(1))=0$ and thus there are no $U(1)$ instantons to worry about.

In the case of topologically more complicated spacetimes additional gravitationally-induced CP violating effects may be present. In particular, terms proportional to $f_{\mu\nu}f^{\ast\mu\nu}$ where $f$ is the electromagnetic field strength tensor and $R^{ab}_{\mu\nu} R^{\ast\mu\nu}_{ab}$ where $R$ is the spacetime curvature tensor can be present. These contributions are not usually considered part of the "strong CP problem", although it is very interesting that these terms are not obviously suppressed by powers of the Planck scale. In the case of the $R^{ab}_{\mu\nu} R^{\ast\mu\nu}_{ab}$ the spacetimes involved clearly are not of the form considered here since corresponding instantons do not refer to topologically nontrivial gauge fields *over* spacetime but rather topologically nontrivial spacetimes. Whether the arguments made here can be extended to this case is not obvious, but I hope to be able to return to this interersting question in future work. Of related interest is also in which the suggestion is made that the usual instanton sums may need to be modified in some theories.

As this paper was being completed, I became aware of a related paper by Etesi. This author considers both black and "white" holes (which it is not clear exist), finding results for $\left[M,SU(3)\right]$ which agree with those here. The claim in that paper however is not that the Pontryagin term integrates to zero, but rather that one should consider a sort of "effective homotopy" which takes into account the causal structure of the relevant spacetime and for which the corresponding homotopy classes are not trivial and the strong CP problem remains. The idea is that one should only consider homotopies whose initial and final stages can be compared by an observer in finite time. This leads to a re-appearance of the $\theta$ vacuum structure which we just got rid of, and the solution of the strong CP problem is based an additional assumption which is certainly not required in the usual formulation of the strong CP problem. In fact $\theta$ arises in a quantum mechanical superposition of states of all instanton numbers making even the meaning of a suitable observer unclear at best. Indeed the term "instanton" refers to the fact that one considers field configurations which can be thought of as at least approximately localized in time. This leads to that paper missing the key point I make here which is that even a single black hole (no "white holes" needed) suffices to make all the $SU(3)$ field configurations topologically trivial. In this way $TrF_{\mu\nu}F^{\ast\mu\nu}$ can still be nonzero locally to solve the $U(1)$ problem, while globally all the corresponding gauge field configurations are topologically trivial.

In contrast to essentially all other attempts to solve the strong CP problem, the approach presented here requires no modification of the standard treatment of the problem other than to include the presence of normal (indeed classical) black holes as part of the structure of spacetime. No undiscovered exotica need be invoked.

It may seem surprising that the existence of even one singular object - in this case a black hole - could have implications for elementary particle physics, but there is actually a rather old analog. Long ago in 1948, Dirac had used topological arguments to show that the presence of just *one* magnetic monopole would require electric charge everywhere to be quantized. Here we see that, similarly, the presence of just one black hole can resolve the strong CP problem.

# Acknowledgements

I would like to thank Kaća Bradonjić and Tom Paul for reading early drafts of this paper. I would also like to thank Egil Lillestol, Danielle Metral, Nick Ellis, and the hospitality of the CERN Latin American School of High Energy Physics 2009 in Colombia where this work was started, Luis Alvarez-Gaumé since it was during his lecture to the students that I started to think about this problem again, and of course all the students who made the school the great success that it was! I would also like to thank Gábor Etesi for email correspondence on the first draft of this paper, as well as for pointing out that while his argument in reference allowed for white holes it does not require them, and for emphasizing the possible differences in considering Euclidean and Minkowsi spacetimes in calcuations. Thanks are also due to Michael Duff and Stanley Deser for email correspondence and bringing references and to my attention. This work was supported in part by the US National Science Foundation under grant NSF0855388.


author: B. M. Swinyard; P. Ade; J-P. Baluteau; H. Aussel ; M. J. Barlow; G. J. Bendo; D. Benielli; J. Bock; D. Brisbin; A. Conley; L. Conversi; A. Dowell; D. Dowell; M. Ferlet; T. Fulton; J. Glenn; A. Glauser; D. Griffin; M. Griffin; S. Guest; P. Imhof; K. Isaak; S. Jones; K. King; S. Leeks; L. Levenson; T. L. Lim; N. Lu; G. Makiwa; D. Naylor; H. Nguyen; S. Oliver; P. Panuzzo; A. Papageorgiou; C. Pearson; M. Pohlen; E. Polehampton; D. Pouliquen; D. Rigopoulou; S. Ronayette; H. Roussel; A. Rykala; G. Savini; B. Schulz; A. Schwartz; D. Shupe; B. Sibthorpe; S. Sidher; A. J. Smith; L. Spencer; M. Trichas; H. Triou; I. Valtchanov; R. Wesson; A. Woodcraft; C. K. Xu; M. Zemcov; L. Zhang
date: Received 2024-10-02; accepted
title: In-flight calibration of the $Herschel$-SPIRE instrument[^1]

# Introduction

The in-flight calibration of any astronomical instrument relies on a combination of previous measurements using instrumentation with known or traceable calibration accuracy, accurate models of standard astronomical sources and calibration sources internal to the instrument or facility. In this paper we discuss the methods and source models employed to convert SPIRE data to physical units and discuss the accuracy of the calibration and any caveats that must be placed on the data.

The SPIRE instrument and its overall calibration scheme are described in detail in Griffin et al. (). For the purposes of discussion of the calibration and performance estimation we treat SPIRE as two separate instruments that are spatially separated at the focal plane of the *Herschel  telescope: a three band imaging photometric camera with bands centred nominally on 250 $\mu$m (PSW), 350 $\mu$m (PMW), and 500 $\mu$m (PLW) with detector arrays of 139, 88 and 43 respectively feedhorn coupled NTD-bolometers (Turner et al ), and a two band imaging Fourier transform spectrometer (FTS) covering 194-313$\mu$m (SSW) and 303-671$\mu$m (SLW) with arrays of 37 and 19 detectors. SPIRE contains two internal calibration sources which were designed to provide a rapidly modulated signal for relative gain calibration (the PCAL source) and a stable thermal source for calibration of the spectrometer and to balance the power from the telescope (the SCAL source). PCAL is placed at an image of the telescope secondary in the common light path of the two instruments and therefore can be used to stimulate all detectors in the photometer and spectrometer. The SCAL source is placed in the second input port of the spectrometer (Swinyard et al. ). It is therefore constantly viewed during observations and provides a reference against which the spectrum from the telescope port is measured. The intention before flight was that the SCAL source would be heated to a temperature sufficient to entirely "null" the spectrum from the telescope. We discuss in Sect. why this has proved to be unnecessary for standard observing conditions. The internal calibrators are only viewed through part of the optics chain and the overall astronomical calibration and performance of the instrument and telescope together can only be established by observation of astronomical sources as we describe in this paper.*

The structure of the paper is as follows: in Sect. we describe the initial in-flight performance evaluation of the instrument before going on to cover the basic conversion of the data from digital encoder values to physical units for the detectors, the beam steering mirror (BSM) and the spectrometer mechanism (SMEC) in Sect. . In Sect. we describe the characterisation of the time response of the photometer bolometers using ionising radiation hits and in Sects. and we discuss how the photometric calibration was established for the photometer and spectrometer. In Sect. we present the results of comparing SPIRE observations to solar system objects with known or well constrained model flux densities. The beamsize and astrometric accuracy are discussed in Sect. and the wavelength accuracy in Sect. . We summarise the present calibration status of SPIRE and draw some conclusions for future calibration activities in Sect. .

# Initial performance evaluation

The initial performance verification and comparison between in-flight and ground performance was carried out using the PCAL source whilst the *Herschel cryostat lid was closed. This allowed a direct measurement of the detector response under conditions similar to those seen during ground testing. The response measured using PCAL before opening of the cryostat lid showed a close correlation between ground and flight in all detector channels in both the photometer and spectrometer. The overall in-flight measurement showed a slightly higher response in all arrays as the bolometers are approximately 10 mK colder in flight than during spacecraft level ground tests. The time response of the detector channels is a combination of the thermal response of the bolometers and the low pass filters in the readout electronics. In the photometer channels the time response is dominated by the low pass filters which have a nominal 5 Hz 3 dB cut off. Here the PCAL flashes, which have switch-on time constant of  90 msec and a switch-off time constant of  50 msec, can be used to make a first order assessment of the response of the system and look for anomalous behaviour. This test is not possible for the faster (25 Hz) spectrometer filters and bolometers although analysis of the interferograms does provide a first order check. This is achieved by measuring the difference in recorded position of the peak in signal at zero path difference when the spectrometer mechanism is scanned in its "forward" and "reverse" directions. In both the photometer and spectrometer all channels that were working nominally before flight showed the same or similar time response. A more detailed investigation of the photometer time response was carried out using an analysis of ionising radiation glitches and is discussed in Sect. .*

Once the lid was opened we were able to observe both astronomical sources with known or well constrained model flux densities and the emission from the *Herschel  telescope. The telescope primary and secondary mirror temperatures are measured using several thermistors with a notional 0.1% accuracy. This, in addition to measurements of the mirror emissivity (Fischer et al. ), allows us to model the direct emission from the telescope which fills the field of view of all detectors. By varying the electrical power applied to the bolometers through changing the bias voltage - i.e. taking a so called "loadcurve" - it is possible to make a direct estimate of the optical power absorbed by the bolometer. Comparing this measurement to the expected absorbed power using the modelled and ground based measurement of the instrument throughput and transmission showed that a) there is little or no straylight in the SPIRE beam due to sources on the satellite and b) that the instrument throughput and transmission are as good or better than expected before flight in both the photometer and spectrometer.*

To assess the optical performance of the system and to calibrate the response of the photometer to point sources two separate types of scanned observations were carried out. The first takes a scan map observation using the standard Astronomical Observation Template (AOT) of a point source of known or modelled flux density. In reconstructing the map from these data the time series from many different detectors in the array are combined to make a single image of the beam. The second approach was to scan a bright point source along the length of the array in a raster pattern with fine steps between the scan legs allowing the beam map, relative position and response of each individual bolometer to be determined. The results from both the photometer and spectrometer arrays confirmed that the beam sizes are as expected and that no significant change to the optical alignment of the instrument had occurred during launch. These measurements are discussed in further detail in Sect. .

# Data processing and calibration

The data processing pipeline for the SPIRE photometer and spectrometer data are discussed in detail in Griffin et al. () and Fulton et al. (). In this section we highlight those elements of the processing that require the application of calibration data and discuss the origin and accuracy of the applied conversion factors.

The basic conversion from digitised output from the SPIRE electronics to voltage (in the case of the detectors) or position (in the case of the BSM and SMEC) relies on the accurate measurement of the gains and offsets of the amplification chains and encoders in the SPIRE warm electronics. These were extensively tested during the ground calibration and again in flight. Fixed resistors mounted on the bolometer arrays were used to check the detector conditioning electronics gains and no change was found compared to the ground measurement. The values of the gains applied in pipeline processing can be assumed to much more accurate than the overall calibration based on astronomical flux densities and need not be considered in estimating the calibration accuracy.

The SMEC position encoder accuracy was tested by comparing the measured and expected frequencies of known lines (see Sect. ). The position encoder of the BSM has been recalibrated in flight by using it to move a bright point source across the field of view and using the known positions of the detectors (see Sect. ) to determine its position. Using this method, combined with an end to end optical model, we can determine the position of the BSM to within an equivalent of 0.5 arcsec on the sky, much better than the absolute pointing error of the *Herschel  satellite which is estimated as 2 arcsec (Pilbratt et al. ).*

# Photometer bolometer time constants

The basic time response of the photometer detector channels is dominated by the low pass filtering in the signal conditioning electronics which has been accurately modelled and verified during ground testing. However, the thermal response of the bolometers themselves has a second order effect that varies with bolometer temperature and it is important that this is understood as correcting for it is critical in both the astrometric and photometric accuracy achievable in SPIRE maps. The first assumption has been that the bolometer thermal response is represented by an 'RC' type response characterised by a single time constant. Investigation into the presence of any slower bolometer time constants on scan map and spectrometer time domain data has so far revealed no evidence for a slow response to optical power and only the correction for a first order fast time response is currently being performed by the pipeline processing.

High energy particles hitting the bolometers cause signal glitches as the energy is absorbed and the temperature of the bolometer increases. The response to these glitches gives a direct measurement of the bolometer and electronics time constants. To measure the photometer bolometer time-constant the data from a large number of photometer glitches have been collected and compared to the theoretical impulse response modelled as a delta-function convolved with the photometer electronics and bolometer transfer functions. Applying this process to sets of individual (single event) glitches, reveals that the photometer glitch response for most photometer bolometers can be best described by a two-component transfer function, with a fast time constant $\tau_1$ $\approx$ 6 ms and a second slow component with a time constant of $\tau_2$ $\approx$`<!-- -->`{=html}105 ms and amplitude about 6.6 % (Fig. 1). A second type of common, concurrent glitches are also seen which are most likely due to interaction between particles and the bolometer array substrate. Applying the same process to these reveals that the common glitch response is best described by a first-order-only bolometer transfer function with a time constant of about 6 ms. This disparity in the system response between single-event and common glitches, and how this relates to the absence of any slow response to optical power, is not yet understood and is currently being investigated. The frequency response correction of both the electronics and the bolometers amounts to no more than 5% of the signal for standard scan map data and given our knowledge of the electronics and bolometers any uncertainties in this correction can, again, be safely ignored in the overall calibration uncertainty.

# Photometer photometric calibration

The conversion from voltage to flux density in the photometer is based on the assumption that the small signal responsivity can be related to the bolometer operating point voltage nearly linearly i.e. that we can write the inverse responsivity in the form $$\frac{dS}{dV} = K_1 + \frac{K_2}{V-K_3}
\label{equation1}$$ where S is the flux density, $K_1,K_2,K_3$ are constants and $V$ is the measured voltage across the bolometer. $K_1$ is in units of Jy/V, $K_2$ is in Jy and $K_3$ is in V. Integrating this equation gives $$S = K_1(V-V_0) + K_2 \ln \frac{[V-K_3]}{[V_0-K_3]}
\label{equation_calibration}$$ where now $V_0$ is another arbitrary constant that we chose to represent the bolometer voltage measured under a standard observing condition on a dark region of the sky. The various parameters have been determined from measurements in-flight using PCAL to measure the relative responsivity while viewing sources of different brightnesses in order to change the bolometer voltage.

The calibration curves were measured for each detector in two steps. First, the change in voltage $\Delta V_{PCAL}$ observed for PCAL flashes were measured when the telescope was pointed at regions of the sky with different flux densities ranging from dark sky to the peak of Sgr A (for the nominal voltage settings) or Sgr B2 (for the lower gain bright source mode settings). The plot of $1/\Delta V_{PCAL}$ versus $V$ gives the unscaled version of Eq. ; an example is shown in Fig.  where the telescope was pointed at the dark sky and different locations in the region of SgrA to vary the flux onto the detectors. The calibration curve itself is then scaled using data from observations in which Neptune (or a source calibrated using Neptune) is scanned across each detector in each array. The feedhorns and bolometers are not inherently sensitive to polarisation and any instrinsic polarisation in the calibration targets will have negligible effect on the overall photometric calibration.

The initial absolute calibration scale was set using the asteroid Ceres and a model of its submillimetre flux density based on the Standard Thermal Model (Lebofsky et al ). We have since observed Neptune and correction factors between the output of the photometer pipeline and the model of Neptune based on the work of Moreno (, ) have been established (see Griffin et al. ). We discuss the accuracy of the Moreno model further in Sect. and the accuracy of the calibration constants in Sect. of the present paper.

# Spectrometer pipeline photometric calibration

The calibration of the spectrometer follows a different method compared to the photometer. The signal that is measured by the spectrometer detectors is not a direct measurement of the flux density integrated over the passband as in the photometer but rather the Fourier component of the spectral content. Therefore Eq.   is not directly applicable but an analogous equation can be used to correct for any non-linearity between absorbed power and bolometer voltage before transforming into the frequency domain. The parameters in this scheme are derived from a model response of the bolometer and have different values depending on the bolometer bias that is set.

Once a linearised timeline in volts has been obtained the signal versus optical path difference is calculated using the mechanism position and further corrections for phase error are made (see Fulton et al. ). The source, telescope, SCAL and the instrument self emission are always measured together and it is necessary to subtract a reference interferogram taken on dark sky to obtain the source on its own. This difference interferogram is transformed into spectral space and finally converted to flux density using a relative spectral response function (RSRF). The RSRF is derived by taking the interferogram of an astronomical source with a well modelled continuum, subtracting the reference interferogram, transforming this into spectral space and dividing by a model of the source spectrum. In the present pipeline we use the asteroid Vesta as the calibration source with a continuum model provided by T. Mueller (Mueller and Lagerros, ).

Interferograms taken on Vesta, Neptune and Uranus showed that, with the SCAL source off, the signal at the central peak does not saturate, or at most only a few samples are saturated, once the detector bias is correctly set. The fact that we do not require SCAL to null the telescope emission is a consequence of the lower total emission from the telescope and straylight compared to the expected values used in the initial design of the SPIRE instrument. Given this, and that using SCAL adds to the photon noise in the measurement, we have decided not to use the SCAL source in routine observations. An additional benefit of this mode of operation is that, with SCAL off it, and the rest of the instrument, are at a temperature between 4.5 and 5 K and the thermal emission from these components is limited to frequencies only detectable in the SLW band.

The standard calibration for SPIRE spectra is based on a point source. An alternative calibration can be derived using the telescope which is appropriate for a source that entirely fills the detector field of view. A correction based on measurements of the spectrometer beam shape is required to convert from one to the other (see Sect. ). It should also be noted that instrument dependent spectral features in the passband of the spectrometer will change as a function of source extent.

# Testing the Astronomical Models

To test the accuracy of the models used for both the spectrometer and photometer we have derived the calibration using an alternative method where we use the telescope itself plus knowledge of the instrument throughput (see Sect. ). As discussed in Sect. the telescope is represented by a blackbody with a well known relative dependence between emissivity and wavelength even if the absolute overall emissivity may have an uncertainty up to 30% (see Fischer et al. ). We can therefore generate an RSRF from the measured telescope spectrum on dark sky and use this to calibrate the spectra obtained on Vesta, Neptune and Uranus and compare to the models to check for self consistency. The results together with continuum models are shown in Fig. . The basic comparison is with the Neptune model of Moreno which shows an excellent agreement across all wavelengths giving confidence in both the relative and absolute level of the telescope emissivity. The Uranus model shown is also that of Moreno, but here we have increased the continuum by 10 Jy at all frequencies to give a better fit to the spectrum calibrated from the telescope. Reapplying a calibration file derived from this model of Uranus to Neptune shows that the Neptune measurement agrees with the model to within 5-10% across the full SPIRE band. The Vesta model as provided by Mueller was shown to be  15-20% lower than the measurement calibrated using Neptune and so we have increased the model spectrum that is used in the derivation of the spectrometer pipeline calibration by 1.2. We discuss the accuracy of the pipeline derived spectra further in Sect. .

# Beam and astrometric calibration

## Photometer

The measurement of the beam size of the photometer has been discussed in Sect. . The average beam derived for the three arrays has a near Gaussian core with a width of 18.1, 25.2 and 36.6 arcsec respectively all with an uncertainty of $\pm$`<!-- -->`{=html}5%. However, the full PSF outside of this core has much structure due to the diffraction from the secondary mirror and its support structure and it is recommended to use equivalent beam areas of 501, 944 and 1924 arcsec$^2$ respectively to calculate surface brightness when observing extended sources (SPIRE Observer's Manual ). These values are quoted as the mean for all detectors in an array, there is some minor variation across the arrays as demonstrated in Fig. which shows the centroid positions and the relative beam sizes from a point source scanned across the PSW array. The variation in beam size has been found to vary by no more than 10% for the majority of detectors across any of the three arrays. In the example shown here the maximum and minimum beam sizes of the nominally operational detectors are  17.3 and  19.5 arcsec. The apparently large beamsizes in two of the detectors are caused by their very slow time response and these detectors are excluded from the map making. Measurement of the beam width versus wavelength within the bands is not possible in flight but characterisation was performed during ground test campaigns (Ferlet et al. ).

The measured positions of the detectors in each array are folded into the map reconstruction algorithms together with the knowledge of the spacecraft pointing. Comparison between fields with sources seen at other wavelengths (mostly from *Spitzer and radio catalogues) shows that, with the present knowledge and algorithms, the typical astrometric accuracy of SPIRE maps processed with the current version of the *Herschel processing pipeline is around 4 arcsec. The cause of part of the inaccuracy has already been identified and subsequent versions of the pipeline have an astrometric accuracy of 2 arcsec or better. This should be compared to the telescope relative pointing error of 0.2 arcsec 1-sigma over 1 minute and the spatial relative pointing error of 1.5 arcec (Pilbratt et al. ).**

## Spectrometer

The shape and extent of the beam in the spectrometer has been subject to much study (Ferlet et al. ) both on the ground and in flight. The beam size versus frequency has been measured directly by taking medium resolution spectra on a point source of Neptune placed at different locations with respect to the central detectors by stepping the satellite position. The results are shown in Fig. . The highly structured nature of the variation in beam size with frequency is expected from the multi-moded feedhorns used for the spectrometer arrays. The beams are only Gaussian at the low frequency end of each band and methods for conversion between calibration valid for point and extended sources require further investigation.

# Wavelength calibration

## Photometer

The central wavelengths for the three photometer bands are derived from pre-flight measurements of the passbands of the detectors and the instrument optical filter chain (Spencer , SPIRE Observers Manual ). As described in Griffin et al. () these are integrated over the waveband assuming a flat spectrum to give band centres at 250, 352 and 504 $\mu$m. The edges of the bands as defined at 50% of the average in band transmission are 211-290 $\mu$m (PSW), 297-405 $\mu$m (PMW), and 409-611 $\mu$m (PLW). The filter profiles themselves are available as standard calibration products within the *Herschel  processing environment.*

## Spectrometer

The accuracy of the wavelength calibration was derived from line fits to $^{12}$CO lines in the spectra of 5 sources: Orion Bar, CRL618, NGC7023, S106 and DR21. The spacecraft velocity and the source velocity were subtracted from the fitted line positions to give the residual velocity offset. Using this method we found a calibration accuracy within 1/10th of a resolution element across the spectrometer band as defined by the FWHM of the Sinc profile - i.e. 1.207$\times 0.04$ cm$^{-1}$ - and that the average line centre offset was 27 km s$^{-1}$. Much of this offset can be explained by the detailed physics of the interferometer (Spencer et al. ) which is now understood and we expect that all instrumental offsets will be removed for all detectors across the field of view in future pipeline processing. The two spectrometer channels' transmission profiles fall rapidly at the edges and we define the trustworthy band edges for the central detectors of the two channels to be 14.9-33.0 cm$^{-1}$ and 32.0-51.5 cm$^{-1}$. These edges vary somewhat depending on the position of detectors across the field of view but nowhere do they change by more than 0.2 cm$^{-1}$.

# Discussion and conclusions

The most significant part of calibrating an astronomical instrument is the accurate conversion from measured signal to source flux density. In this paper we have given an overview of how this is currently carried out for the *Herschel-SPIRE instrument. We have demonstrated an inter-comparison between the models used for astronomical sources and the *Herschel telescope emission spectrum that gives a final agreement between the Neptune model used as the primary calibration standard and our measured spectrum of between 5-10%.**

The calibrator source model is only one component of the final uncertainty and we have tested the conversion parameters for the photometer Volts to Jy using a Monte Carlo approach to determine the variance in the calibration curve as a function of voltage. This method involved creating multiple calibration curves of $1/\Delta V_{PCAL}$ versus $V$ with addition random noise representative of the measurement errors and measuring the effect on the uncertainty in derived calibration coefficients. This translates into uncertainties that are $\leq$ 0.1% for all detectors. This is much smaller than the other two sources of error, which are the uncertainty in the model for Neptune (5-10%) and the uncertainty in peak flux density measurements for Neptune derived from the maps after processing. Taking these factors into account for data processed with the current processing software the overall photometric error can be assumed to be 15%.

For the spectrometer photometric calibration a comparison of the pipeline derived Neptune spectrum based on Vesta to the Neptune model shows a typical variation of 15-20% in the SSW band and 20-30% in the SLW band above 20 cm$^{-1}$ and up to 50% below 20 cm$^{-1}$. It should be noted that variations in the temperature of the instrument itself between the measurement of the dark spectrum and the source can cause extra flux to be observed in the SLW band. The net result is an additive term in the spectrum below 25 cm$^{-1}$ of the order of a few Jy which can leave significant positive or negative extraneous flux density levels in the source spectrum especially for faint objects. With careful interactive processing it is possible to remove this signature and obtain accurate source flux densities (Ivison et al. , Bocklee-Morvan et al. ). Further investigation on this aspect of the spectrometer calibration is ongoing and will be reported in future work.

The calibration of the data from the SPIRE instrument processed through the current data reduction pipeline software is already of a standard more than sufficient to allow significant scientific conclusions to be drawn from the observations. Future work will concentrate on improving the photometry using other astronomical sources and source models, including cross calibration with other facilities such as the *Planck satellite. We will also work on improvements to the data processing tools, the astrometry of the maps through better attitude reconstruction and the wavelength calibration of the spectrometer through more sophisticated instrument performance models.*

[^1]: Herschel is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with participation from NASA.


abstract: We briefly summarize the work and discussions that occurred during the Working Group 8 sessions of the EAAC 2017, dedicated to advanced and novel accelerators for high energy physics applications.
address: LPGP CNRS, Université Paris-Sud, Université Paris Saclay, Orsay, France; Max Planck Institute for Physics, München, Germany; CERN, Geneva, Switzerland; Lawrence Berkeley National Laboratory, Berkeley, California, USA; Tsinghua University, Beijing China
author: B. Cros; P. Muggli; C. B. Schroeder; C. Tang
title: Summary of Working Group 8: Advanced and Novel Accelerators for High Energy Physics

# Introduction

Working Group 8 (WG8), Advanced and Novel Accelerators for High Energy Physics, was created as a result of the Advanced and Novel Accelerators for High Energy Physics Roadmap (ANAR) Workshop 2017 that took place at CERN, April 25-28, 2017 . The workshop was organized at the initiative of the ICFA subpanel on Advanced and Novel Accelerators chaired by B. Cros. The workshop was a first step towards federating the fragmented advanced and novel accelerator (ANA) community, to define an international roadmap towards colliders based on advanced accelerator concepts, including intermediate milestones, and to discuss the needs for international coordination. A preliminary report summarizing the state-of-the-art of ANAs and the proposed future steps can be found at .

Working Group 8 was organized in two specific parallel sessions and one joint session with Working Group 1 - Electron beams from plasmas. The WG8 charge was to examine key challenges, discuss suitable concepts, and identify topics for future R&D or innovation, including electron and positron sources, damping rings, optics between stages, acceleration of positrons, luminosity goals, final focus, and overall efficiency.

# General Discussion

The work started with a general discussion led by A. Seryi who emphasized the need for our ANA community to take advantage of the vast knowledge of the conventional accelerators community. That community has built working accelerators, but has also developed Conceptual Design Reports (CDRs) and Technical Design Reports (TDRs) for future colliders that are ready to be built, such as the International Linear Collider (ILC) and the Compact Linear Collider (CLIC).

In the ILC case, the enterprise started with a group of motivated scientists who produced more and more refined designs that were periodically and critically reviewed by the community over many years. These designs were organized around systems (injector, damping rings, linac, beam dynamics, beam delivery system, etc.). Each system was analyzed, and for each system issues were identified and ranked. Required R&D programs relevant to collider criteria such as feasibility, reliability, production, and technical and cost optimization were defined and carried out. The work started with a number of options for each system, e.g., warm and superconducting technology for the linac. From the many designs that were generated (Tesla, JLC, NLC, etc.), one was eventually selected to become the ILC.

The ILC process toward a collider could serve as model for the design of an ANA-based advanced linear collider (ALC). Indeed, at this time the community includes four major ANA R&D efforts: Laser-driven plasma wakefield accelerators (LWFA), particle-driven plasma wakefield accelerators (PWFA), structure wakefield accelerators (SWFA), and dielectric laser accelerators (DLA), all with potentials and challenges, as summarized, for example, in the ANAR report .

This discussion helped define the structure of the next workshop organized by ICFA-ANA and the ALEGRO (Advanced LinEar collider study GROup) group , to be hosted by the University of Oxford and the John Adams Institute March 26-29, 2018. The upcoming ALEGRO workshop will have seven main working groups: Physics Case (PC), Collider machine design/definitions (CMD), Theory, Modeling, Simulations (TMS), LWFA, PWFA, SWFA, DLA, and a joint sub-WG on positron acceleration (PAC).

# Plasma-based ANAs

In this section, we report on the conclusions of the presentations discussing acceleration in plasma.

In plasma-based ANAs operating in the bubble regime, the transverse fields are of the same order as longitudinal fields. Whereas it can be seen as an advantage to channel bunches over long distances, it also means that the particles emit betatron radiation. The radiated power depends on the particle energy and amplitude of oscillation about the beam axis. Calculations performed by V. Shpakov show that, per unit length, the radiation loss is always smaller than the energy gain and betatron radiation loss is therefore not a limitation to energy gain. As expected, emission betatron radiation together with acceleration lead to some cooling of the bunch emittance. However, radiation loss causes increase of the energy spread, which decreases the overall quality of the accelerated bunch.

Under-dense plasma lenses have very strong focusing gradients compared to magnetic lenses. The focusing gradient, proportional to the local plasma electron density, can in principle be tailored along the lens. Increasing the focusing gradient *adiabatically* along the length allows for minimization of synchrotron radiation and overcoming of the Oide limit. Calculations, presented by F. Filippi, show that for an experiment at the INFN SPARC-LAB, an exponential increase of the plasma density by a factor of 400 over 6 cm and starting with a density of 5$\times$`<!-- -->`{=html}10$^{15}$ cm$^{-3}$ would be needed to focus a 110 MeV electron bunch from an initial 75 $\mu$m radius to $\sim
5\mu$m. Using a capillary discharge with a tapered entrance tube could be used for the density tapering. The size of the focused bunch could be inferred from the spectrum of the radiation emitted by the electrons in the constant density plasma inside the capillary.

Nikolay Andreev discussed external injection and acceleration of short electron bunches in plasma waves generated by intense femtosecond lasers in plasma channels. The dependence of the bunch length and energy spread was studied as a function of the phase of the plasma wave in the quasi-linear regime where the electrons are injected. For the parameters studied, when an electron bunch is injected with a size of the order of 20% the plasma wave length, at the maximum of the wave potential, the energy gain is or the order of 2 GeV over 20 cm of propagation with an rms energy spread smaller than $0.01$. This energy spread is achieved for charges of the order of a few pC; for larger charges, beam loading leads to increased energy spread. Mechanisms relevant to multi-stage accelerators, such as the preservation of emittance and spin polarization of electron beams between plasma stages were examined by D. Pugacheva. Multi-stage LWFA was shown to preserve beam polarization, as the theoretical prediction gives 0.02% for depolarization after acceleration up to 1 TeV, which is suitable for high-energy-physics applications.

Driving plasma wakefields with a positively charged particle bunch raises a number of challenges related to the nonlinearity of the transverse forces acting on the bunch when reaching the non-linear regime. Some of the challenges can be overcome by using a hollow plasma channel rather than a uniform plasma. High-energy proton bunches are long, which means that one has to rely on self-modulation to produce a train of bunches that can then drive large amplitude wakefields. However, the self-modulation process operates by de-focusing a large fraction of the protons. Y. Li proposed to deal with these two issues by driving wakefields with a pre-formed train of short proton bunches in a hollow plasma channel. Initial estimates also show that this scheme is relatively weakly sensitive to variations of the plasma channel.

# DLA

The dielectric laser accelerator (DLA) scheme has shown accelerating gradients exceeding GV/m, however over rather short accelerating distances and with small energy gains (sub-MeV). One of the advantages of the DLA is that optical elements can be used for all the functions necessary for an accelerator (power source, focusing, diagnostics, etc.), potentially allowing for an "on-chip" accelerator. As with other ANA concepts, one of the next steps for DLA is a global concept that includes an injector and a number of acceleration stages to produce high-energy, high-quality bunches. The DLA community, gathered in the ACHIP collaboration, is addressing challenges related to developing a high-energy accelerator on a chip.

J. England summarized the ACHIP program, including issues such as acceleration from rest to relativistic energies, micro-bunching at the laser wavelength scale, powering of multiple structures, compensation of globally de-focusing transverse forces during the weakly-relativistic part of the acceleration process, etc. The results of these studies are continuously integrated into updated versions of the "straw-man" concept collider based on DLA.

# Simulations

A key component for the progress towards an Advanced Linear Collider based on one of the ANAs is the availability of computationally efficient simulation tools. Simulations are particularly challenging for plasma-based ANAs.

Most simulation codes use explicit particle-in-cell (PIC) algorithms that are particularly consumming in terms of computer memory and time. Reduced models (quasi-static, hybrid, etc.) and numerical methods (boosted frame, azimuthal decomposition, etc.) can provide significant time and memory savings, without loss of physical fidelity. However, at present, the full modeling of a single multi-GeV accelerator stage in 3D remains extremely challenging. J.-L. Vay presented a very ambitious program, a seven-year, exascale modeling project, based on WarpX. The goal of the project aims at a full numerical description (3D) of 100 plasma-based accelerator stages, toward a 1 TeV collider design. The project plan includes modeling of a single plasma accelerator stage with static mesh refinement by 2017, modeling a single GeV-scale plasma accelerator stage by 2018, and modeling of multiple GeV-scale plasma accelerator stages by 2020. WarpX is an up-graded version of the WARP, coupled to the latest adaptive mesh refinement algorthms. WARP is an open-access code, and WarpX is scheduled to be released to the public by 2018. The plans also include visualization tools as well as interfaces for analysis with various programs. This is a very important and welcome initiative.

# Facilities

During the joint session with Working Group 1 (Electron beams from plasmas) of some of the present and planned ANA experimental facilities were described.

EuPRAXIA is a design study for a 5 GeV electron beam research facility for light source (FEL) and high-energy physics applications. The major goal for EuPRAXIA is to design and, in a second phase, build a compact European advanced accelerator with superior beam quality. P. A. Walker presented a designed and preliminary global parameters of a generic facility that uses ANA technology to drive an FEL. The facility would generate 5 GeV electron beams and include two users areas: one for the FEL and one for high-energy physics applications. The design accommodates options for the electron source (RF-gun with x-band linac, LWFA), and for the accelerator (LWFA, PWFA). Drive lasers and klystrons would be located in an upstairs mezzanine, whereas accelerators, undulators and experimental areas would be located on the ground floor.

At the moment there are five European facilities that could candidate to host the construction of a facility resulting from Eupraxia design study.

M. Ferrario presented plans for a new compact FEL facility at INFN-LNF that would be compatible with the EuPRAXIA design. The plans include all the options envisaged by EuPRAXIA, but also the option to have a 1 GeV x-band linac that could drive the FEL without plasma acceleration. The TDR is expected to be completed by the end of 2017. The facility would be an upgrade of the current SPARC-LAB that already operates with an s- and c-band accelerator to drive an FEL. In addition, the high-power laser FLAME, presently used for LWFA studies, would be upgraded with the aim of externally injecting and accelerating the bunch from the RF-gun. Plasma sources are already available for PWFA and plasma lens studies. Preliminary simulation results indicate that beams produced in the current SPARC-LAB, and thus its upgraded version, from plasma accelerators would lead to lasing of the FEL in the water window.

FLASHForward is the future-oriented wakefield accelerator research and development facility at FLASH, DESY. The main goal of the project is to explore PWFA. At FLASHForward, a PWFA witness bunch can be generated by a plasma cathode or from the FLASH linac. Single bunches can in principle be tailored for large transformer ratio and optimal beam loading experiments. Studies will also focus on the hosing instability, a effect that could disrupt beams in very long acceleration sections. In a second phase (starting 2020), undulators will be added to the beamline to drive an drive an FEL using an electron bunch produced in the PWFA. The work is supported by a strong simulation program, including start-to-end simulations. In his presentation J. Osterhoff described how many of the critical scientific challenges faced by the PWFA, and other ANAs, can be addressed at FLASHForward.

ARES is another facility at DESY studying LWFA and possibly DLA for FELs and other applications. The beam produced by the ARES linac would be injected into the LWFA. Towards this end, simulations including deleterious effects such as CSR, presented by J. Zhu, showed that after compression the ARES bunch could reach sub-fs length, which would make the bunch suitable for injection into an LWFA. These short bunches contained 0.8 pC, with 600 A current. Longer bunches with larger currents can also be produced and, perhaps, are more suitable to drive the FEL. Maintaining a low energy spread suitable for FEL lasing is a challenge, as in all LWFA designs. In the ARES case, simulations showed space charge effectscan limit the bunch current and the beam matching into the plasma.

The MAX IV is a 3 GeV light source storage ring located in Lund (Sweden). The 3 GeV beam from the linac is available for ANA experiments, when not used for filling of the ring. Simulation results presented by O. Lundh show that the 100 pC bunch can be compressed to reach more than 10 kA of peak current. However, in this case emittance growth becomes a serious issue. Less aggressive compression to 1.5 kA produces a bunch with low emittance and relative energy spread. Driving the gun photo-cathode with two laser pulses could lead, after two stages of compression to a drive/witness bunch train that would be very interesting for PWFA experiments. In particular, simulations showed that a loaded accelerating field could reach the 10 GV/m level. Plans for PWFA experiments are awaiting endorsement by the MAX IV Laboratory.

# Conclusions

Discussions during the working group re-enforced the need for active collaborations, community-gathering behind a global high-energy physics application of ANAs, continuous critical analysis of the various "straw-man" designs for an ALC, and need for leadership in this endeavor. It is clear that there is much to learn from the conventional accelerator community's experiences builing accelerators and developing TDRs for future linear colliders (e.g., CLIC and ILC). It is also clear that there is an urgent need to involve this community in the development of ANAs and also in the development of a global collider concept starting at the injector and ending at the detector and physics case. This model will be used for the up-coming ALEGRO workshop and for the generation of a document as an input for the next European Strategy on High-Energy Physics.

From the scientific point of view, presentations in WG8 showed that significant progress is being made towards high-energy accelerator concepts. The most attractive feature of ANAs remains the large accelerating gradients.

There are now roadmaps for ANA development that have been drawn in the US and in Europe . There is a European facility design study (EuPRAXIA) focusing on a plasma-based ANA to produce high-quality beams. A number of facilities are either operating, coming on-line, in the building phase, or in the planing phase. All of these facilities are geared toward light sources or high-energy physics applications. They will all contribute to addressing the major scientific challenges that have been identified by the community, most recently at the ANAR Workshop. ALEGRO will provide a forum for organizing the global ANA efforts towards producing high-quality and high energy beams, as well as deliver a document for the next European Strategy input and for an ALC.

# References


abstract: This letter presents a revised radiative transfer model for the infrared (IR) emission of active galactic nuclei (AGN). While current models assume that the IR is emitted from a dusty torus in the equatorial plane of the AGN, spatially resolved observations indicate that the majority of the IR emission from $\la$`<!-- -->`{=html}100 pc in many AGN originates from the polar region, contradicting classical torus models. The new model *CAT3D-WIND* builds upon the suggestion that the dusty gas around the AGN consists of an inflowing disk and an outflowing wind. Here, it is demonstrated that (1) such disk+wind models cover overall a similar parameter range of observed spectral features in the IR as classical clumpy torus models, e.g. the silicate feature strengths and mid-IR spectral slopes, (2) they reproduce the $3-5\,\micron$ bump observed in many type 1 AGN unlike torus models, and (3) they are able to explain polar emission features seen in IR interferometry, even for type 1 AGN at relatively low inclination, as demonstrated for NGC3783. These characteristics make it possible to reconcile radiative transfer models with observations and provide further evidence of a two-component parsec-scaled dusty medium around AGN: the disk gives rise to the $3-5\,\micron$ near-IR component, while the wind produces the mid-IR emission. The model SEDs will be made available for download.
author: Sebastian F. Hönig; Makoto Kishimoto
title: Dusty winds in active galactic nuclei: reconciling observations with models

# Introduction

It was held that the infrared (IR) emission from active galactic nuclei (AGN) originates from a "dusty torus" around the central black hole and accretion disk . This geometrically thick structure would explain the relative strength of the IR emission as well as the angle-dependent obscuration of the accretion disk and broad line region, leading to the well-established type 1/type 2 dichotomy of unobscured and obscured objects. Radiative transfer models based on such a toroidal geometry have proven successful in reproducing the IR spectral energy distribution (SED) of AGN . However, given the milliarcsecond sizes of the IR emission in most AGN, the true geometric distribution of the emission was never scrutinised up until recently.

Advances in IR interferometry made it possible to resolve the near- and mid-IR emission in about three dozen AGN . A growing number of nearby Seyfert galaxies have a sufficient interferometric measurements to determine the wavelength- and angular-dependent size of the emission. It has been shown that contrary to the classical torus picture, the bulk of the mid-IR emission originates from the polar region of the AGN rather than the equatorial plane . Detailed studies of the type 1 AGN NGC3783 and the type 2 sources NGC 1068 and the Circinus galaxy imply that the IR emission emerges from at least two distinct parsec-scale structures: (1) a compact, geometrically thin disk in the equatorial region of the AGN, and (2) an extended, elongated polar structure, which is cospatial with the outflow region of the AGN on larger scales . The disk seems to be contributing a relatively larger fraction to shorter wavelengths, while the polar region dominates the overall IR emission energetically, which requires a significant covering fraction. Individual aspects of this behaviour may be explained with radiative transfer effects within a torus; however, the combination of extension, direction and emitted radiation strongly disfavours current torus models, either smooth or clumpy.

One of the hypotheses put forward to explain this structure leans on the idea that matter accreting onto the AGN is subject to radiative and hydrodynamic pressure producing strong outflows . In the specific picture outlined in , the disk component represents the inflowing dusty gas. Near the sublimation radius, the radiation pressure from both the AGN and the dust itself causes material to lift up and form a hollow outflow cone. Both outflow and disk contribute to the obscuration. This paper presents the radiative transfer model *CAT3D-WIND* of a clumpy disk and outflow with the geometry and dust composition built upon the hypothesis outlined in . It will be shown that such a configuration reproduces key observational characteristics of AGN, including the overall IR SEDs, the type 1/type 2 dichotomy, the observed disk plus polar emission structure, and the $3-5\,\micron$ bump seen in many unobscured AGN.

# CAT3D-WIND: disk + wind model description

The *CAT3D-WIND* model is based on the parametrisation of a clumpy dust torus used in the underlying *CAT3D* model . It distributes dust clumps according to a radial power law $\propto r^a$, where $a$ is the power-law index and $r$ denotes the distance from the AGN in units of the sublimation radius $r_\mathrm{sub}$. The sublimation radius forms the inner boundary of the dust distribution where the dust becomes too hot to survive. The vertical distribution is characterised by a Gaussian distribution $\propto \exp(-z^2/(2(hr)^2))$ with dimensionless scale height $h$ and vertical distance from the mid-plane $z$ in units of $r_\mathrm{sub}$. The model space is limited from $r_\mathrm{sub}$ as the inner boundary to the outer radius $R_\mathrm{out}$. Please note that $R_\mathrm{out}$ should not be considered a free model parameter, but rather a minimum boundary that has to be chosen to not cause an artificial brightness cut-off in the wavelength range of interest . Obscuration is defined by the average number of dust clouds $N_0$ along an equatorial line-of-sight. The clouds are characterised by their optical depth $\tau_\mathrm{cl}$ and radius $R_\mathrm{cl}$, but the size of the clouds only controls the total number of clouds and not the typical model SED . Moreover, as long as the cloud is optically thick in the IR, which most clumpy models assume, the dominating emission from the hot side of the clouds is insensitive to the exact choice of $\tau_\mathrm{cl}$. In summary, these 6 parameters fully define typical clumpy torus models.

Classical torus models do not account for the observed predominance of polar emission in the mid-IR emission of AGN. Thus, a second component is added to the standard model in the form of a polar outflow. This structure is modelled as a hollow cone and characterised by three parameters: (1) the radial distribution of dust clouds in the wind $a_w$, (2) the half-opening angle of the wind $\theta_w$, and (3) its angular width $\sigma_\theta$. Finally, a wind-to-disk ratio $f_{wd}$ defines the ratio between number of clouds along the cone and $N_0$. Note that converting $f_{wd}$ to a mass ratio would require knowledge of $\tau_\mathrm{cl}$ and $R_\mathrm{cl}$. However, these are not constrained by the SEDs. For each set of model parameters, the emission as seen from inclinations of $0\degr-90\degr$ is simulated, in steps of $15\degr$. An illustration of the dust distribution and typical view of a disk+wind is shown in Fig. .

The new model makes use of a more physical description of dust sublimation near the AGN than the commonly assumed single sublimation radius. It accounts for the facts that for any given density, the silicate sublimation temperature is lower than the graphite sublimation temperature and that small grains are hotter for a fixed distance from the AGN. When the dust heats above 1200 K, silicates are removed from the dust composition, leaving only graphites that can heat up to 1900 K. In the intermediate stage from 1200 K leading up to 1900 K, the smallest graphite grains are removed, so that at the innermost radius where dust can survive only grains with a size between 0.075$\,\micron$ and 1$\,\micron$ are present. In this framework, it is possible to reconcile the observed small near-IR reverberation mapping sizes and near-IR interferometry sizes with dust sublimation physics .

One specific aspect of the disk+wind hypothesis put forward in is that the wind is launched near the sublimation zone of the dusty disk. Thus, the chemical composition of the dust in the wind is expected to be very similar to the one seen at about the sublimation radius. Accordingly, the model dust clouds in the polar region are devoid of silicates and small grains. This will also help near edge-on views of the disk+wind structure to display silicate absorption features at $\sim$`<!-- -->`{=html}10 $\micron$ from the disk despite the direct view to hotter dust in the outflow region. Wind-only models with standard silicate and graphite mixtures predominantly show silicate emission features . This is confirmed with the presented models when using normal ISM dust for the wind clouds. Such models show exclusively silicate emission features for the range of observed SED mid-IR slopes of both type 1 and type 2 AGN (see Sect. ).

| disk parameters |  |
|:--:|:---|
| $a$ | $\left[-0.5,-1.0,-1.5,-2.0,-2.5,-3.0\right]$ |
| $N_0$ | $\left[5,7.5,10\right]$ |
| $h$ | $\left[0.1,0.2,0.3,0.4,0.5\right]$ |
| wind parameters |  |
| $a_w$ | $\left[-0.5,-1.0,-1.5,-2.0,-2.5\right]$ |
| $\theta_w$ | $\left[45\degr,30\degr\right]$ |
| $\sigma_\theta$ | $\left[7.5,10,15\right]$ |
| global parameters |  |
| $f_{wd}$ | $\left[0.15,0.3,0.45,0.6,0.75\right]$ |
|  | additionally $\left[1,1.25,1.5,1.75\right]$ for $\theta_w=30\degr$ |
| $R_\mathrm{out}$ $^{(1,2)}$ | $500$ (large graphite) / $151$ (ISM) |
| $R_\mathrm{cl}$ $^{(1)}$ | $0.035 \cdot r$ |
| $\tau_V$ | $50$ |
| line of sight |  |
| inclination | $\left[0\degr,15\degr,30\degr,45\degr,60\degr,75\degr,90\degr\right]$ |

Range of model parameters used to create the parameter grid.

*— Notes:* For parameter definition see main text. $^{(1)}$ in units of $r_\mathrm{sub}$. $^{(2)}$ The radius where large graphite grains sublimate at 1900 K is a factor of 3.3 smaller than the classical definition of the sublimation radius. $R_\mathrm{out}$ is set to match the classical *CAT3D* value that does not result in artificial cut-offs in the mid-IR .

# Results and discussion

The parameter space of the model has been explored in the range of parameters as listed in Table . A total of 132,300 model SEDs have been simulated. It should be noted that the range of $f_{wd}$ is larger for $\theta_w = 30\degr$ than for $45\degr$. In general, the covered surface for a cone with half-opening angle closer to the polar axis is smaller than for a cone with a larger half-opening angle. To counter this effect, the effectively filled area in the cone for the models with smaller half-opening angle has been increased.

## SED shapes and silicate features

To obtain a swift overview of the typical SED shapes of all the models, each model SED has been broken down into three parameters: (1) the spectral index $\alpha_\mathrm{MIR}$ of the mid-IR emission, expressed by the flux ratio between $14\,\micron$ and $8\,\micron$, (2) the spectral index $\alpha_\mathrm{IR}$ from the near-IR to the mid-IR, expressed by the flux ratio between $6\,\micron$ and $3\,\micron$, and (3) the silicate feature optical depth $\tau_\mathrm{Si}$. Here, the spectral index is defined as $$\alpha = -\frac{\log\left(F_\nu(\lambda_j)/F_\nu(\lambda_i)\right)}{\log\left(\lambda_j/\lambda_i\right)}$$ where $F_\nu(\lambda)$ refers to the flux at wavelength $\lambda$ and $\lambda_j > \lambda_i$. Fig.  shows an overview of the range of SED shapes obtained by the *CAT3D-WIND* model based on the parameter grid in Table . In addition, a corresponding grid for the classical clumpy torus based on *CAT3D* is shown for comparison by setting $f_{wd}=0$. Greater spectral indices relate to bluer SEDs. For $\tau_\mathrm{Si}$, blue colours indicate silicate emission features and red colours silicate absorption features. The range of $\alpha_\mathrm{MIR}$ has been studied in the past for a variety of AGN classes . In general, Seyfert galaxies and quasars fall in the range of $-0.5 > \alpha_\mathrm{MIR} > -3$, with type 1 AGN showing slightly bluer SEDs than type 2 AGN. This range of observed parameters is well covered by both model grids, illustrating the general viability of the disk+wind models as well as the classical torus models.

One interesting aspect of the covered model space is the relation of $\alpha_\mathrm{MIR}$ and $\alpha_\mathrm{IR}$. Fig.  shows a black dashed line that marks equal spectral index in the near- and mid-IR. Generally, classic clumpy torus models cover the parameter space above this line, where the near-IR rises relatively steeply towards the mid-IR . Then, the SEDs flatten in the mid-IR and turn over longward of $\sim20-30\,\micron$. Such models struggle to reproduce the $3-5\,\micron$ bump that is seen in several type 1 AGN, including quasars . Those SEDs are characterised by a first, local emission peak in the near-IR regime of the SED before they become steeper again towards the global $20-30\,\micron$ peak. SEDs with a $3-5\,\micron$ bump are found below the dashed line where the near-to-mid-IR spectral index $\alpha_\mathrm{IR}$ is bluer than the mid-IR spectral index $\alpha_\mathrm{IR}$. Indeed, the disk+wind models of *CAT3D-WIND* naturally produce such SEDs for a significant number of parameter combinations. Most notably, the $3-5\,\micron$ bump emerge from models where the disk has a relatively steep dust cloud distribution of $-2 > a > -3$ while the distribution of clouds in the wind is much shallower with $-0.5 > a > -1.5$. In these cases, the near-IR emission is dominated by the hot dust in the disk and the wind primarily contributes at longer wavelengths, which is the scenario implied by interferometry observations of NGC3783 .

The silicate emission features in most of the *CAT3D-WIND* model SEDs are relatively shallow while the absorption features can be deeper. Qualitatively, this is consistent with observational constraints for both space-based and ground-based AGN samples . The distribution of observed feature strengths peak in the region of about $0.1 \la \tau_\mathrm{Si} \la -0.4$ for quasars, $0.3 \la \tau_\mathrm{Si} \la -0.3$ for Seyfert 1 galaxies, and $1 \la \tau_\mathrm{Si} \la -0.1$ for Seyfert 2 AGN. The model space of the presented parameter grid covers these values.

In summary, the model SEDs of the disk+wind model cover the observed SED shapes and silicate feature depths quite well, very similar to classical clumpy torus models. In addition and beyond the capabilities of the torus models, the disk+wind model is able to reproduce the observed $3-5\,\micron$ bump in type 1 AGN from a compact disk, while preserving the mid-IR bump produced by the wind. This makes it superior to classical clumpy torus models on the criterium of reproducing SEDs alone.

## Polar elongations in type 1 AGN: NGC3783 as a case study

Type 1 AGN NGC3783 displays a pronounced $3-5\,\micron$ bump (see Fig. ). Decomposing the SED of NGC3783 similarly to the *CAT3D-WIND* model SEDs results in $\alpha_\mathrm{MIR} = -2.05\pm0.20$, $\alpha_\mathrm{IR} = -1.0\pm0.2$, and $\tau_\mathrm{Si} = -0.07 \pm 0.06$. This puts NGC3783 in a range that is only covered by the disk+wind models, not the classical clumpy torus models (see grey contours in Fig.  left & right). In Fig.  we plot the range of model SEDs that fall within the uncertainty interval of these parameters as a grey band. The majority of these models are characterised by a compact, relatively flat disk ($a \sim -2.0\ldots-3$; $h \la 0.3$) and an extended wind ($a_w \ga -1$). Inclinations range from $0\degr \le i \le 60\degr$ indicating a moderately inclined AGN, as can be expected for a type 1. The SED-only models are inconclusive with regards to the half-opening angle of the outflow.

The major motivation for the disk+wind model stems from IR interferometry of nearby AGN. Therefore, viability of the new models rests on the ability to not only reproduce the SEDs but also the interferometry measurements. In general, it can be expected that a type 2 AGN seen edge on will indeed show elongation towards the polar region in these models. On the other hand, radiative transfer effects might also produce some polar emission in obscured AGN in classical torus models . The more challenging objects are moderately-inclined type 1 AGN where classical models would produce an elongated, disk-like structure with orientation along the disk plane.

NGC3783 can be considered a prototypical inclined type 1 AGN with dominating polar IR emission. It has been extensively covered with near- and mid-IR interferometry. report simple modelling of 41 VLTI MIDI observations in the mid-IR and 6 VLTI AMBER measurements in the near-IR . The IR SED shows the characteristic $3-5\,\micron$ bump with an estimated hot-dust covering factor of about 30%. The near-to-mid-IR photometric and interferometric data is best reproduced by two IR emitting components, characterised by two perpendicular elongated structures, with the polar-oriented dominating in the mid-IR. The mid-IR emission is very elongated with a major-to-minor axis ratio of $\sim$`<!-- -->`{=html}3:1, which is challenging to achieve in classical torus models for an unobscured type 1 AGN at moderate inclinations.

The best SED model fits (see grey contour in Fig. , right) were selected to simulate multi-wavelength images, which is computationally expensive. Interferometric visibilities were obtained for the model images at $12\,\micron$ and compared to the corresponding interferometric visibilities of NGC3783 . It is immediately found that the strong elongation seen in NGC3783 disfavour low inclinations at $0\degr-15\degr$, as expected. The strongest elongations are seen either at higher inclination, where the cone is seen from the side, or when the line-of-sight into the cone is close to its edge. Stronger polar elongations are found for $\theta_w = 30\degr$ than for $\theta_w = 45\degr$.

Fig.  shows model visibilities for several position angles for a representative model with $a=-3$, $a_w = -0.5$, $N_0=10$, $\theta_w = 38\degr$, $\sigma_\theta = 10\degr$ and $f_{wd}=1.2$ seen under an inclination of $i=30\degr$. This orientation and cone shape implied by this model are in between the suggestions of and , who infer $i=60\degr/\theta_w = 30\degr$ and $i=15\degr/\theta_w=45\degr$, respectively. While argue that kinematics disfavour the high-inclination solution, the strong elongation seen in IR interferometry might indicate a higher inclination. On the other hand, it is shown here that a configuration where the line-of-sight is close to the cone edge produces strong elongations as well.

## Anisotropy of the mid-IR emission

An optically thick dusty torus causes significant anisotropy of the emerging IR emission depending on inclination: More IR radiation escapes towards low inclinations than towards highly inclined line-of-sights. However, when taking the intrinsic X-ray emission as an isotropic tracer for the AGN luminosity, it has been shown that the mid-IR emission is surprisingly isotropic , unlike expected from smooth torus models. While this small scatter in the mid-IR-to-X-ray luminosity correlation can be partly explained in the framework of clumpy torus models, the range of torus parameters leading to low anisotropy is very restricted, which would imply a physical mechanism that restricts those parameters .

In the disk+wind picture, the hollow cone is visible to the observer at any inclination, naturally leading to low anisotropy of the mid-IR emission. Quantitatively, the peak of the luminosity-normalised flux distributions at $12\,\micron$ shifts by about 0.4 dex from $0\degr$ to $90\degr$. When considering typical inclinations of $30\degr$ for type 1 AGN and $75\degr$ for type 2s, the anisotropy is $0.24^{+0.11}_{-0.16}$ dex. This is consistent with estimates from the relation between X-ray and mid-IR luminosities in AGN, finding $<0.3$ dex and samples isotropically selected in the radio at $\sim0.22$ dex .

# Summary and Conclusions

In this letter, a new radiative transfer model *CAT3D-WIND* is presented to reconcile spatially resolved observations of AGN in the IR with model predictions. The model consists of a geometrically thin disk of optically-thick dust clumps and a hollow cone of putatively outflowing dust clouds. It is phenomenologically motivated by recent results in IR interferometry that hint towards such a two-component structure. In this picture, dust clouds are accreted in the plane of the disk. As dust partially sublimates near the sublimation radius (leading to a dearth of silicates), some dust clouds are lifted up by radiation pressure and flow out into the polar region. Based on a parametrisation of this picture, the following has been found:

-   The new disk+wind models cover a similar parameter space in terms of mid-IR spectral slope and strength of the silicate feature as classical clumpy torus models.

-   It is found that the new models are able to explain the $3-5\,\micron$ SED bump seen in many type 1 AGN. This local emission peak appears from the disk when it has a more compact distribution of dust clouds than the wind.

-   Many of the disk+wind models do indeed show polar elongation for inclinations typically associated with both type 1 and type 2 AGN. Indeed, a model has been found to simultaneously reproduce the IR SED and IR interferometry of the type 1 AGN NGC3783.

-   Exposure of the cone to the observer at all inclinations naturally results in a low degree of anisotropy of the mid-IR emission with respect to the inclination or viewing angle. This is consistent with the observed low scatter of the mid-IR/X-ray luminosity correlation.

The presented models are radiative transfer models of an empirical hypothesis. However, the physical mechanisms to drive the wind are not clear yet. It is possible that a combination of radiation pressure on the dust and neutral gas, both from the AGN radiation and from the IR-emitting medium itself, will play a role. Furthermore, hydrodynamic pressure close to the sublimation radius may help in launching such a wind. Work on radiative-hydrodynamical simulations will be needed to find a physical foundation for the observed polar structure . The SEDs presented in this paper will be made available for download at <http://cat3d.sungrazer.org>.

We thank the anonymous referee for helpful comments that improved the manuscript. SFH acknowledges support for this work from the European Research Council Horizon 2020 grant DUST-IN-THE-WIND (677117). Part of the numerical computations were done on the Sciama High Performance Compute (HPC) cluster which is supported by the ICG, SEPNet and the University of Portsmouth. Special thanks go to D. Asmus, P. Gandhi, C. Ricci, M. Stalevski, and K. Tristram for helpful discussions.