{'fraction_non_alphanumeric': 0.04286669983017994, 'fraction_numerical': 0.034341389779920066, 'mean_word_length': 4.259179070816419, 'pattern_counts': {'":': 0, '<': 6, '<?xml version=': 0, '>': 12, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
NuSTAR and Suzaku X-ray Spectroscopy of NGC 4151: Evidence for Reflection from the Inner Accretion Disk

We present X-ray timing and spectral analyses of simultaneous 150 ks Nuclear Spectroscopic Telescope Array (NuSTAR) and Suzaku X-ray observations of the Seyfert 1.5 galaxy NGC 4151. We disentangle the continuum emission, absorption, and reflection properties of the active galactic nucleus (AGN) by applying inner accretion disk reflection and absorption-dominated models. With a time-averaged spectral analysis, we find strong evidence for relativistic reflection from the inner accretion disk. We find that relativistic emission arises from a highly ionized inner accretion disk with a steep emissivity profile, which suggests an intense, compact illuminating source. We find a preliminary, near-maximal black hole spin a>0.9 accounting for statistical and systematic modeling errors. We find a relatively moderate reflection fraction with respect to predictions for the lamp post geometry, in which the illuminating corona is modeled as a point source. Through a time-resolved spectral analysis, we find that modest coronal and inner disk reflection flux variation drives the spectral variability during the observations. We discuss various physical scenarios for the inner disk reflection model, and we find that a compact corona is consistent with the observed features.


INTRODUCTION
The dominant features of the X-ray spectra of Seyfert galaxies are a Comptonized continuum, reflection from the accretion disk as well as distant material, and line of sight absorption. Uniquely disentangling these spectral features has been challenging with most Xray observatories, which primarily have high sensitivity in the E ∼ 0.1 − 10 keV band. High sensitivity, broadband (0.5 − 79 keV) observations of Seyfert galaxies enabled by the Nuclear Spectroscopic Telescope Array (NuSTAR) ) combined with observatories with high sensitivity and energy resolution at lower energies, including Suzaku (Mitsuda et al. 2007), allow the continuum, reflection, and absorption features to be definitively separated in Seyfert galaxy X-ray spectra. Deconvolving these features can allow for fundamental properties of active galactic nucleus (AGN) coronae, supermassive black holes (SMBHs), and AGN absorption structures to be determined (e.g., NGC 1365, Risaliti et al. 2013, Walton et al. 2014; SWIFT J2127.4+5654, Marinucci et al. 2014a).
The dominant X-ray continuum emission in Seyfert galaxies is commonly thought to be produced by Comptonization of black hole accretion disk thermal emission in a coronal electron plasma (e.g. Haardt et al. 1994, Reynolds & Nowak 2003. The corona could be produced by an accretion disk atmosphere or the base of a jet (e.g., Markoff et al. 2005). Seyfert galaxies can produce weak jets even though they are radio quiet (Ghisellini et al. 2004). For some AGNs, the corona is compact with a characteristic distance from the accretion disk D 10 r g (where r g ≡ GM/c 2 ) as inferred from X-ray micro-lensing, reverberation, and eclipse measurements (e.g., Chartas et al. 2009;Zoghbi et al. 2012;Reis & Miller 2013;Sanfrutos et al. 2013 and references therein). Because of the compact sizes inferred for AGN coronae, their detailed geometry has been difficult to determine. Coronal emission is reprocessed in the black hole accretion disk, which primarily produces relativistically skewed fluorescent Fe Kα line emission (rest-frame energy E = 6.4 keV for neutral Fe) and a Compton scattering 'hump' at energies E 10 keV (e.g., Fabian et al. 1989, George & Fabian 1991. Through measuring the relativistic distortion of the inner accretion disk reflection, important properties of the inner accretion disk and black hole can be determined including the black hole spin (Brenneman 2013, Reynolds 2013, and references therein). SMBH spin (which has dimensionless form a ≡ cJ/GM 2 where J and M are the black hole angular momentum and mass, respectively) encodes the accretion and merger history of a SMBH (Berti & Volonteri 2008). Building upon the census of the ∼ 22 SMBH spin measurements currently in the literature (e.g., Brenneman 2013, Reynolds 2013, Walton et al. 2013, and references therein) is crucial for understanding SMBH evolution. SMBH spin also likely plays a critical role in driving AGN outflows and jets that deposit matter and energy in their environment (e.g., Fabian 2012).
Nuclear emission is absorbed in several regions commonly associated with ionized outflows (e.g., Elvis 2000), the broad-line region (Antonucci 1993;Urry & Padovani 1995), and/or a molecular torus (e.g, Krolik & Begelman 1988). Absorption and relativistic reflection can produce similar spectral features in the Fe K band. This has led some authors to suggest that the relativistic reflection features actually arise purely from absorption and Compton scattering in material relatively far from the black hole where relativistic effects are negligible (e.g., Miller et al. 2009).
The Seyfert 1.5 galaxy NGC 4151 (z=0. 00332, de Vaucouleurs et al. 1991), which is sometimes considered to be the archetypical Seyfert 1 galaxy, is a source in which many AGN phenomena were first characterized (for a review, see Ulrich 2000). The source has bolometric luminosity L bol ∼ 5 × 10 43 erg s −1 (Woo & Urry 2002), SMBH mass M = 4.57 +0.57 −0.47 × 10 7 M ⊙ from optical and UV reverberation (Bentz et al. 2006), and corresponding Eddington ratio L bol /L Edd ∼ 0.01. Despite being one of the brightest and most studied AGN, few constraints on the coronal geometry and the black hole spin have been found for the source.
NGC 4151 has a harder X-ray continuum relative to the average Seyfert AGN. If approximated with a powerlaw, it has a photon index historically measured to be in the range Γ ∼ 1.3−1.9 (e.g., Ives et al. 1976;Perola et al. 1986;Yaqoob & Warwick 1991;Zdziarski et al. 1996;Piro et al. 1999;Beckmann et al. 2005). NGC 4151 was also the first AGN in which a high-energy cutoff was clearly detected (Jourdain et al. 1992, SIGMA;Maisack et al. 1993, CGRO OSSE). From an analysis of all INTEGRAL data of the source and data from other X-ray observatories during the period 2003 January to 2009 June, Lubiński et al. (2010) found that the coronal emission had an approximately constant X-ray spectral index and Compton y parameter and inferred that the corona had an approximately constant geometry. The X-ray continuum conceivably could arise in the base of a jet, as NGC 4151 has displayed a non-relativistic jet as observed using radio interferometry (Wilson & Ulvestad 1982;Mundell et al. 2003;Ulvestad et al. 2005) and characterized with optical (e.g., Storchi-Bergmann et al. 2009) and X-ray observations (e.g., Wang et al. 2011a, Chandra ACIS).
NGC 4151 has shown inconsistent evidence for reflection from the inner accretion disk from X-ray spectroscopy. Spectral analyses of observations taken with instruments with a E ∼ 0.1 − 10 keV bandpass and CCD resolution have indicated the presence of relativistic reflection in NGC 4151 (e.g., Yaqoob et al. 1995, ASCA;Nandra et al. 2007, XMM-Newton). However, other analyses find no such evidence for it (e.g., Schurch et al. 2003, XMM-Newton;Patrick et al. 2012, Suzaku+Swift ). The lack of simultaneous broadband (E ∼ 1 − 100 keV), high-sensitivity X-ray observations of NGC 4151 combined with the complex absorption (total N H ∼ 10 23 cm −2 ) and soft emission (e.g., Armentrout et al. 2007, Wang et al. 2011b in the source have also made the unambiguous detection of relativistic reflection challenging. Recent X-ray reverberation studies of eight archival XMM-Newton observations have shown that NGC 4151 has a compact corona surrounding a near-maximally rotating black hole. The discovery of Fe Kα reverberation in NGC 4151 by Zoghbi et al. (2012) indicates the presence of broad Fe Kα line emission that responds to coronal emission originating from a height above the accretion disk on the order of a few r g . A more recent Fe Kα reverberation analysis by Cackett et al. (2014) shows evidence supporting a maximally spinning (a = 0.998) black hole compared to a low spin black hole (a = 0.1) and an illuminating source characterized by the lamp post geometry (in which the corona is modeled as a point source on the spin axis; Martocchia et al. 2002, Miniutti et al. 2003 with height h = 7.0 +2.9 −2.6 r g (1σ errors are given). On larger scales, NGC 4151 exhibits partial-covering, neutral material (e.g., Holt et al. 1980) and ionized outflows (e.g., Wang et al. 2011b and references therein) viewed in absorption. The source has displayed absorption column density variability on time-scales as short as ∆t ∼ 2 days that indicates that the absorbing material is located in the broad-line region (Puccetti et al. 2007, de Rosa et al. 2007). The absorption has also been characterized as having an ionized 18 (ξ ∼ 500 erg cm s −1 ), outflowing (v out ∼ 500 km s −1 ) component suggestive of an accretion disk wind (Chandra HETG; Kraemer et al. 2005) as well as a highly ionized (ξ ∼ 10 4 − 10 5 erg cm s −1 ) and blue-shifted (v out ∼ 0.1c) 'ultra-fast outflow' (UFO) component that might be associated with an accretion disk wind or a weak jet (Tombesi et al. 2011).
In this paper, we present timing and spectral analyses of simultaneous 150 ks NuSTAR and Suzaku X-ray observations of NGC 4151 in a moderately bright state focused on deconvolving the continuum, reflection, and absorption features. In a separate paper, we will present a reverberation analysis of these observations (Zoghbi et al., in prep.). This paper is structured as follows: §2 summarizes the observations and data reduction; §3  Table 1. Source and background spectra provide a source-dominated view of the 2.5 − 79 keV emission from NGC 4151 as shown in Figure 1. We reprocessed all data using HEAsoft version 6.16, which includes XSPEC v12.8.2 (Arnaud 1996) and xronos v5.22.

Suzaku
We produced Suzaku data from the XIS (Koyama et al. 2007) (XIS 0, XIS 1, and XIS 3) using the aepipeline script following the Suzaku ABC Guide 19 with version 22 of the Suzaku software and the XIS CALDB release of 2015 January 5, which includes updated XIS gain files. We extracted Suzaku XIS data products (i.e. spectra, response files, and light curves) using xselect and xisresp from circular, 170 ′′ radius regions centered on the source. For each observation, background spectra were extracted from as much of the chip as possible. Source and background regions excluded all contaminating sources in the field of view (primarily the BL Lac object 1E1207.9+3945 and the LINER NGC 4156) and calibration sources in the corners of the chip.
We rebinned Suzaku XIS spectra and response files from 4096 to 2048 channels to expedite spectral fitting without compromising the energy resolution of the data. We then co-added data from the XIS frontilluminated (XIS-FI) instruments, XIS 0 and XIS 3, using the addascaspec script to maximize signal-to-noise (S/N). We rebinned XIS source spectra into channels with a minimum of 25 counts per bin using grppha to allow for robust χ 2 fitting. In the spectral analysis, we 19 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ apply a cross-normalization constant to all detectors relative to the co-added XIS-FI data. XIS-1 fits with a crossnormalization of 0.96±0.01, which is lower than the nominal XIS-1/XIS-FI cross-normalization 20 , 1.019 ± 0.010. However, our value is reasonable given that the crossnormalization depends on the exact spectral shape of the source and exactly where the source is on the chip (K. Mukai, private communication).
Because Suzaku and NuSTAR observed NGC 4151 in a relatively bright state, we assessed the influence of pileup in the XIS data (note that NuSTAR does not suffer from pileup) using the pileest script. We found that, for XIS 0, 1, and 3, the central ∼ 30 ′′ has a pileup fraction 3% f pl 10% assuming a grade migration parameter of 0.5. We found that this minor pileup fraction does not influence our conclusions in the timing and spectral analyses on the spectral properties of NGC 4151 as described in Appendix A.
We reprocessed the Suzaku Hard X-ray Detector (HXD) PIN instrument (Takahashi et al. 2007) data using the aepipeline script with the 2011 September 15 HXD calibration files. We reduced and extracted it following the Suzaku ABC Guide 19 . We estimated the PIN non-X-ray background using the version 2.2 'tuned' non-X-ray background (NXB) event file for November 2012 as well as response and flat-field files for epoch 11 data. We estimated the PIN cosmic X-ray background with XSPEC simulations using the model of Boldt (1987). We rebinned PIN data to a S/N of 5 using ISIS (Houck & Denicola 2000) to facilitate χ 2 fitting. Three percent systematic errors were added to PIN data as recommended by the Suzaku team to account for uncertainties present in the NXB data. We allow the crossnormalization factor of PIN relative to the XIS-FI instruments to vary in the spectral analysis to improve upon the agreement between the PIN and FPMA/FPMB spectra found with the nominal value 20 , 1.164 ± 0.014. It fits to 1.26 ± 0.01.

NuSTAR
We reprocessed NuSTAR focal plane module A and B (FPMA and FPMB) data using the NuSTAR Data Analysis Software (NuSTARDAS) version 1.4.1, using the NuSTAR CALDB from 2014 July 1. NuSTARDAS performs standard data calibration and data screening, including the elimination of data collected when the Earth occulted NGC 4151 and when NuSTAR was in the South Atlantic Anomaly. For each of the NuSTAR observations, we extracted source data products from a 70 ′′ radius circular region centered on NGC 4151, and we extracted background spectra from 70 ′′ radius circular regions on the same detector chosen to avoid contamination and detector edges. Background regions were separated by at least 4 ′ from NGC 4151 and 7 ′ from 1E1207.9+3945 and NGC 4156, which exclude 98% and 99% of the enclosed counts for a point source, respectively . Spectra from the three NuSTAR observations were co-added and analyzed separately for FPMA and FPMB. We rebinned FPMA and FPMB spectra into channels with a minimum S/N of 5 to allow for χ 2 fitting. We allow the FPMA and FPMB cross-normalization constants relative to XIS-FI to fit freely in the spectral analysis, and they fit to 1.04 ± 0.01 and 1.08 ± 0.01, respectively.

TIMING ANALYSIS
The Suzaku XIS and NuSTAR light curves and hardness ratios have small but significant variability as shown in Figure 2. Notably, a drop in flux coincident with an increase in spectral hardness is seen during a ∆t ∼ 15 ks period labelled interval 3. During interval 3, the 3−5 keV flux drops by a factor of ∼ 1.4 and the 10 − 30 keV flux decreases by a factor of ∼ 1.2. At the same time, the hardness ratio (HR) between the 6 − 10 keV flux and the 3 − 5 keV flux increases by ∼ 5σ. Between intervals 4 and 6 (∆t ∼ 45 ks), the 3 − 5 keV flux rises by a factor of ∼ 1.6, the 10 − 30 keV flux increases by a factor of ∼ 1.3, and the F (6 − 10 keV)/F (3 − 5 keV) HR decreases by ∼ 4σ. From the beginning of interval 6 to the end of interval 7 (∆t ∼ 40 ks), the 3−5 keV flux drops by a factor of ∼ 1.4, the 30 − 79 keV flux decreases by a factor of ∼ 1.3, and the F (6 − 10 keV)/F (3 − 5 keV) HR remains approximately constant.
The root-mean-square fractional variability (RMS F var ) is shown in Figure 3. The overall source variability drops dramatically below E ∼ 3 keV, where circumnuclear and scattered emission dominate the absorbed emission from the AGN (e.g., Wang et al. 2011b). There is also an overall drop in variability with increasing energy above E ∼ 3 keV, which is typical of Seyfert AGN, where the more variable model components tend to dominate at softer energies. The sharp, relatively narrow dip at 6.4 keV is an indication that the neutral, narrow component of the Fe Kα line is not varying much relative to the surrounding continuum during the observation. This is expected of the neutral, narrow reflection component because it likely originates relatively far from the nucleus.

SPECTRAL ANALYSIS
We perform time-averaged and time-resolved analyses to characterize the broadband X-ray spectral features and variability of NGC 4151. We first undertake a time-averaged spectral analysis (i.e., jointly fitting the full XIS-FI, XIS-1, PIN, FPMA, and FPMB spectra) to understand the components present in the source spectrum using the high signal-to-noise (S/N) of the full data. We then carry out a time-resolved analysis to characterize spectral variability that is averaged out in the full observations, which may include variations in the coronal emission (which can vary on hour-long time scales in AGN; e.g., Markoff et al. 2005), the inner accretion disk reflection emission (which can vary on the order of hours in NGC 4151; e.g., Zoghbi et al. 2012) and absorption features (which can vary on the order of days in NGC 4151; e.g., Puccetti et al. 2007).
To carry out the spectral analyses, we use XSPEC with a standard cosmology of H 0 = 70 km s −1 Mpc −1 , Ω M = 0.27, and Ω Λ = 0.73 (Komatsu et al. 2011). We use X-ray cross-section values from Verner et al. (1996) and solar abundance values from Wilms et al. (2000). We include data in the energy ranges indicated in Table 1.
For the spectral analysis, we apply an inner disk reflection model and an absorption-dominated model. In brief, the inner disk reflection model includes a partialcovering, neutral absorption component (which we refer to as PC1), a distant reflection component (DRC), a cut-off power-law component (PLC), and an inner disk reflection component (IDR). The absorption-dominated model is identical to the inner disk reflection model except that the IDR component is replaced with another partial-covering absorber (PC2). For the time-averaged analysis, we link all parameters between data sets, and we fix redshifts of all model components to the systemic redshift. We model Milky way absorption with Tbabs (Wilms et al. 2000), and we fix the Milky Way absorption column density to the weighted-average value for NGC 4151 from the LAB survey (Kalberla et al. 2005), N H = 2.3 × 10 20 cm −2 . We determine 90% confidence intervals for one interesting parameter for all parameter values presented in the time-averaged and time-resolved analyses using the XSPEC error algorithm unless noted otherwise.

Time-Averaged Analysis
To develop the time-averaged model, we first characterize the continuum absorption and emission features in the 2.5 − 4.0 keV, 7.5 − 10 keV, and 50 − 79 keV bands. We model the coronal continuum emission with a cutoff power-law with photon index Γ and cut-off fixed at E cut = 1000 keV. Data/model ratios in Figure 4 illustrate the presence of a strong absorption cut-off at energy E ∼ 2 − 3 keV indicative of a relatively high lineof-sight column density (N H ∼ 10 23 cm −2 ) as well as a prominent, narrow Fe Kα emission line at 6.4 keV and Compton hump indicative of reflection.
We apply a partial-covering, neutral absorption component (partcov*zTbabs; Wilms et al. 2000). This absorber greatly improves the fit (see Table 2) and satisfactorily accounts for the low-energy continuum ratios as shown in Figure 4. We then model the prominent, narrow Fe Kα emission line and Compton hump as arising by reflection of cut-off power-law emission from Compton-thick, plane-parallel, and relatively neutral (i.e., ξ = 1 erg cm s −1 ) material with the xillver  Note. -Each added model component provides a significant improvement to the fit. a Values includes a change ∆χ 2 /∆ν = −39/ − 4 from allowing the instrument cross-normalizations to fit freely. reflection model (García et al. 2013). We specifically employ xillver-Ec3 21 , which allows for a variable cut-off energy in the range E = 20 − 1000 keV. We fix the incident cut-off power-law photon index and cut-off energy of this component self-consistently to those of the cut-off power-law component, and we additionally fix its inclination to i = 60 deg. As shown in the bottom panels of Figure 4, the added xillver component significantly reduces the ratios towards unity. The remaining residuals indicate the presence of blended Fe XXV and Fe XXVI Kα (E = 6.67 keV and E = 6.97 keV, respectively) and Kβ (E = 7.88 keV and E = 8.25 keV, respectively) absorption lines as well as broad emission in the Fe K band and convex curvature above E ∼ 10 keV.
To account for the unmodeled absorption lines in the Fe band, we additionally apply two separate XSTAR absorption table models, which we refer to as warm absorbers one (WA1) and two (WA2). WA1 is parameterized by its column density in the range N H = 10 22 −10 24 cm −2 and ionization parameter ξ in the range log(ξ/ erg cm s −1 ) = 2.5 − 4.5. The absorber has turbu-21 available at http://hea-www.cfa.harvard.edu/~javier/xillver/ lent velocity v t = 3000 km s −1 , which is appropriate for a highly ionized absorber (Tombesi et al. 2011). We fix the Fe abundance to one. WA2 is parameterized by its column density in the range N H = 10 20 − 10 24 cm −2 and ξ in the range log(ξ/ erg cm s −1 ) = 0.0 − 4.0. It has a turbulent velocity v t = 200 km s −1 . The Fe abundance is hard-wired at one. We limit ξ to be in the range log(ξ/ erg cm s −1 ) = 2.5 − 4.0. For both WA1 and WA2, ξ is sampled with resolution ∆log(ξ/ erg cm s −1 ) = 0.2, which is important for accurately modeling the absorption opacities and, consequently, determining the black hole spin . Having accounted for these spectral features, we develop separate inner disk reflection and absorption-dominated models to account for the remaining broad residuals.

Inner disk reflection model
We first characterize the broad residuals in the Fe K-band phenomenologically. We begin by applying a broad Gaussian. The line fits with equivalent width EW = 66 +16 −15 eV, σ = 0.39 ± 0.04 keV, and line energy E = 6.20 ± 0.05 keV. The fit improves by ∆χ 2 /∆ν = −341/ − 3. The measured line energy and width are suggestive of broadened Fe Kα emission. To test for an origin from the inner accretion disk, we replace the Gaussian with a model of line emission from the inner accretion disk (relline; Dauser et al. 2010) with an emissivity profile ǫ ∝ r −q with q = 3 fixed. The relline component fits with line energy E = 6.52 ± 0.01 keV, EW = 111 +16 −22 eV, inclination i = 6 ± 1 deg, and spin a > 0.99. The addition of the broad Fe line is statistically motivated by a reduction in χ 2 by ∆χ 2 /∆ν = −354/ − 4.
To account for emission from the inner accretion disk self-consistently (i.e., by including the full reflection spectrum, including the Compton hump, Fe Kα and Kβ emission lines, and the Fe K absorption edge), we include an inner disk reflection component (IDR) composed of the xillver model convolved with the relconv relativistic convolution model (Dauser et al. 2010). Specifically, we use an optimized version relconvf (Jeremy Sanders, private communication). relconv allows the black hole spin, a, to fit freely in the range −0.998 ≤ a ≤ 0.998. We further allow the inner disk reflection component to have a broken power-law emissivity profile that transitions from inner emissivity index q 1 to outer emissivity index q 2 at a radius r br . relconv additionally assumes an inner accretion disk that is geometrically thin, optically thick, and aligned perpendicular to the spin axis.
We assume that the inner disk radius (r in ) is fixed at the innermost stable circular orbit (ISCO) radius and that the outer disk radius is fixed at r out = 400 r g . This value should be unconstrained assuming q 2 3 as expected for rapidly spinning black holes (e.g., Wilkins & Fabian 2011). We assume an isotropic emission distribution as discussed in Svoboda et al. (2009). We further assume that the AGN is chemically homogeneous as discussed in, e.g., Reynolds et al. (2012) and Walton et al. (2013), so we fix the Fe abundances of the distant reflector and IDR together. We additionally limit the Fe abundance to be in the range 0.5 < A Fe < 5.0 in order to keep the Fe abundance within the range measured for relatively neutral gas in the source (Yaqoob et al. 1993;Schurch et al. 2003; Nandra et al.  2007). We are also motivated to do this to limit the influence of the degeneracy between the black hole spin and the Fe abundance . As shown in Table 2, the inclusion of the inner disk reflection component is motivated by a reduction in χ 2 /ν by ∆χ 2 /∆ν = −1060/ − 7.
Best-fit parameter values are shown in Table 3, and the best-fit model and data/model ratios are shown in Figure 5.
Applying an inner disk reflection model, we find strong evidence for relativistic reflection and a preliminary SMBH spin constraint a = 0.98 ± 0.01 as shown in Figure 6. If we calculate the reflection fraction, R, as the ratio of the inner disk reflector flux to the power-law component flux in the 20 − 40 keV band, we find that R = 1.3 ± 0.2. This reflection fraction is relatively low given the near-maximal spin and relatively steep inner emissivity index, q 1 = 10.0 +0.0p −0.4 , with respect to the maximum expected reflection fraction in the lamp post geometry . We discuss this further in §5.
We find a partial covering fraction that is consistent with a scattered fraction. Parameters largely do not show degeneracies as illustrated with confidence intervals in Figure 7. We do find a degeneracy between the inner disk component spin and inclination although we find values for these parameters that are consistent with previous studies as discussed in §5. For the cut-off power-law component, we find only a lower limit E cut > 600 keV if we allow the cut-off energy to fit freely. For simplicity, we fix E cut = 1000 keV in the best-fit model and find Γ = 1.75 +0.01 −0.02 . To illustrate the spectral features that are modeled by the inner disk reflection component, we remove this component from the best-fit model giving the data/model ratios shown in Figure 8. We then refit this model, which gives a significantly worse fit than the best-fit model (∆χ 2 /∆ν = +1060/ + 7). This increase in χ 2 is driven by residuals in the Fe K-band and the energy range E 10 keV identical to those motivating its inclusion.
We note that non-unity ratios remain at the < 10% level in the best-fit data/model ratios at energies E 30 keV, as χ 2 /ν = 1304/1219 (1.07) ignoring data below 30 keV. We attempt to account for these non-unity ratios with a phenomenological hard power-law component, which have been observed in radio-loud AGNs (e.g., 3C 273; Grandi & Palumbo 2004), in order to test whether this may be emission from the dim, extended radio jet in the source (Mundell et al. 2003;Wang et al. 2011a). The hard power-law component fits with Γ hard = 1.23 +0.05 −0.10 , and the coronal power-law Γ increases to Γ = 1.83 +0.03 −0.02 . This component gives a small but significant reduction in the fit statistic of ∆χ 2 /∆ν = −39/ − 2. The presence of such a hard power-law component, which could perhaps correspond to extended jet emission, would be interesting although we caution that the component is mainly significant in an energy band where there is more uncertainty in the NuSTAR calibration.
In attempt to constrain the outer extent of the corona, we try replacing the relconv component in the best-fit model with a relativistic convolution component with a twice broken power-law emissivity profile, kdblur3 (Wilkins & Fabian 2011). We allow the inner disk radius to fit freely because the spin is fixed at a = 0.998. The outer index is fixed to q 3 = 3 as expected from basic geometry at radii in the disk where r ≫ h. We limit the outer break radius to values r br,2 > 10 r g , and we find r br,2 = 10 +6 −0p r g . We also find r in = 1.5 ± 0.1 r g , q 1 = 10.0 +0.0p −0.4 , r br,1 = 3.5±0.1 r g , and q 2 = 2.5 +0.3 −0.2 . This model provides only a slightly better fit than the best-fit model (∆χ 2 /∆ν = −4/−1), so a twice-broken power-law emissivity profile is not strongly warranted statistically.
If we apply a lamp post model, relconv lp , in place of relconv in the bestfit model, we find a lamp post height h = 1.3 +0.1 −0.0p r g and a reflection fraction consistent with the best-fit model (R = 1.2 ± 0.1) although we also find a significantly worse fit (∆χ 2 /∆ν = +493/ + 2 relative to the bestfit model). We note that the worse fit with relconv lp results from the inner disk emissivity profile being determined solely with the lamp post height compared to the broken power-law profile used in relconv.
We attempt to add a second lamp post component by including a second inner disk component because coronae could be extended, such as in a scenario where a jet produces the coronal emission or a corona has multisite activity. We fix all of its parameters to those of the first inner disk component except for the normalization and lamp post height. We find a similar reflection fraction and an improved fit relative to the single lamp post model (∆χ 2 /∆ν = −240/ − 2) where the second lamp post component has height h = 17 ± 3 r g and low flux relative to the first lamp post component, which again fits with a low height h = 1.3 +0.3 −0.0p r g . This fit along with our best-fit model indicate that a model more complicated than a single point source on the spin axis provides a statistically favored description of the reflection spectrum. This fit also suggests that the corona may have an extended component.
In attempt to test the extended corona scenario self-consistently, we apply the relconv LP ext model . We model the disk-illuminating source as an outflow on the spin axis with outflow velocity v out , base height h base , and maximum height h top . We fix h base to the event horizon radius as motivated by the dual lamp post model described above. We find v out /c = 0.01 +0.08 −0.01p and h top = 1.3 +0.1 −0.0p r g . This model (which gives ∆χ 2 /∆ν = +477/ + 1 relative to the bestfit model) is consistent with the single lamp post model presented above within errors. We find a better fit with the dual lamp post model than with this extended corona model, which may result from describing the corona with two physically discrete illuminating components rather than a continuous illuminating source.
We try applying the relxill model in place of the IDR component in the best-fit model to more accurately account for the emission angle as a function of radius. The inner disk reflection models we employ in the best-fit model assume a single emission angle for the inner accretion disk emission and also assume an isotropic emission distribution. While these are valid assumptions for constraining properties of relativistic reflection, the emis-  sion angle of photons leaving the disk is not constant and depends upon radius, especially in the innermost regions of the disk as a consequence of relativistic lightbending (García et al. 2014). We find parameter values consistent within errors with those for the best-fit model including a reflection fraction in the 20 − 40 keV band R = 1.3 ± 0.2 and a similar fit to the best-fit model (∆χ 2 /∆ν = −5/ + 0). Consequently, the assumption of a single emission angle for the inner accretion disk does not introduce significant systematic error into our results.
We find a best-fit model shown in Table 4 and Figure 9. This model provides a statistically less satisfactory description of the time-averaged data (∆χ 2 /∆ν = +86/ + 5) relative to the best-fit model.
The absorption-dominated model also has positive residuals at the < 10% level at energies E > 30 keV, so we also attempt adding a hard power-law component to it. We find a similar hard power-law component as for the inner disk reflection model with Γ hard = 1.25 +0.08 −0.07 can account for the high energy non-unity ratios. We also find a moderate reduction in the fit statistic of ∆χ 2 /∆ν = −24/ − 2.

Time-Resolved Analysis
To perform the time-resolved analysis, we produce strictly simultaneous time-resolved data products (i.e., the data products are time-filtered to only include times when the two observatories were simultaneously viewing NGC 4151) in seven time-intervals as shown in Figure 2. Intervals are chosen to maximize S/N and have approximately constant hardness ratios. We omit making PIN time-resolved products due to NuSTAR's superior sensitivity at E > 10 keV. Information for the time-resolved data is shown in Table 5. We note that similar timeresolved analyses have been performed on joint NuSTAR observations with Suzaku or XMM-Newton of the Seyfert galaxies IC 4329A (Brenneman et al. 2014), NGC 1365 (Walton et al. 2014), and MCG 6-30-15 (Marinucci et al. 2014b).
To first assess the spectral variability in a model independent fashion, we analyze difference spectra between high flux states (intervals 2 and 6) and a low flux state (interval 4) as shown in Figure 10. The difference spectra peak at E ∼ 3 keV and gradually decrease with increasing energy. Additionally, the flux of the narrow, neutral Fe Kα line at 6.4 keV varies less than the continuum.
To fully assess the spectral variability, we perform a joint fit of the time-resolved spectra using a strategy similar to our time-averaged analysis. We initially set the parameter values to their time-averaged model values (which are shown in Tables 3 and 4). We require many parameters to fit freely to the same value across all intervals on physical grounds. We assume the instrument cross-normalizations are constant during the observations. The distant reflection component, which likely arises from material at distances in the range r ∼ 7 lightdays (which corresponds to the approximate distance to the optical broad-line region in the source; Bentz et al. 2006) to r ∼ 4.2×10 4 light-days (which is the upper-limit on the size of a dusty torus in the source; Radomski et al. 2003), likely does not vary on the time-scale of the observations as also suggested by the RMS F var of the data. The Fe abundance and inclination of the inner accretion disk also can be assumed to be constant during the observations. Similarly, the black hole spin can be safely  assumed to be constant, as SMBH angular momentum evolves only on cosmic time-scales from accretion and/or mergers (e.g., Berti & Volonteri 2008). We additionally fix the WA2 ionization parameter to its time-averaged best-fit value.

Inner disk reflection model
We systematically assess which parameters should be variable and constant in the time-resolved analysis by carrying out fits for all combinations of variable parameters for the following seven parameters: neutral absorber  Note. -Counts correspond to the 2.5−9 keV energy bandpass for XIS-FI, 2.5 − 7.5 keV for XIS-1, and 5 − 79 keV for FPMA and FPMB. The relatively low exposure times for the interval 6 data result from relatively poor overlap between the NuSTAR and Suzaku observations during this interval.
N H , f cov , DRC normalization (K DRC ), PLC Γ, PLC 20 − 40 keV flux (F P LC ), IDR ξ, and IDR 20 − 40 keV flux (F IDR ). We parameterize the normalizations of the PLC and IDR in terms of their 20 − 40 keV flux for the time-resolved analysis using the cflux model in XSPEC. We allow parameters to be variable until we no longer find an improvement in the fit based on the F-test. We find a best-fit time-resolved fit with four variable parameters (F P LC , F IDR , N H , and Γ). If we additionally try allowing f cov to be variable, it provides a marginally significant improvement to the fit according to the F-test (∆χ 2 /∆ν = −30/ − 6). We find that this parameter is correlated with the power-law normalization, which is likely unphysical whether the covering fraction represents a scattered fraction or a true covering fraction for The ratios of the interval 2 (black, top), interval 4 (red, middle), and difference (blue, bottom) spectra to the best-fit model to interval 2. Right: The ratios of the interval 6 (black, top), interval 4 (red, middle), and difference (blue, bottom) spectra to the best-fit model to interval 6. Notice the narrow dip in both difference spectra at the energy of neutral Fe Kα (6.4 keV as indicated with a dashed, light blue line). XIS-FI, XIS-1, FPMA, and FPMB spectra are shown although only NuSTAR data is shown above 9 keV and only Suzaku data is below 9 keV for clarity. Spectra have been rebinned for plotting purposes only. The solid green line indicates a data/model ratio of unity. (A color version of this figure is available in the online journal.) a neutral absorber. Thus, we elect to keep the covering fraction constant.
In summary, the best-fit time-resolved model values are shown in Tables 6 and 7. Time-series of the parameter values allowed to vary between intervals (see Figure 11 Specifically, for small lamp post heights, the reflection fraction can increase when the lamp post moves to smaller heights due to a decrease in the continuum flux that reaches the observer (Miniutti & Fabian 2004;Dauser et al. 2014).

Absorption-dominated model
For the absorption-dominated model, we carry out a similar systematic analysis as followed for the inner disk reflection model. We perform fits for all combinations of variable parameters for the following seven parameters: the column densities and covering fractions of both neutral absorbers, DRC normalization (K DRC ), PLC Γ, and the PLC 20 − 40 keV flux (F P LC ). We allow parameters to be variable until we no longer find an improvement in the fit based on the F-test. We find a best-fit timeresolved fit with four variable parameters (F P LC , Γ, and N H and f cov of PC2) as shown in Figure 12. For the time-resolved analysis, we allow these four parameters to fit independently in each interval (which are shown in Table 8) and allow all other free parameters (which are shown in Table 9) to fit freely but linked between each interval.
We find a best-fit absorption-dominated model that provides a slightly worse fit of the time-resolved data (∆χ 2 /∆ν = +29/+5) relative to the inner disk reflection model. Within the context of the absorption-dominated model, variation in the covering fraction of PC2 and the PLC flux drives the variability during the observations. The PC2 N H and f cov parameters appear to be correlated with variations in the coronal continuum on the short timescales seen in the top panel of Figure 12. We show correlation plots in the bottom panel of Figure 12 to show that f cov and F P LC have closed confidence contours.

DISCUSSION
Using timing and spectral analyses, we have characterized the constant and variable spectral behavior of NGC 4151 as observed simultaneously with NuSTAR and Suzaku. Through a spectral analysis, we apply inner accretion disk reflection and absorption-dominated models. We find that the inner disk reflection model provides a better fit to the time-averaged and time-resolved spectra. Within the context of the inner disk reflection model, we find evidence for relativistic reflection from the inner accretion disk and a preliminary spin measurement for the supermassive black hole.
It is unlikely that a putative neutral absorber (which are at radii of r ∼ 10 4 − 10 5 r g ; Antonucci 1993, Urry & Padovani 1995 can partially cover AGN coronae (which have been found to have size D 10 r g for several AGNs including NGC 4151; Reis & Miller 2013). Furthermore, there is no physical reason to expect a neutral absorber to be correlated with variations in the coronal continuum on the short timescales seen in §4.2.2 because the neutral absorbing gas likely is located relatively far from continuum source. A reasonable range of distances for the absorbing material might extend as close as the broad-line region (r = 6.6 +1.1 −0.8 lightdays; Bentz et al. 2006) and as far as the upper-limit on the size of a dusty torus in the source (r ∼ 4.2 × 10 4 light-days; Radomski et al. 2003). In contrast, a correlation between the inner disk reflection and the coronal  Confidence contours indicate 99% (black), 90% (red), and 68% (green) confidence intervals for two interesting parameters.  continuum is physically possible based on light crossing time arguments and has been observed. The distance between the corona and the disk inferred by Cackett et al. (2014) and the size scale of the inner disk, which are both 10 r g , correspond to a light-travel time 2 ks assuming the black hole mass M = 4.57 +0.57 −0.47 ×10 7 M ⊙ (Bentz et al. 2006). The reverberation analysis of Zoghbi et al. (2012) shows that the inner disk reflection flux can lag behind the coronal flux on time-scales ∼ 2 ks.
Thus, on statistical and physical grounds, we discuss only the inner disk reflection model. Here, we discuss assumptions that may introduce systematic errors in the preliminary black hole spin value and physical scenarios for the disk-corona system. For calculations that use the black hole mass, we use M = 4.57 +0.57 −0.47 × 10 7 M ⊙ (Bentz et al. 2006) measured from optical and UV reverberation.

Systematic Uncertainties in the Black Hole Spin
Assumptions in our model introduce systematic uncertainty into the derived confidence intervals in the dimensionless black hole spin. A primary assumption of modeling the inner disk reflection with the relconv kernel is that the accretion disk is radiatively efficient, optically thick, and geometrically thin at radii greater than or equal to the ISCO radius. These assumptions generally should be applicable to NGC 4151 considering the modest accretion rate implied by its bolometric luminosity L bol ∼ 0.01L Edd . The assumption of a geometrically thin disk breaks down where the disk thickness, t, is on the same order as the radius in the disk, r, or t/r ∼ 1. For Seyfert galaxies with Eddington ratios in the range L bol /L Edd ∼ 0.01 − 0.3, we expect disk thicknesses in the range t ∼ (0.1 − 1.0) r g for a Shakura & Sunyaev (1973) disk. Thus, this assumption may be invalid for r in 2 r g , which corresponds to spins a 0.94. Thus, considering the finite thickness of the disk close to the black hole, we estimate that this assumption relaxes our spin constraint to a > 0.94.
We assume that the inner accretion disk radius is the ISCO radius and that a negligible level of emission originates within that radius. The latter assumption is supported by the near-maximal spin and relatively low Eddington ratio of the source, which may indicate that some of the accretion flow within the ISCO radius is optically thin (Reynolds & Fabian 2008). Additionally, when we allow the inner edge of the accretion disk to fit freely, it fits to a value consistent with the ISCO radius although we restrict it to fit to values r in ≥ r ISCO . Based on MHD simulations (Reynolds & Fabian 2008), if significant emission comes from within the ISCO radius, the fit will be biased to higher spins systematically with error 2% for rapidly spinning, prograde black holes (a ≥ 0.9). Thus, we estimate this assumption adds ∼ 2% systematic uncertainty to our spin value. Including the systematic uncertainties introduced by our assumptions about the thickness and extent of the inner disk, we find a constraint a > 0.9.
Black hole spin has some degeneracy with the iron abundance of the inner disk, as both properties affect the strength of the broad Fe line (e.g., Reynolds et al. 2012). We find that NGC 4151 has a super-solar Fe abundance. There is no apparent degeneracy between our derived Fe abundance and spin values as shown in Figure 7. Fe abundance is constant on long time-scales. There is evidence for super-solar Fe abundance in material distant from the SMBH. From several Ginga measurements, Yaqoob et al. (1993) report an average value A Fe ∼ 2.5 from modeling an Fe K absorption edge. From analyzing the strength of a neutral Fe absorption edge in XMM-Newton observations, Schurch et al. (2003) report an Fe abundance A Fe ∼ 3.6, and they report A Fe ∼ 2.0 if they associate the Fe abundance of the neutral absorber with that of a Compton reflection component.
For the Fe abundance of the inner disk, Nandra et al. (2007) measure an Fe abundance in the range A Fe = 0.3 − 0.5 from analyzing three XMM-Newton observations. They assume the distant reflection component and inner disk reflection component (which are modeled with pexmon and kdblur2*pexmon, respectively) are chemically homogeneous in the same way as in the present work. Nandra et al. (2007) caution that their XMM-Newton values require high sensitivity, broad-band observations that can simultaneously constrain the Fe Kα line and Compton hump reflection features. Our analysis simultaneously constrains both features, and our derived Fe abundance ( A Fe = 5.0 +0.0p −0.1 ) is moderately high in comparison with measurements of super-solar Fe abundance in material distant from the SMBH. This may introduce some bias in our spin constraint as a result of the positive correlation between spin and the Fe abundance ). However, we conclude that our derived Fe abundance negligibly biases our spin con- Note. -90% confidence intervals are given for the variable parameters for the absorption-dominated model, which are the PC2 covering fraction, PC2 column density, power-law component flux in the 20 − 40 keV band, and power-law component photon index. straint of a > 0.9.
Black hole spin can show a degeneracy with the disk inclination if the blue wing of the broad Fe line profile is not significantly detected, as both parameters affect the width of the broad Fe line profile. A degeneracy is seen between the disk inclination and spin shown in Figure 7. The inclination of the inner disk should be constant on long time-scales, so we compare our determined inclination angle with various inclination angles previously reported for the source.
The inner disk inclination angle we find is consistent with most values determined for the inner disk. From the analysis of three XMM-Newton observations, Nandra et al. (2007) find an inclination of the inner disk of i = 17 +12 −17 deg, i = 21 +69 −21 deg, and i = 33 +1 −3 deg, respectively, where values are given for the best-fit kdblur2*pexmon model presented for each observation. From an ASCA observation, Yaqoob et al. (1995) find i = 0 +19 −0 deg using the diskline model (Fabian et al. 1989). Cackett et al. (2014) report i < 30 deg for the inner disk from Fe Kα reverberation measurements.
Comparing the inner disk inclination angle to that of more distant structures, our determined inner disk inclination angle is not consistent with the inclination of the narrow-line region, 45 ± 5 deg, as determined from Hubble Space Telescope observations (Das et al. 2005). This suggests that the inner disk and narrow-line region are misaligned. The inclination we find is consistent with the galactic disk inclination, i ∼ 21 deg, as inferred from optical photometry (Simkin 1975) and Very Large Array observations of neutral hydrogen (Pedlar et al. 1992).
It is possible that some of the broad Fe K line profile may arise from an alternate mechanism than reflection from the inner accretion disk. For example, scattering in a Compton thick wind near the disk can produce broadened Fe line profiles (Sim et al. 2008). However, the strength of this feature significantly depends upon the outflow mass loss rate. A low mass outflow rate is expected for the low mass accretion rate of the source implied by its Eddington ratio L bol /L Edd ∼ 0.01. Thus, we estimate that any contribution of scattering in a Compton thick wind near the accretion disk to the broad Fe line is relatively weak compared to the broad line profile from inner accretion disk reflection.

Physical Scenarios for the Corona and Inner
Accretion Disk Our derived high inner disk emissivity index (q 1 = 10.0 +0.0p −0.4 ) indicates that the corona that illuminates the inner disk is compact. The inner emissivity index and photon index, Γ = 1.75 +0.01 −0.02 , together indicate a lamp post height h < 10 r g assuming a lamp post geometry, as the maximum inner index for Γ ∼ 1.8, h = 10 r g , and a = 0.99 is q 1 ∼ 3 . A lamp post height h 10 r g is expected in order to produce strong reflection features from which black hole spin can be well constrained Fabian et al. 2014). Furthermore, the high ionization parameter we find for the inner disk atmosphere is expected for a close illu-minating source. We note that a high emissivity index (q 1 3) is only physically consistent with a nearmaximal black hole spin.
We calculate theoretical emissivity profiles for the lamp post model of Dauser et al. (2013) for our best-fit Γ and a values. For this calculation, we assume a minimum height in order to give a rough estimate of the highest emissivity index expected. We estimate a value of q 1 ∼ 7 if we average the emissivity from r ISCO to the emissivity break radius that we find for the best-fit model (r br ). So, we find an inner emissivity that is steeper than expected for the lamp post model of Dauser et al. (2013). This partially illustrates why the broken power-law emissivity profile for the inner disk provides a better fit to the data than the lamp post emissivity profile. However, this estimate does not necessarily indicate that our best-fit model is unphysical.
We find a relatively low reflection fraction compared to the maximum expected reflection fraction, which is R ∼ 5 − 10 for the lamp post model described in §4.1.1 with a = 0.98 ± 0.01 and h = 1.3 +0.1 −0.0p r g ). This value is a factor ∼ 4 − 8 higher than what we find. There are several possible physical scenarios for the corona and inner accretion disk that may cause this apparent discrepancy.
If the corona results from a mildly relativistic outflow rather than a static region, then the reflection fraction is expected to be reduced because of relativistic aberration (e.g., Beloborodov 1999, Malzac et al. 2001. A photon index Γ ∼ 1.75 indicates v/c ∼ 0.4 near the base of an AGN jet and a reduction in the strength of reflection, Ω/2π, by a factor of ∼ 10 using the disk inclination i = 18 deg (Beloborodov 1999). This can help explain the apparently low reflection fraction we find. We note that VLBI observations have shown that the radio jet in NGC 4151 has had a component with velocity v ≤ 0.050c at a distance r ∼ 0.16 pc from the nucleus (Ulvestad et al. 2005). Thus, if the corona corresponds to the base of a mildly relativistic jet, then a mechanism is required to decelerate the jet to non-relativistic speeds on ∼ 0.1 pc scales. Also, an outflow tends to flatten the emissivity profile . Keeping these facts in mind, an outflowing corona may help explain the low reflection fraction we find.
Truncation of the inner edge of the accretion disk at radius r > r ISCO can similarly reduce the expected reflection fraction. Lubiński et al. (2010) conclude that the relatively weak strength of reflection from the disk they measure (Ω/2π ≃ 0.3) from all INTEGRAL observations of NGC 4151 from 2003 January to 2009 June can be explained by an inner hot accretion flow surrounded by a cold disk truncated at r in ∼ 15 r g or, alternatively, a mildly relativistic coronal outflow. When we allow the inner radius of the accretion disk to fit freely using the kdblur3 model as described previously, it fits to a value consistent with the ISCO radius of a black hole with the near-maximal, prograde spin that we find for the best fit model. Thus, we do not find evidence for a truncated disk from spectral fitting of our observations. Dilution of the reflection features may also reduce the measured reflection fraction. For a close diskilluminating source suggested by the steep emissivity profile, it is plausible that the inner disk sees an intense ir-radiating flux and has a tenuous, hot atmosphere. This may have a density profile roughly similar to that of an accretion disk in hydrostatic equilibrium, for which an intense irradiating flux can lead to disk thermal instability (Nayakshin & Dove 2001;Ballantyne et al. 2001). This instability can produce a highly ionized 'skin', which could both further broaden the disk reflection features and reduce the observed reflection fraction. Ballantyne et al. (2001) showed that simulated RXTE and ASCA data of hydrostatic reflection models (e.g., Ballantyne et al. 2001) fit with constant density models (Magdziarz & Zdziarski 1995;Ross & Fabian 1993) tend to fit with systematically low reflection fraction values because of dilution. Although NuSTAR provides superior sensitivity in the 20 − 40 keV energy band, where the reflection spectrum is almost completely independent of ionization parameter, our fit is biased towards best matching the spectral shape at lower energies (E 20 keV), where the reflection spectrum is sensitive to the ionization parameter.
With the xillver models, we assume that the disk is illuminated at 45 deg. However, for a compact coronal geometry, the disk illumination may be incident on the disk at closer to grazing angles, especially within a few r g of the black hole ). In the innermost regions of the disk, this would enhance the influence of dilution and thus further reduce the reflection fraction.
In summary, both dilution from a highly ionized disk skin and a mildly relativistic outflow may contribute to the reduced reflection fraction that we find relative to predictions for the lamp post geometry. Because an outflow tends to flatten the emissivity profile, the formation of a highly ionized disk 'skin' may be the dominant cause of the reduced reflection fraction. We do not find evidence for truncation of the inner accretion disk.

Coronal Properties
We model the coronal emission with a cut-off powerlaw with E Cut = 1000 keV fixed and find Γ = 1.75 +0.01 −0.02 . This value is consistent within 3σ with the average Seyfert photon index, Γ = 1.84 ± 0.03, measured for 105 Seyfert galaxies with redshift z ≤ 0.1 with Bepposax (Dadina 2008). It is relatively hard compared to the average photon index, Γ = 1.93 ± 0.01, for 144 Seyfert galaxies in the second INTEGRAL AGN catalogue (Beckmann et al. 2009). The value we find agrees within 1σ with the correlation between the photon index and Eddington ratio for a sample of 69 radio-quiet AGNs out to z ∼ 2 (Brightman et al. 2013), which gives Γ = 1.63 ± 0.16 for L bol /L Edd = 0.01. We note that the significant dispersion in this correlation as well as uncertainty in estimating the Eddington ratio should be noted when applied to a single source.
The photon index we find is consistent with the range in photon index (Γ ∼ 1.7 − 1.86) that Lubiński et al. (2010) measure for the source in a large range of Xray flux states as observed with INTEGRAL during the period 2003 January to 2009 June. The energy cut-off we find is corroborated by the energy cut-offs found by Lubiński et al. (2010) for the source in high and medium flux states (E Cut = 264 +48 −26 keV and E Cut > 1025 keV, respectively), which serve as upper and lower bounds, respectively, to the 20 − 100 keV flux of the source in our data.

CONCLUSIONS
We present X-ray timing and spectral analyses of contemporaneous, 150 ks NuSTAR and Suzaku X-ray observations of the Seyfert 1.5 galaxy NGC 4151. We test separate inner disk reflection and absorption-dominated models in order to determine the properties of the innermost regions of the AGN. Our primary conclusions are enumerated below: 1. We find that the inner disk reflection model provides a better statistical and physical description of the data than the absorption-dominated model for both the timeaveraged and time-resolved spectral analyses.
2. Within the context of this inner disk reflection model, we find that relativistic emission originates in a highly ionized inner accretion disk with a steep inner emissivity profile, which suggests a bright, compact inner disk-illuminating source. We find a relatively moderate reflection fraction with respect to predictions for the lamp post geometry. We find a preliminary, nearmaximal black hole spin a > 0.9 accounting for statistical and systematic modeling errors.
3. Through a time-resolved spectral analysis, we find that modest coronal and inner disk reflection variation drive the spectral variability during the observations. 4. We discuss physical scenarios with respect to our derived inner accretion disk and coronal properties. We find that a compact coronal geometry can reproduce the observed features.
This work further displays the ability of NuSTAR in conjunction with X-ray observatories such as Suzaku, XMM-Newton, Chandra, and Swift to robustly constrain coronal and accretion disk properties in Seyfert galaxies.
This work was supported under NASA Contract No. NNG08FD60C, and made use of data from the NuS-TAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. This research made use of the NuS-TAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center (ASDC, Italy) and the California Institute of Technology (Caltech, USA). MLK gratefully acknowledges support through NASA grant #NNX13AE90G. GM and AM acknowledge financial support from Italian Space Agency under grant ASI/INAF I/037/12/0-011/13 and from the European Union Seventh Framework Programme (FP7/2007(FP7/ -2013 under grant agreement n.312789. CSR acknowledges support from the NASA-ADAP program under grants NNX14AF86G and NNX14AF89G. We thank Alan Marscher for helpful discussions. We thank the anonymous referee for their comments, which have improved this manuscript. This research made use of Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013), and Matplotlib (Hunter 2007).

INFLUENCE OF THE MINOR PILEUP IN SUZAKU XIS DATA
To investigate the influence of the mild pileup in the XIS data on our results, we extract XIS data from an annular region with inner radius 30 ′′ and outer radius 170 ′′ . We fit our best-fit model to this data, and we find consistent parameter values within 90% confidence in one interesting parameter except for instrument cross-normalizations relative to XIS-FI. We find cross-normalizations for XIS-1 of 0.93±0.01, PIN of 1.21±0.01, FPMA of 1.00±0.01, and FPMB of 1.03 ± 0.01. Thus, we find that the minor pileup fraction does not influence our conclusions on the spectral properties of NGC 4151.