Improved Measurement of the Evolution of the Reactor Antineutrino Flux and Spectrum at Daya Bay.

Reactor neutrino experiments play a crucial role in advancing our knowledge of neutrinos. In this Letter, the evolution of the flux and spectrum as a function of the reactor isotopic content is reported in terms of the inverse-beta-decay yield at Daya Bay with 1958 days of data and improved systematic uncertainties. These measurements are compared with two signature model predictions: the Huber-Mueller model based on the conversion method and the SM2018 model based on the summation method. The measured average flux and spectrum, as well as the flux evolution with the ^{239}Pu isotopic fraction, are inconsistent with the predictions of the Huber-Mueller model. In contrast, the SM2018 model is shown to agree with the average flux and its evolution but fails to describe the energy spectrum. Altering the predicted inverse-beta-decay spectrum from ^{239}Pu fission does not improve the agreement with the measurement for either model. The models can be brought into better agreement with the measurements if either the predicted spectrum due to ^{235}U fission is changed or the predicted ^{235}U, ^{238}U, ^{239}Pu, and ^{241}Pu spectra are changed in equal measure.

Precision Measurement of Reactor Antineutrino Oscillation at Kilometer-Scale Baselines by Daya Bay.

We present a new determination of the smallest neutrino mixing angle θ_{13} and the mass-squared difference Δm_{32}^{2} using a final sample of 5.55×10^{6} inverse beta-decay (IBD) candidates with the final-state neutron captured on gadolinium. This sample is selected from the complete dataset obtained by the Daya Bay reactor neutrino experiment in 3158 days of operation. Compared to the previous Daya Bay results, selection of IBD candidates has been optimized, energy calibration refined, and treatment of backgrounds further improved. The resulting oscillation parameters are sin^{2}2θ_{13}=0.0851±0.0024, Δm_{32}^{2}=(2.466±0.060)×10^{-3} eV^{2} for the normal mass ordering or Δm_{32}^{2}=-(2.571±0.060)×10^{-3} eV^{2} for the inverted mass ordering.

First Measurement of High-Energy Reactor Antineutrinos at Daya Bay.

This Letter reports the first measurement of high-energy reactor antineutrinos at Daya Bay, with nearly 9000 inverse beta decay candidates in the prompt energy region of 8-12 MeV observed over 1958 days of data collection. A multivariate analysis is used to separate 2500 signal events from background statistically. The hypothesis of no reactor antineutrinos with neutrino energy above 10 MeV is rejected with a significance of 6.2 standard deviations. A 29% antineutrino flux deficit in the prompt energy region of 8-11 MeV is observed compared to a recent model prediction. We provide the unfolded antineutrino spectrum above 7 MeV as a data-based reference for other experiments. This result provides the first direct observation of the production of antineutrinos from several high-Q_{ß} isotopes in commercial reactors.

Joint Determination of Reactor Antineutrino Spectra from

A joint determination of the reactor antineutrino spectra resulting from the fission of ^{235}U and ^{239}Pu has been carried out by the Daya Bay and PROSPECT Collaborations. This Letter reports the level of consistency of ^{235}U spectrum measurements from the two experiments and presents new results from a joint analysis of both data sets. The measurements are found to be consistent. The combined analysis reduces the degeneracy between the dominant ^{235}U and ^{239}Pu isotopes and improves the uncertainty of the ^{235}U spectral shape to about 3%. The ^{235}U and ^{239}Pu antineutrino energy spectra are unfolded from the jointly deconvolved reactor spectra using the Wiener-SVD unfolding method, providing a data-based reference for other reactor antineutrino experiments and other applications. This is the first measurement of the ^{235}U and ^{239}Pu spectra based on the combination of experiments at low- and highly enriched uranium reactors.

Improved Constraints on Sterile Neutrino Mixing from Disappearance Searches in the MINOS, MINOS+, Daya Bay, and Bugey-3 Experiments.

Searches for electron antineutrino, muon neutrino, and muon antineutrino disappearance driven by sterile neutrino mixing have been carried out by the Daya Bay and MINOS+ collaborations. This Letter presents the combined results of these searches, along with exclusion results from the Bugey-3 reactor experiment, framed in a minimally extended four-neutrino scenario. Significantly improved constraints on the θ_{µe} mixing angle are derived that constitute the most constraining limits to date over five orders of magnitude in the mass-squared splitting Δm_{41}^{2}, excluding the 90% C.L. sterile-neutrino parameter space allowed by the LSND and MiniBooNE observations at 90% CL_{s} for Δm_{41}^{2}<13 eV^{2}. Furthermore, the LSND and MiniBooNE 99% C.L. allowed regions are excluded at 99% CL_{s} for Δm_{41}^{2}<1.6 eV^{2}.

Extraction of the

This Letter reports the first extraction of individual antineutrino spectra from ^{235}U and ^{239}Pu fission and an improved measurement of the prompt energy spectrum of reactor antineutrinos at Daya Bay. The analysis uses 3.5×10^{6} inverse beta-decay candidates in four near antineutrino detectors in 1958 days. The individual antineutrino spectra of the two dominant isotopes, ^{235}U and ^{239}Pu, are extracted using the evolution of the prompt spectrum as a function of the isotope fission fractions. In the energy window of 4-6 MeV, a 7% (9%) excess of events is observed for the ^{235}U (^{239}Pu) spectrum compared with the normalized Huber-Mueller model prediction. The significance of discrepancy is 4.0σ for ^{235}U spectral shape compared with the Huber-Mueller model prediction. The shape of the measured inverse beta-decay prompt energy spectrum disagrees with the prediction of the Huber-Mueller model at 5.3σ. In the energy range of 4-6 MeV, a maximal local discrepancy of 6.3σ is observed.

Measurement of the Electron Antineutrino Oscillation with 1958 Days of Operation at Daya Bay.

We report a measurement of electron antineutrino oscillation from the Daya Bay Reactor Neutrino Experiment with nearly 4 million reactor ν[over ¯]_{e} inverse ß decay candidates observed over 1958 days of data collection. The installation of a flash analog-to-digital converter readout system and a special calibration campaign using different source enclosures reduce uncertainties in the absolute energy calibration to less than 0.5% for visible energies larger than 2 MeV. The uncertainty in the cosmogenic ^{9}Li and ^{8}He background is reduced from 45% to 30% in the near detectors. A detailed investigation of the spent nuclear fuel history improves its uncertainty from 100% to 30%. Analysis of the relative ν[over ¯]_{e} rates and energy spectra among detectors yields sin^{2}2θ_{13}=0.0856±0.0029 and Δm_{32}^{2}=(2.471_{-0.070}^{+0.068})×10^{-3} eV^{2} assuming the normal hierarchy, and Δm_{32}^{2}=-(2.575_{-0.070}^{+0.068})×10^{-3} eV^{2} assuming the inverted hierarchy.

Evolution of the Reactor Antineutrino Flux and Spectrum at Daya Bay.

The Daya Bay experiment has observed correlations between reactor core fuel evolution and changes in the reactor antineutrino flux and energy spectrum. Four antineutrino detectors in two experimental halls were used to identify 2.2 million inverse beta decays (IBDs) over 1230 days spanning multiple fuel cycles for each of six 2.9 GW_{th} reactor cores at the Daya Bay and Ling Ao nuclear power plants. Using detector data spanning effective ^{239}Pu fission fractions F_{239} from 0.25 to 0.35, Daya Bay measures an average IBD yield σ[over ¯]_{f} of (5.90±0.13)×10^{-43} cm^{2}/fission and a fuel-dependent variation in the IBD yield, dσ_{f}/dF_{239}, of (-1.86±0.18)×10^{-43} cm^{2}/fission. This observation rejects the hypothesis of a constant antineutrino flux as a function of the ^{239}Pu fission fraction at 10 standard deviations. The variation in IBD yield is found to be energy dependent, rejecting the hypothesis of a constant antineutrino energy spectrum at 5.1 standard deviations. While measurements of the evolution in the IBD spectrum show general agreement with predictions from recent reactor models, the measured evolution in total IBD yield disagrees with recent predictions at 3.1σ. This discrepancy indicates that an overall deficit in the measured flux with respect to predictions does not result from equal fractional deficits from the primary fission isotopes ^{235}U, ^{239}Pu, ^{238}U, and ^{241}Pu. Based on measured IBD yield variations, yields of (6.17±0.17) and (4.27±0.26)×10^{-43} cm^{2}/fission have been determined for the two dominant fission parent isotopes ^{235}U and ^{239}Pu. A 7.8% discrepancy between the observed and predicted ^{235}U yields suggests that this isotope may be the primary contributor to the reactor antineutrino anomaly.

Improved Search for a Light Sterile Neutrino with the Full Configuration of the Daya Bay Experiment.

This Letter reports an improved search for light sterile neutrino mixing in the electron antineutrino disappearance channel with the full configuration of the Daya Bay Reactor Neutrino Experiment. With an additional 404 days of data collected in eight antineutrino detectors, this search benefits from 3.6 times the statistics available to the previous publication, as well as from improvements in energy calibration and background reduction. A relative comparison of the rate and energy spectrum of reactor antineutrinos in the three experimental halls yields no evidence of sterile neutrino mixing in the 2×10^{-4}≲|Δm_{41}^{2}|≲0.3 eV^{2} mass range. The resulting limits on sin^{2}2θ_{14} are improved by approx imately a factor of 2 over previous results and constitute the most stringent constraints to date in the |Δm_{41}^{2}|≲0.2 eV^{2} region.

Measurement of the Reactor Antineutrino Flux and Spectrum at Daya Bay.

This Letter reports a measurement of the flux and energy spectrum of electron antineutrinos from six 2.9 GWth nuclear reactors with six detectors deployed in two near (effective baselines 512 and 561 m) and one far (1579 m) underground experimental halls in the Daya Bay experiment. Using 217 days of data, 296 721 and 41 589 inverse ß decay (IBD) candidates were detected in the near and far halls, respectively. The measured IBD yield is (1.55±0.04) ×10(-18) cm(2) GW(-1) day(-1) or (5.92±0.14) ×10(-43) cm(2) fission(-1). This flux measurement is consistent with previous short-baseline reactor antineutrino experiments and is 0.946±0.022 (0.991±0.023) relative to the flux predicted with the Huber-Mueller (ILL-Vogel) fissile antineutrino model. The measured IBD positron energy spectrum deviates from both spectral predictions by more than 2σ over the full energy range with a local significance of up to ∼4σ between 4-6 MeV. A reactor antineutrino spectrum of IBD reactions is extracted from the measured positron energy spectrum for model-independent predictions.

New measurement of antineutrino oscillation with the full detector configuration at Daya Bay.

We report a new measurement of electron antineutrino disappearance using the fully constructed Daya Bay Reactor Neutrino Experiment. The final two of eight antineutrino detectors were installed in the summer of 2012. Including the 404 days of data collected from October 2012 to November 2013 resulted in a total exposure of 6.9×10^{5} GW_{th} ton days, a 3.6 times increase over our previous results. Improvements in energy calibration limited variations between detectors to 0.2%. Removal of six ^{241}Am-^{13}C radioactive calibration sources reduced the background by a factor of 2 for the detectors in the experimental hall furthest from the reactors. Direct prediction of the antineutrino signal in the far detectors based on the measurements in the near detectors explicitly minimized the dependence of the measurement on models of reactor antineutrino emission. The uncertainties in our estimates of sin^{2}2θ_{13} and |Δm_{ee}^{2}| were halved as a result of these improvements. An analysis of the relative antineutrino rates and energy spectra between detectors gave sin^{2}2θ_{13}=0.084±0.005 and |Δm_{ee}^{2}|=(2.42±0.11)×10^{-3} eV^{2} in the three-neutrino framework.

Spectral structure of electron antineutrinos from nuclear reactors.

Recent measurements of the positron energy spectrum obtained from inverse beta decay interactions of reactor electron antineutrinos show an excess in the 4 to 6 MeV region relative to current predictions. First-principles calculations of fission and beta decay processes within a typical pressurized water reactor core identify prominent fission daughter isotopes as a possible origin for this excess. These calculations also predict percent-level substructures in the antineutrino spectrum due to Coulomb effects in beta decay. Precise measurement of these substructures can elucidate the nuclear processes occurring within reactors. These substructures can be a systematic issue for measurements utilizing the detailed spectral shape.

Search for a light sterile neutrino at Daya Bay.

A search for light sterile neutrino mixing was performed with the first 217 days of data from the Daya Bay Reactor Antineutrino Experiment. The experiment's unique configuration of multiple baselines from six 2.9 GW(th) nuclear reactors to six antineutrino detectors deployed in two near (effective baselines 512 m and 561 m) and one far (1579 m) underground experimental halls makes it possible to test for oscillations to a fourth (sterile) neutrino in the 10(-3) eV(2)<|Δm(41)(2) |< 0.3 eV(2) range. The relative spectral distortion due to the disappearance of electron antineutrinos was found to be consistent with that of the three-flavor oscillation model. The derived limits on sin(2) 2θ(14) cover the 10(-3) eV(2) ≲ |Δm(41)(2)| ≲ 0.1 eV(2) region, which was largely unexplored.

Spectral measurement of electron antineutrino oscillation amplitude and frequency at Daya Bay.

A measurement of the energy dependence of antineutrino disappearance at the Daya Bay reactor neutrino experiment is reported. Electron antineutrinos (ν¯(e)) from six 2.9 GW(th) reactors were detected with six detectors deployed in two near (effective baselines 512 and 561 m) and one far (1579 m) underground experimental halls. Using 217 days of data, 41 589 (203 809 and 92 912) antineutrino candidates were detected in the far hall (near halls). An improved measurement of the oscillation amplitude sin(2)2θ(13)=0.090(-0.009)(+0.008) and the first direct measurement of the ν¯(e) mass-squared difference |Δm(ee)2|=(2.59(-0.20)(+0.19))×10(-3) eV2 is obtained using the observed ν¯(e) rates and energy spectra in a three-neutrino framework. This value of |Δm(ee)2| is consistent with |Δm(µµ)2| measured by muon neutrino disappearance, supporting the three-flavor oscillation model.

Observation of electron-antineutrino disappearance at Daya Bay.

The Daya Bay Reactor Neutrino Experiment has measured a nonzero value for the neutrino mixing angle θ(13) with a significance of 5.2 standard deviations. Antineutrinos from six 2.9 GWth reactors were detected in six antineutrino detectors deployed in two near (flux-weighted baseline 470 m and 576 m) and one far (1648 m) underground experimental halls. With a 43,000 ton-GWth-day live-time exposure in 55 days, 10,416 (80,376) electron-antineutrino candidates were detected at the far hall (near halls). The ratio of the observed to expected number of antineutrinos at the far hall is R=0.940±0.011(stat.)±0.004(syst.). A rate-only analysis finds sin(2)2θ(13)=0.092±0.016(stat.)±0.005(syst.) in a three-neutrino framework.

Impact of fatty acid composition on the accuracy of mid-infrared fat analysis of farm milks.

Our objective was to determine whether data from a previous study using model milk emulsions to characterize the influence of variation in fatty acid chain length and unsaturation on mid-infrared (MIR) fat predictions could be used to identify a strategy to improve the accuracy of MIR fat predictions on a population of farm milks with a wide variation in fatty acid chain length and unsaturation. The mean fatty acid chain length for 45 farm milks was 14.417 carbons, and the mean unsaturation was 0.337 double bonds per fatty acid. The range of fatty acid chain lengths across the 45 farm milks was 1.23 carbons, and the range in unsaturation was 0.167 double bonds per fatty acid. Fat B (absorbance by the carbon-hydrogen stretch) MIR predictions increased and fat A MIR (absorbance by the ester carbonyl stretch) predictions decreased relative to reference chemistry with increasing fatty acid chain length. When the fat B MIR fat predictions were corrected for sample-to-sample variation in unsaturation, the positive correlation between fat B and fatty acid chain length increased from a coefficient of determination of 0.42 to 0.89. A 45:55 ratio of fat B corrected for unsaturation and fat A gave a smaller standard deviation of the difference between MIR prediction and reference chemistry than any ratio of the fat B (without correction for unsaturation) and fat A or either fat B or fat A alone. This demonstrates the technical feasibility of this approach to improve MIR testing accuracy for fat, if a simple procedure could be developed to determine the unsaturation of fat in milk rapidly and to correct the fat B reading for the effect of unsaturation before being combined with fat A.

Uptake and utilization of trans octadecenoic acids in lactating dairy cows.

Trans fatty acids (FA) arise in ruminant-derived foods as a consequence of rumen biohydrogenation and are of interest because of their biological effects and potential role in chronic human diseases. Our objective was to compare 2 trans FA, elaidic acid (EA; trans-9 18:1) and vaccenic acid (VA; trans-11 18:1), with oleic acid (OA; cis-9 18:1) relative to plasma lipid transport and mammary utilization for milk fat synthesis. Three ruminally cannulated, Holstein dairy cows, 259 +/- 6 DIM (mean +/- SEM), were randomly assigned in a 3 x 3 Latin square design. Treatments were a 4-d abomasal infusion of 1) OA (45.5 g/d), 2) EA (41.7 g/d), and 3) VA (41.4 g/d). Milk samples were collected at each milking and blood samples were collected at the start and end of each treatment period. The proportions of total plasma FA associated with each plasma lipid fraction at baseline (pretreatment) were 62.6 +/- 0.6% phospholipids, 26.1 +/- 0.6% cholesterol esters, 9.8 +/- 0.4% triglycerides, and 1.5 +/- 0.1% nonesterified fatty acids; these values were unaffected by treatment. There were striking differences in the FA composition of the individual plasma lipid fractions and in the distribution of specific 18-carbon FA among the lipid fractions. Infusion of treatment isomers caused their specific increase in the various plasma lipid fractions but had no effect on milk production variables, including milk fat yield and content. Transfer efficiency of infused OA, EA, and VA to milk fat averaged 65.5 +/- 3.0%, 59.7 +/- 1.5%, and 54.3 +/- 0.6%, respectively. For the VA infusion, 24.6 +/- 1.1% of the transfer was accounted for by the increased yield of cis-9, trans-11 conjugated linoleic acid in milk fat, consistent with its endogenous synthesis from VA via the mammary enzyme Delta(9)-desaturase. Notably, linoleic acid (18:2n-6) and linolenic acid (18:3n-3) accounted for 47.7% of total plasma FA, but only 2.6% of FA in milk. Overall, results demonstrate clear differences in plasma transport and mammary uptake and utilization of 18-carbon FA, and these relate to the location, orientation, and number of double bonds.

Effect of supplementation with calcium salts of fish oil on n-3 fatty acids in milk fat.

Enrichment of milk fat with n-3 fatty acids, in particular eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA), may be advantageous because of their beneficial effects on human health. In addition, these fatty acids play an important role in reproductive processes in dairy cows. Our objective was to evaluate the protection of EPA and DHA against rumen biohydrogenation provided by Ca salts of fish oil. Four Holstein cows were assigned in a Latin square design to the following treatments: 1) ruminal infusion of Ca salts of fish oil and palm fatty acid distillate low dose (CaFO-1), 2) ruminal infusion of Ca salts of fish oil and palm fatty acid distillate high dose (CaFO-2), 3) ruminal infusion of fish oil high dose (RFO), and 4) abomasal infusion of fish oil high dose (AFO). The high dose of fish oil provided approximately 16 and approximately 21 g/d of EPA and DHA, respectively, whereas the low dose (CaFO-1) provided 50% of these amounts. A 10-d pretreatment period was used as a baseline, followed by 9-d treatment periods with interceding intervals of 10 d. Supplements were infused every 6 h, milk samples were taken the last 3 d, and plasma samples were collected the last day of baseline and treatment periods. Milk fat content of EPA and DHA were 5 to 6 times greater with AFO, but did not differ among other treatments. Milk and milk protein yield were unaffected by treatment, but milk fat yield and DM intake were reduced by 20 and 15%, respectively, by RFO. Overall, results indicate rumen biohydrogenation of long chain n-3 fatty acids was extensive, averaging >85% for EPA and >75% for DHA for the Ca salts and unprotected fish oil supplements. Thus, Ca salts of fish oil offered no protection against the biohydrogenation of EPA and DHA beyond that observed with unprotected fish oil; however, the Ca salts did provide rumen inertness by preventing the negative effects on DM intake and milk fat yield observed with unprotected fish oil.

Trans-9, cis-11 conjugated linoleic acid reduces milk fat synthesis in lactating dairy cows.

Under certain dietary situations, rumen biohydrogenation results in the production of unique fatty acids that inhibit milk fat synthesis. The first of these to be identified was trans-10, cis-12 conjugated linoleic acid (CLA), but others are postulated to contribute to diet-induced milk fat depression (MFD). Our objective was to examine the potential role of trans-9, cis-11 CLA in the regulation of milk fat. In a preliminary study, we used gas-liquid and high-performance liquid chromatography techniques to examine milk fat samples from a diet-induced MFD study and found that an increase in trans-9, cis-11 CLA corresponded to the decrease in milk fat yield. We investigated this further using a CLA enrichment of 9, 11 isomers to examine the biological effect of trans-9, cis-11 CLA on milk fat synthesis. Four rumen-fistulated Holstein cows were randomly assigned in a 4 x 4 Latin square experiment involving 5-d treatment periods and abomasal infusion of 1) ethanol (control), 2) a 9, 11 CLA mix (containing 32% trans-9, cis-11, 29% cis-9, trans-11, and 17% trans-9, trans-11), 3) a trans-9, trans-11 CLA supplement, and 4) a trans-10, cis-12 CLA supplement (positive control). The trans-9, trans-11 CLA and trans-10, cis-12 CLA supplements were of high purity (>90%), and all supplements were infused at a rate to provide 5 g/d of the CLA isomer of interest. Milk yield and dry matter intake did not differ among treatments. Compared with the control treatment, milk fat yield was reduced by 15% for the 9, 11 CLA mixture and by 27% for the trans-10, cis-12 CLA treatment. We also found that trans-9, trans-11 CLA had no effect on milk fat yield, and previous research has shown that milk fat yield is unaltered when cows are infused with cis-9, trans-11 CLA. When all treatments were considered, results suggested that trans-9, cis-11 was the CLA isomer in the 9, 11 CLA mix responsible for the reduction in milk fat synthesis, although the magnitude was less than that observed for trans-10, cis-12 CLA. Interestingly, trans-9, trans-11 CLA altered the milk fat desaturase index, further demonstrating that alterations in desaturase can occur independently of effects on milk fat synthesis. Overall, our investigations identified that an increase in milk fat content of trans-9, cis-11 CLA was associated with diet-induced MFD and provided evidence of a role for this isomer in MFD based on the 15% reduction in milk fat yield with abomasal infusion of a CLA enrichment that supplied 5 g/d of trans-9, cis-11 CLA.

Response to conjugated linoleic acid in dairy cows differing in energy and protein status.

The trans-10, cis-12 conjugated linoleic acid (CLA) isomer inhibits milk fat synthesis, whereas milk yield and synthesis of other milk components generally remain unchanged in established lactation. However, in some CLA studies increases in milk yield, milk protein yield, or both have been observed in cows limited in energy, either in early lactation or when grazing pasture. Our objective was to evaluate the performance and monitor peripheral tissue responses to homeostatic signals regulating lipolysis and glucose uptake with CLA supplementation when cows were limited in metabolizable energy in combination with moderate or excess metabolizable protein supply. Holstein cows (n = 48; 112 +/- 5 d in milk; mean +/- SE) were provided ad libitum access to a diet that met energy and protein requirements for a 16-d standardization interval. Based on performance during this interval, the Cornell Net Carbohydrate and Protein System was used to design energy-limiting rations that provided 80% of metabolizable energy requirements, and these were fed throughout the treatment periods. Cows were randomly allocated to 4 treatments, in a 2-period crossover design. Treatments were 1) moderate metabolizable protein (MP) supply, 2) moderate MP supply + CLA, 3) excess MP supply, and 4) excess MP supply + CLA. Moderate and excess MP supply were at 88 and 117%, respectively, of the MP requirement established during the standardization period, as estimated by the Cornell Net Carbohydrate and Protein System. Each experimental period comprised 16 d, with crossover of CLA within each protein level. The lipid-encapsulated CLA supplement provided 12 g/d of trans-10, cis-12 CLA. Conjugated linoleic acid treatment reduced milk fat yield by 21% but increased milk yield and milk protein yield by 2.6 and 2.8%, respectively. Milk yield and content and yield of both milk protein and fat were unaltered by either protein treatment alone or in combination with CLA. Basal concentrations of glucose, insulin, and nonesterified fatty acids were unaffected by CLA supplementation. The fractional rate of glucose clearance in response to an insulin challenge and the nonesterified fatty acid response to an epinephrine challenge were also not altered by either CLA treatment or MP supply. Overall, the results demonstrate that CLA supplementation when cows are energy-limited allows for repartitioning of nutrients, resulting in increased yields of milk and milk protein, and this can occur without changes in whole-body glucose homeostasis and adipose tissue response to lipolytic stimuli.