The late stellar evolutionary phases of low and intermediate-mass stars are strongly constrained by their mass-loss rates. The wind surrounding cool evolved stars frequently shows non-spherical features, thought to be due to an unseen companion orbiting the donor star. We study the morphology of the circumbinary envelope, in particular around oxygen-rich asymptotic giant branch (AGB) stars. We run a grid of 70 3D hydrodynamics simulations of a progressively accelerating wind propagating in the Roche potential formed by a mass-loosing evolved star in orbit with a main sequence companion. We resolve the flow structure both in the immediate vicinity of the secondary, where bow shocks, outflows and wind-captured disks form, and up to 40 orbital separations, where spiral arms, arcs and equatorial density enhancements develop. When the companion is deeply engulfed in the wind, the lower terminal wind speeds and more progressive wind acceleration around oxygen-rich AGB stars make them more prone than carbon-rich AGB stars to display more disturbed outflows, a disk-like structure around the companion and a wind concentrated in the orbital plane. In these configurations, a large fraction of the wind is captured by the companion which leads to a significant shrinking of the orbit over the mass-loss timescale, if the donor star is at least a few times more massive than its companion. Provided the companion has a mass of at least a tenth of the mass of the donor star, it can compress the wind in the orbital plane up to large distances. Our grid of models covers a wide scope of configurations function of the dust chemical content, the terminal wind speed relative to the orbital speed, the extension of the dust condensation region around the cool evolved star and the mass ratio. It provides a frame of reference to interpret high-resolution maps of the outflows surrounding cool evolved stars.

Radiative transfer is a key component in almost all astrophysical and cosmological simulations. We present Magritte: a modern open-source software library for 3D radiative transfer. It uses a deterministic ray-tracer and formal solver, i.e. it computes the radiation field by tracing rays through the model and solving the radiative transfer equation in its second-order form along a fixed set of rays originating from each point. Magritte can handle structured and unstructured input meshes, as well as smoothed-particle hydrodynamics (SPH) particle data. In this first paper, we describe the numerical implementation, semi-analytic tests and cross-code benchmarks for the non-LTE line radiative transfer module of Magritte. This module uses the radiative transfer solver to self-consistently determine the populations of the quantised energy levels of atoms and molecules using an accelerated Lambda iteration (ALI) scheme. We compare Magritte with the established radiative transfer solvers Ratran (1D) and Lime (3D) on the van Zadelhoff benchmark and present a first application to a simple Keplerian disc model. Comparing with Lime, we conclude that Magritte produces more accurate and more precise results, especially at high optical depth, and that it is faster.

Red supergiant stars are surrounded by a gaseous and dusty circumstellar environment created by their mass loss which spreads heavy elements into the interstellar medium. The structure and the dynamics of this envelope are crucial to understand the processes driving the red supergiant mass loss and the shaping of the pre-supernova ejecta. We have observed the emission from the CO $J = 2-1$ line from the red supergiant star $\mu$~Cep with the NOEMA interferometer. In the line the synthesized beam was $0.92 \times 0.72$~arcsec ($590 \times 462$~au at 641~pc). The continuum map shows only the unresolved contribution of the free-free emission of the star chromosphere. The continuum-subtracted channel maps reveal a very inhomogeneous and clumpy circumstellar environment. In particular, we detected a bright CO clump, as bright as the central source in the line, at 1.80~arcsec south-west from the star, in the blue channel maps. After a deprojection of the radial velocity assuming two different constant wind velocities, the observations were modelled using the 3D radiative transfer code \textsc{lime} to derive the characteristics of the different structures. We determine that the gaseous clumps observed around $\mu$~Cep are responsible for a mass loss rate of $(4.9 \pm 1.0) \times 10^{-7}~{\rm M}_\odot\,{\rm yr}^{-1}$, in addition to a spatially unresolved wind component with an estimated mass-loss rate of $2.0 \times 10^{-6}~{\rm M}_\odot\,{\rm yr}^{-1}$. Therefore, the clumps have a significant role in $\mu$~Cep's mass loss ($\ge 25 \%$). We cannot exclude that the unresolved central outflow may be made of smaller unresolved clumps.

In 1981, the idea of a superwind that ends the life of cool giant stars was proposed. Extreme OH/IR-stars develop superwinds with the highest mass-loss rates known so far, up to a few 10^(-4) Msun/yr, informing our understanding of the maximum mass-loss rate achieved during the Asymptotic Giant Branch (AGB) phase. A condundrum arises whereby the observationally determined duration of the superwind phase is too short for these stars to become white dwarfs. Here, we report on the detection of spiral structures around two cornerstone extreme OH/IR-stars, OH26.5+0.6 and OH30.1-0.7, identifying them as wide binary systems. Hydrodynamical simulations show that the companion's gravitational attraction creates an equatorial density enhancement mimicking a short extreme superwind phase, thereby solving the decades-old conundrum. This discovery restricts the maximum mass-loss rate of AGB stars around the single-scattering radiation-pressure limit of a few 10^(-5) Msun/yr. This brings about crucial implications for nucleosynthetic yields, planet survival, and the wind-driving mechanism.

High-resolution interferometric observations of the circumstellar environments of AGB stars show a variety of morphologies. Guided by the unusual carbon monoxide line profile of the AGB star EP Aquarii, we have observed its circumstellar environment with ALMA band 6 in cycle 4. We describe the morphological complexity of the CO, SiO and SO2 molecular emission. The CO emission exhibits the characteristics of a bi-conical wind with a bright nearly face-on spiral feature around the systemic velocity. This is the first convincing detection of a spiral morphology in an O-rich wind. Based on the offsets of the centers of the two bi-conical wind hemispheres, we deduce the position angle of the inclination axis to be ~150 deg measured counterclockwise from north. Based on the velocity width of the spiral signature we estimate the inclination angle of the system to be between 4 deg and 18 deg. The central emission zone exhibits a morphology that resembles simulations modelling the spiral-inducing wind-Roche-lobe-overflow mechanism. Though the spiral may be a (companion-induced) density enhancement in the stellar outflow, the extremely narrow width of the spiral signature in velocity space suggests that it may be a hydrodynamical perturbation in a face-on differentially rotating disk. The SiO emission does not show the spiral, but exhibits a local emission void approximately 0.5'' west of the continuum brightness peak. We hypothesise that this may be a local environment caused by the presence of a (stellar) companion with a mass of at most 0.1 Msol, based on its non-detection in the continuum...

A spectral scan of the circumstellar environment of the asymptotic giant branch (AGB) star R~Doradus was taken with ALMA in cycle 2 at frequencies between 335 and 362 GHz and with a spatial resolution of $\sim$150 milliarcseconds. Many molecular lines show a spatial offset between the blue and red shifted emission in the innermost regions of the wind. The position-velocity diagrams of this feature, in combination with previous SPHERE data and theoretical work point towards the presence of a compact differentially rotating disk, orientated nearly edge-on. We model the $^{\rm 28}$SiO ($v=1,~J=8\to7$) emission with a disk model. We estimate the disk mass and angular momentum to be $3 \times 10^{-6}$ Solar masses and $5 \times 10^{40}\ {\rm m^2 kg/s}$. The latter presents an `angular momentum problem' that may be solved by assuming that the disk is the result of wind-companion interactions with a companion of at least 2.5 earth masses, located at 6 AU, the tentatively determined location of the disk's inner rim. An isolated clump of emission is detected to the south-east with a velocity that is high compared to the previously determined terminal velocity of the wind. Its position and mean velocity suggest that it may be associated with a companion planet, located at the disk's inner rim.

Equatorial density enhancements (EDEs) are a very common astronomical phenomenon. Studies of the circumstellar environments (CSE) of young stellar objects and of evolved stars have shown that these objects often possess these features. These are believed to originate from different mechanisms, ranging from binary interactions to the gravitational collapse of interstellar material. Quantifying the effect of the presence of this type of EDE on the observables is essential for a correct interpretation of high-resolution data. We seek to investigate the manifestation in the observables of a circumstellar EDE, to assess which properties can be constrained, and to provide an intuitive bedrock on which to compare and interpret upcoming high-resolution data (e.g. \emph{ALMA} data) using 3D models. We develop a simplified analytical parametrised description of a 3D EDE, with a possible substructure such as warps, gaps, and spiral instabilities. In addition, different velocity fields (Keplerian, radial, super-Keplerian, sub-Keplerian and rigid rotation) are considered. The effect of a bipolar outflow is also investigated. The geometrical models are fed into the 3D radiative transfer code {\tt LIME}, that produces 3D intensity maps throughout velocity space. We investigate the spectral signature of the $J$=3$-$2 up to $J$=7$-$6 rotational transitions of CO in the models, as well as the spatial aspect of this emission by means of channel maps, wide-slit position-velocity (PV) diagrams, stereograms, and spectral lines. Additionally, we discuss methods of constraining the geometry of the EDE, the inclination, the mass-contrast between the EDE and the bipolar outflow, and the global velocity field. Finally, we simulated \emph{ALMA} observations...

The circumstellar environment of the AGB star L$_{\rm 2}$ Puppis was observed with ALMA in cycle 3, with a resolution of $15 \times 18 \rm\ mas$. The molecular emission shows a differentially rotating disk, inclined to a nearly edge-on position. In the first paper in this series (paper I) the molecular emission was analysed to accurately deduce the motion of the gas in the equatorial regions of the disk. In this work we model the optically thick $^{\rm 12}$CO $J=$3$-$2 and the optically thin $^{\rm 13}$CO $J=$3$-$2 rotational transition to constrain the physical conditions in the disk. To realise this effort we make use of the 3D NLTE radiative transfer code {\tt LIME}. The temperature structure and velocity structure show a high degree of complexity, both radially and vertically. The radial H$_{\rm 2}$ density profile in the disk plane is characterised by a power law with a slope of $-3.1$. We find a $^{\rm 12}$CO over $^{\rm 13}$CO abundance ratio of 10 inside the disk. Finally, estimations of the angular momentum in the disk surpass the expected available angular momentum of the star, strongly supporting the indirect detection of a compact binary companion reported in paper I. We estimate the mass of the companion to be around 1 Jupiter mass.

(abridged) In the inner regions of AGB outflows, several molecules have been detected with abundances much higher than those predicted from thermodynamic equilibrium (TE) chemical models. The presence of the majority of these species can be explained by shock-induced non-TE chemical models, where shocks caused by the pulsating star take the chemistry out of TE in the inner region. Moreover, a non-uniform density structure has been detected in several AGB outflows. A detailed parameter study on the quantitative effects of a non-homogeneous outflow has so far not been performed. We implement a porosity formalism for treating the increased leakage of light associated with radiation transport through a clumpy, porous medium. The effects from the altered UV radiation field penetration on the chemistry, accounting also for the increased reaction rates of two-body processes in the overdense clumps, are examined. We present a parameter study of the effect of clumping and porosity on the chemistry throughout the outflow. Both the higher density within the clumps and the increased UV radiation field penetration have an important impact on the chemistry, as they both alter the chemical pathways. The increased amount of UV radiation in the inner region leads to photodissociation of parent species, releasing the otherwise deficient elements. We find an increased abundance in the inner region of all species not expected to be present assuming TE chemistry, such as HCN in O-rich outflows, H$_2$O in C-rich outflows, and NH$_3$ in both. Outflows whose clumps have a large overdensity and that are very porous to the interstellar UV radiation field yield abundances comparable to those observed in O- and C-rich outflows for most of the unexpected species investigated. The inner wind abundances of H$_2$O in C-rich outflows and of NH$_3$ in O- and C-rich outflows are however underpredicted.

A spectral line and imaging survey of the low mass-loss rate AGB star R Dor (Mdot ~ 1e-7 Msun/yr) and the high mass-loss rate AGB star IK Tau (Mdot ~5e-6 Msun/yr) was made with ALMA between 335 and 362 GHz at a spatial resolution of ~150 mas, corresponding to the locus of the main dust formation region of both targets. Some 200 spectral features from 15 molecules (and their isotopologues) were observed, including rotational lines in both the ground and vibrationally excited states. Detected species include the gaseous precursors of dust grains such as SiO, AlO, AlOH, TiO, and TiO2. We present a spectral atlas for both stars and the parameters of all detected spectral features. A clear dichotomy for the sulphur chemistry is seen: while CS, SiS, SO, and SO2 are abundantly present in IK Tau, only SO and SO2 are detected in R Dor. Also other species such as NaCl, NS, AlO, and AlOH display a completely different behaviour. From some selected species, the minor isotopologues can be used to assess the isotopic ratios. The channel maps of many species prove that both large and small-scale inhomogeneities persist in the inner wind of both stars in the form of blobs, arcs, and/or a disk. The high sensitivity of ALMA allows us to spot the impact of these correlated density structures in the spectral line profiles. The spectral lines often display a half width at zero intensity much larger than expected from the terminal velocity, v_inf, previously derived for both objects (36 km/s versus v_inf ~17.7 km/s for IK Tau and 23 km/s versus v_inf ~5.5 km/s for R Dor). Both a more complex 3D morphology and a more forceful wind acceleration of the (underlying) isotropic wind can explain this trend. The formation of fractal grains in the region beyond ~400 mas can potentially account for the latter scenario. From the continuum map, we deduce a dust mass of ~3.7e-7 Msun for IK Tau and ~2e-8 Msun for R Dor.

We observed Betelgeuse using ALMA's extended configuration in band 7 (f~340 GHz, {\lambda}~0.88 mm), resulting in a very high angular resolution of 18 mas. Using a solid body rotation model of the 28SiO(v=2,J=8-7) line emission, we show that the supergiant is rotating with a projected equatorial velocity of v_eq sin i = 5.47 +/- 0.25 km/s at the equivalent continuum angular radius R_star = 29.50 +/- 0.14 mas. This corresponds to an angular rotation velocity of {\omega} sin i = (5.6 +/- 1.3) x 10^(-9) rad/s. The position angle of its north pole is PA = 48.0 +/- 3.5{\deg}. The rotation period of Betelgeuse is estimated to P/sin i = 36 +/- 8 years. The combination of our velocity measurement with previous observations in the ultraviolet shows that the chromosphere is co-rotating with the star up to a radius of ~10 au (45 mas or 1.5x the ALMA continuum radius). The coincidence of the position angle of the polar axis of Betelgeuse with that of the major ALMA continuum hot spot, a molecular plume, and a partial dust shell (from previous observations) suggests that focused mass loss is currently taking place in the polar region of the star. We propose that this hot spot corresponds to the location of a particularly strong "rogue" convection cell, which emits a focused molecular plume that subsequently condenses into dust at a few stellar radii. Rogue convection cells therefore appear to be an important factor shaping the anisotropic mass loss of red supergiants.

We study the sub-Keplerian rotation and dust content of the circumstellar material around the asymptotic giant branch (AGB) star L$_2$ Puppis. We find that the thermal pressure gradient alone cannot explain the observed rotation profile. We find that there is a family of possible dust populations for which radiation pressure can drive the observed sub-Keplerian rotation. This set of solutions is further constrained by the spectral energy distribution (SED) of the system, and we find that a dust-to-gas mass ratio of $\sim10^{-3}$ and a maximum grain size that decreases radially outwards can satisfy both the rotation curve and SED. These dust populations are dynamically tightly coupled to the gas azimuthally. However grains larger than $\sim0.5\mu$m are driven outward radially by radiation pressure at velocities $\sim5$kms$^{-1}$, which implies a dust replenishment rate of $\sim3\times10^{-9}$M$_\odot$yr$^{-1}$. This replenishment rate is consistent with observational estimates to within uncertainties. Coupling between the radial motion of the dust and gas is weak and hence the gas does not share in this rapid outward motion. Overall we conclude that radiation pressure is a capable and necessary mechanism to explain the observed rotation profile of L$_2$ Pup, and offers other additional constraints on the dust properties.

(abridged) Our aim is to determine the radial abundance profile of SiO and HCN throughout the stellar outflow of R Dor, an oxygen-rich AGB star with a low mass-loss rate. We have analysed molecular transitions of CO, SiO, and HCN measured with the APEX telescope and all three instruments on the Herschel Space Observatory, together with literature data. Photometric data and the infrared spectrum measured by ISO-SWS were used to constrain the dust component of the outflow. Using both continuum and line radiative transfer methods, a physical envelope model of both gas and dust was established. We have performed an analysis of the SiO and HCN molecular transitions in order to calculate their abundances. We have obtained an envelope model that describes the dust and the gas in the outflow, and determined the abundance of SiO and HCN throughout the region of the outflow probed by our molecular data. For SiO, we find that the initial abundance lies between $5.5 \times 10^{-5}$ and $6.0 \times 10^{-5}$ w.r.t. H$_2$. The abundance profile is constant up to $60\ \pm 10\ R_*$, after which it declines following a Gaussian profile with an $e$-folding radius of $3.5 \pm 0.5 \times 10^{13}$ cm. For HCN, we find an initial abundance of $5.0 \times 10^{-7}$ w.r.t. H$_2$. The Gaussian profile that describes the decline starts at the stellar surface and has an $e$-folding radius $r_e$ of $1.85 \pm 0.05 \times 10^{15}$ cm. We cannot to unambiguously identify the mechanism by which SiO is destroyed at $60\ \pm 10\ R_*$. The initial abundances found are larger than previously determined (except for one previous study on SiO), which might be due to the inclusion of higher-$J$ transitions. The difference in abundance for SiO and HCN compared to high mass-loss rate Mira star IK Tau might be due to different pulsation characteristics of the central star and/or a difference in dust condensation physics.

(abbreviated) We aim to constrain the dust formation histories in the winds of oxygen-rich AGB stars. We have obtained ALMA observations with a spatial resolution of 120x150 mas tracing the dust formation region of a low mass-loss rate and a high mass-loss rate AGB star, R Dor and IK Tau. Emission line profiles of AlO, AlOH and AlCl are detected and are used to derive a lower limit of atomic aluminium incorporated in molecules. We show that the gas-phase aluminium chemistry is completely different in both stars, with a remarkable difference in the AlO and AlOH abundance stratification. The amount of aluminium locked up in these 3 molecules is small, <=1.1e-7, for both stars, i.e. only <=2% of the total aluminium budget. This leaves ample of room for aluminium to be incorporated in grains. A fundamental result is that AlO and AlOH, being the direct precursors of alumina grains, are detected well beyond the onset of the dust condensation proving that the aluminium oxide condensation cycle is not fully efficient. The ALMA observations allow us to quantitatively assess the current generation of theoretical dynamical-chemical models for AGB winds. We discuss how the current proposed scenario of aluminium dust condensation for low mass-loss rate AGB stars at a distance of ~1.5 Rstar, in particular for the stars R Dor and W Hya, poses a challenge if one wishes to explain both the dust spectral features in the spectral energy distribution (SED), in interferometric data, and in polarized light signal. In particular, the estimated grain temperature of Al2O3 is too high for the grains to retain their amorphous structure. We propose that large gas-phase (Al2O3)n-clusters (n>34) can be the potential agents of the broad 11 micron feature in the SED and in the interferometric data and we explain how these large clusters can be formed.

The aim of this paper is to investigate the $^{17}$O/$^{18}$O ratio for a sample of AGB stars, containing M-, S- and C-type stars. These ratios are evaluated in relation to fundamental stellar evolution parameters: the stellar initial mass and pulsation period. Circumstellar $^{13}$C$^{16}$O, $^{12}$C$^{17}$O and $^{12}$C$^{18}$O line observations were obtained for a sample of nine stars with various single-dish long-wavelength facilities. Line intensity ratios are shown to relate directly to the surface $^{17}$O/$^{18}$O abundance ratio. Stellar evolution models predict the $^{17}$O/$^{18}$O ratio to be a sensitive function of initial mass and to remain constant throughout the entire TP-AGB phase for stars initially less massive than 5\,$M_{\odot}$. This makes the measured ratio a probe of the initial stellar mass. Observed $^{17}$O/$^{18}$O ratios are found to be well in the range predicted by stellar evolution models that do not consider convective overshooting. From this, accurate initial mass estimates are calculated for seven sources. For the remaining two sources two mass solutions result, though with a larger probability that the low-mass solution is the correct one. Finally, hints at a possible separation between M/S- and C-type stars when comparing the $^{17}$O/$^{18}$O ratio to the stellar pulsation period are presented.

Six billion years from now, while evolving on the asymptotic giant branch (AGB), the Sun will metamorphose from a red giant into a beautiful planetary nebula. This spectacular evolution will impact the Solar System planets, but observational confirmations of the predictions of evolution models are still elusive as no planet orbiting an AGB star has yet been discovered. The nearby AGB red giant L2 Puppis (d = 64 pc) is surrounded by an almost edge-on circumstellar dust disk. We report new observations with ALMA at very high angular resolution (18 x 15 mas) in band 7 (f ~ 350 GHz) that allow us to resolve the velocity profile of the molecular disk. We establish that the gas velocity profile is Keplerian within the central cavity of the dust disk, allowing us to derive the mass of the central star L2 Pup A, mA = 0.659 +/- 0.011 +/- 0.041 Msun (+/- 6.6%). From evolutionary models, we determine that L2 Pup A had a near-solar main sequence mass, and is therefore a close analog of the future Sun in 5 to 6 Gyr. The continuum map reveals the presence of a secondary source (B) at a radius of 2 AU contributing fB/ fA = 1.3 +/- 0.1% of the flux of the AGB star. L2 Pup B is also detected in CO emission lines at a radial velocity of vB = 12.2 +/- 1.0 km/s. The close coincidence of the center of rotation of the gaseous disk with the position of the continuum emission from the AGB star allows us to constrain the mass of the companion to mB = 12 +/- 16 MJup. L2 Pup B is most likely a planet or low mass brown dwarf with an orbital period around 5 years. Its continuum brightness and molecular emission suggest that it may be surrounded by an extended molecular atmosphere or an accretion disk. L2 Pup therefore emerges as a promising vantage point on the distant future of our Solar System.

(abreviated) We aim to study the inner-wind structure (R<250 Rstar) of the well-known red supergiant VY CMa. We analyse high spatial resolution (~0".24x0".13) ALMA Science Verification (SV) data in band 7 in which four thermal emission lines of gaseous sodium chloride (NaCl) are present at high signal-to-noise ratio. For the first time, the NaCl emission in the inner wind region of VY CMa is spatially resolved. The ALMA observations reveal the contribution of up to four different spatial regions. The NaCl emission pattern is different compared to the dust continuum and TiO2 emission already analysed from the ALMA SV data. The emission can be reconciled with an axisymmetric geometry, where the lower density polar/rotation axis has a position angle of ~50 degrees measured from north to east. However, this picture can not capture the full morphological diversity, and discrete mass ejection events need to be invoked to explain localized higher-density regions. The velocity traced by the gaseous NaCl line profiles is significantly lower than the average wind terminal velocity, and much slower than some of the fastest mass ejections, signalling a wide range of characteristic speeds for the mass loss. Gaseous NaCl is detected far beyond the main dust condensation region. Realising the refractory nature of this metal halide, this hints at a chemical process preventing all NaCl from condensing onto dust grains. We show that in the case of the ratio of the surface binding temperature to the grain temperature being ~50, only some 10% of NaCl remains in gaseous form, while for lower values of this ratio thermal desorption efficiently evaporates NaCl. Photodesorption by stellar photons seems not to be a viable explanation for the detection of gaseous NaCl at 220 Rstar from the central star, and instead, we propose shock-induced sputtering driven by localized mass ejection events as alternative.

Recent high-resolution observations have shown stellar winds to harbour complexities which strongly deviate from spherical symmetry, generally assumed as standard wind model. One such morphology is the archimedean spiral, generally believed to be formed by binary interactions, which has been directly observed in multiple sources. We seek to investigate the manifestation in the observables of spiral structures embedded in the spherical outflows of cool stars. We aim to provide an intuitive bedrock with which upcoming ALMA data can be compared and interpreted. By means of an extended parameter study, we model rotational CO emission from the stellar outflow of asymptotic giant branch stars. To this end, we develop a simplified analytical parametrised description of a 3D spiral structure. This model is embedded into a spherical wind, and fed into the 3D radiative transfer code LIME, which produces 3D intensity maps throughout velocity space. Subsequently, we investigate the spectral signature of rotational transitions of CO of the models, as well as the spatial aspect of this emission by means of wide-slit PV diagrams. Additionally, the potential for misinterpretation of the 3D data in a 1D context is quantified. Finally, we simulate ALMA observations to explore the impact of interefrometric noise and artifacts on the emission signatures. The spectral signatures of the CO rotational transition v=0 J=3-2 are very efficient at concealing the dual nature of the outflow. Only a select few parameter combinations allow for the spectral lines to disclose the presence of the spiral structure. The inability to disentangle the spiral from the spherical signal can result in an incorrect interpretation in a 1D context. Consequently, erroneous mass loss rates would be calculated...

Gamma Doradus stars (hereafter gamma Dor stars) are gravity-mode pulsators of spectral type A or F. Such modes probe the deep stellar interior, offering a detailed fingerprint of their structure. Four-year high-precision space-based Kepler photometry of gamma Dor stars has become available, allowing us to study these stars with unprecedented detail. We selected, analysed, and characterized a sample of 67 gamma Dor stars for which we have Kepler observations available. For all the targets in the sample we assembled high-resolution spectroscopy to confirm their F-type nature. We found fourteen binaries, among which four single-lined binaries, five double-lined binaries, two triple systems and three binaries with no detected radial velocity variations. We estimated the orbital parameters whenever possible. For the single stars and the single-lined binaries, fundamental parameter values were determined from spectroscopy. We searched for period spacing patterns in the photometric data and identified this diagnostic for 50 of the stars in the sample, 46 of which are single stars or single-lined binaries. We found a strong correlation between the spectroscopic vsini and the period spacing values, confirming the influence of rotation on gamma Dor-type pulsations as predicted by theory. We also found relations between the dominant g-mode frequency, the longest pulsation period detected in series of prograde modes, vsini, and log Teff.