IPhO 2015 · Theory problem 1
Particles from the Sun
Photons from the surface of the Sun and neutrinos from its core can tell us about solar temperatures and also confirm that the Sun shines because of nuclear reactions.
Throughout this problem, take the mass of the Sun to be \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\), its radius, \(R_\odot=7.00\times10^8\ \mathrm{m}\), its luminosity (radiation energy emitted per unit time), \(L_\odot=3.85\times10^{26}\ \mathrm{W}\), and the Earth-Sun distance, \(d_\odot=1.50\times10^{11}\ \mathrm{m}\).
Note:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{constant} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{constant} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{constant} \]
A Radiation from the sun :
| Part | Question | Marks |
|---|---|---|
| A1 | Assume that the Sun radiates like a perfect blackbody. Use this fact to calculate the temperature, \(T_s\), of the solar surface. | 0.3 |
The spectrum of solar radiation can be approximated well by the Wien distribution law. Accordingly, the solar energy incident on any surface on the Earth per unit time per unit frequency interval, \(u(\nu)\), is given by
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
where \(\nu\) is the frequency and \(A\) is the area of the surface normal to the direction of the incident radiation.
Now, consider a solar cell which consists of a thin disc of semiconducting material of area, \(A\), placed perpendicular to the direction of the Sun’s rays.
| Part | Question | Marks |
|---|---|---|
| A2 | Using the Wien approximation, express the total radiated solar power, \(P_{\mathrm{in}}\), incident on the surface of the solar cell, in terms of \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) and the fundamental constants \(c\), \(h\), \(k_B\). | 0.3 |
| A3 | Express the number of photons, \(n_\gamma(\nu)\), per unit time per unit frequency interval incident on the surface of the solar cell in terms of \(A\), \(R_\odot\), \(d_\odot\), \(T_s\), \(\nu\) and the fundamental constants \(c\), \(h\), \(k_B\). | 0.2 |
The semiconducting material of the solar cell has a “band gap” of energy, \(E_g\). We assume the following model. Every photon of energy \(E\geq E_g\) excites an electron across the band gap. This electron contributes an energy, \(E_g\), as the useful output energy, and any extra energy is dissipated as heat (not converted to useful energy).
| Part | Question | Marks |
|---|---|---|
| A4 | Define \(x_g=h\nu_g/k_BT_s\) where \(E_g=h\nu_g\). Express the useful output power of the cell, \(P_{\mathrm{out}}\), in terms of \(x_g\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) and the fundamental constants \(c\), \(h\), \(k_B\). | 1.0 |
| A5 | Express the efficiency, \(\eta\), of this solar cell in terms of \(x_g\). | 0.2 |
| A6 | Make a qualitative sketch of \(\eta\) versus \(x_g\). The values at \(x_g=0\) and \(x_g\to\infty\) should be clearly shown. What is the slope of \(\eta(x_g)\) at \(x_g=0\) and \(x_g\to\infty\)? | 1.0 |
| A7 | Let \(x_0\) be the value of \(x_g\) for which \(\eta\) is maximum. Obtain the cubic equation that gives \(x_0\). Estimate the value of \(x_0\) within an accuracy of \(\pm0.25\). Hence calculate \(\eta(x_0)\). | 1.0 |
| A8 | The band gap of pure silicon is \(E_g=1.11\ \mathrm{eV}\). Calculate the efficiency, \(\eta_{\mathrm{Si}}\), of a silicon solar cell using this value. | 0.2 |
\(^{1}\) Amol Dighe (TIFR), Anwesh Mazumdar (HBCSE-TIFR) and Vijay A. Singh (ex-National Coordinator, Science Olympiads) were the principal authors of this problem. The contributions of the Academic Committee, Academic Development Group and the International Board are gratefully acknowledged.
In the late nineteenth century, Kelvin and Helmholtz (KH) proposed a hypothesis to explain how the Sun shines. They postulated that starting as a very large cloud of matter of mass, \(M_\odot\), and negligible density, the Sun has been shrinking continuously. The shining of the Sun would then be due to the release of gravitational potential energy through this slow contraction.
| Part | Question | Marks |
|---|---|---|
| A9 | Let us assume that the density of matter is uniform inside the Sun. Find the total gravitational potential energy, \(\Omega\), of the Sun at present, in terms of \(G\), \(M_\odot\) and \(R_\odot\). | 0.3 |
| A10 | Estimate the maximum possible time, \(\tau_{\mathrm{KH}}\) (in years), for which the Sun could have been shining, according to the KH hypothesis. Assume that the luminosity of the Sun has been constant throughout this period. | 0.5 |
The \(\tau_{\mathrm{KH}}\) calculated above does not match the age of the solar system estimated from studies of meteorites. This shows that the energy source of the Sun cannot be purely gravitational.
B Neutrinos from the Sun :
In 1938, Hans Bethe proposed that nuclear fusion of hydrogen into helium in the core of the Sun is the source of its energy. The net nuclear reaction is:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
The “electron neutrinos”, \(\nu_e\), produced in this reaction may be taken to be massless. They escape the Sun and their detection on the Earth confirms the occurrence of nuclear reactions inside the Sun. Energy carried away by the neutrinos can be neglected in this problem.
| Part | Question | Marks |
|---|---|---|
| B1 | Calculate the flux density, \(\Phi_\nu\), of the number of neutrinos arriving at the Earth, in units of \(\mathrm{m^{-2}\,s^{-1}}\). The energy released in the above reaction is \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\). Assume that the energy radiated by the Sun is entirely due to this reaction. | 0.6 |
Travelling from the core of the Sun to the Earth, some of the electron neutrinos, \(\nu_e\), are converted to other types of neutrinos, \(\nu_x\). The efficiency of the detector for detecting \(\nu_x\) is 1/6 of its efficiency for detecting \(\nu_e\). If there is no neutrino conversion, we expect to detect an average of \(N_1\) neutrinos in a year. However, due to the conversion, an average of \(N_2\) neutrinos (\(\nu_e\) and \(\nu_x\) combined) are actually detected per year.
| Part | Question | Marks |
|---|---|---|
| B2 | In terms of \(N_1\) and \(N_2\), calculate what fraction, \(f\), of \(\nu_e\) is converted to \(\nu_x\). | 0.4 |
In order to detect neutrinos, large detectors filled with water are constructed. Although the interactions of neutrinos with matter are very rare, occasionally they knock out electrons from water molecules in the detector. These energetic electrons move through water at high speeds, emitting electromagnetic radiation in the process. As long as the speed of such an electron is greater than the speed of light in water (refractive index, \(n\)), this radiation, called Cherenkov radiation, is emitted in the shape of a cone.
| Part | Question | Marks |
|---|---|---|
| B3 | Assume that an electron knocked out by a neutrino loses energy at a constant rate of \(\alpha\) per unit time, while it travels through water. If this electron emits Cherenkov radiation for a time, \(\Delta t\), determine the energy imparted to this electron (\(E_{\mathrm{imparted}}\)) by the neutrino, in terms of \(\alpha\), \(\Delta t\), \(n\), \(m_e\) and \(c\). (Assume the electron to be at rest before its interaction with the neutrino.) | 2.0 |
The fusion of H into He inside the Sun takes place in several steps. Nucleus of \({}^{7}\mathrm{Be}\) (rest mass, \(m_{\mathrm{Be}}\)) is produced in one of these intermediate steps. Subsequently, it can absorb an electron, producing a \({}^{7}\mathrm{Li}\) nucleus (rest mass, \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)) and emitting a \(\nu_e\). The corresponding nuclear reaction is:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
When a Be nucleus (\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\)) is at rest and absorbs an electron also at rest, the emitted neutrino has energy \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\). However, the Be nuclei are in random thermal motion due to the temperature \(T_c\) at the core of the Sun, and act as moving neutrino sources. As a result, the energy of emitted neutrinos fluctuates with a root mean square (rms) value \(\Delta E_{\mathrm{rms}}\).
| Part | Question | Marks |
|---|---|---|
| B4 | If \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\), calculate the rms speed of the Be nuclei, \(V_{\mathrm{Be}}\), and hence estimate \(T_c\). (Hint: \(\Delta E_{\mathrm{rms}}\) depends on the rms value of the component of velocity along the line of sight). | 2.0 |
Official English text of the problem.
Answers
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 The curve starts at the origin with slope \(1/3\), reaches one maximum near \((2.27,0.457)\), and approaches the horizontal axis from above: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{from below}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
Solution
For numerical evaluation, use the standard constants \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] with \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) and \(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\).
A1
By the Stefan–Boltzmann law, a blackbody emits power \(\sigma T^4\) per unit surface area. The Sun’s emitting area is \(4\pi R_\odot^2\), so \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] Therefore, \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
The total incident power is the integral over all frequencies: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] Set \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] Then \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] The supplied integration formula gives \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] Thus \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] All efficiencies below use this same Wien-approximation input power.
A3
Each photon of frequency \(\nu\) carries energy \(h\nu\). Dividing the spectral power by this energy gives the spectral photon arrival rate: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] Hence \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] Its units are photons per second per hertz.
A4
Only photons with \(\nu\geq\nu_g\) contribute, and each contributes useful energy \(E_g\). Consequently, \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] Using \(E_g=x_gk_BT_s\) and the substitution from A2, \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] Now \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] Therefore, \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
Efficiency is useful output power divided by incident power: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] The common dimensional and geometrical factors cancel, giving \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \]
A6
At zero band gap, \[ \boxed{\eta(0)=0}. \] Although photons can then excite electrons, the useful energy per electron is zero.
For large \(x_g\), \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] so \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] The exponentially small number of photons above the gap makes the useful output vanish.
Differentiating, \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] Thus \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{from below}}. \]
The sketch starts at the origin with tangent \(\eta=x_g/3\), rises to a single maximum near \(x_g=2.27\), and then decreases toward the horizontal asymptote \(\eta=0\), remaining positive. For additional shape information, \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] so the curve is concave downward for \(0<x_g<4\) and concave upward for \(x_g>4\).
Representative points are:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
At the maximum, \(\eta'(x_0)=0\). Since the exponential factor never vanishes, \[ -x_0^3+x_0^2+2x_0+2=0, \] or \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
For \(x_0>0\), this is equivalent to \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] The left-hand side increases strictly with \(x_0\), while the right-hand side decreases strictly. Therefore, there is only one positive stationary point. Together with the endpoint behaviour in A6, this establishes that it is the global maximum.
Let \(F(x)=x^3-x^2-2x-2\). Then \[ F(2)=-2,\qquad F(2.5)=2.375. \] Thus \(2<x_0<2.5\), already giving \(x_0=2.25\pm0.25\). Refinement gives \[ \boxed{x_0\simeq2.27}. \] Substituting into the efficiency, \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] Hence \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
At the surface temperature found in A1, \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] For silicon, \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] Therefore, \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] Thus \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] This is very close to the maximum because the silicon gap is close to the optimal gap \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\).
A9
Build the uniform sphere by bringing successive thin spherical shells from infinity. At radius \(r\), the mass already assembled is \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] and the shell mass is \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] By the shell theorem, the potential energy added is \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] Therefore, \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] Hence \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] The negative sign corresponds to taking zero potential energy at infinite separation.
A10
The initial very extended cloud has approximately zero gravitational potential energy. Its contraction therefore releases at most \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] The maximum possible shining time follows by taking all this released energy to have been radiated. At constant luminosity, \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] so \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] Numerically, \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] Any released energy retained as internal energy would reduce the duration supported by gravitational contraction, so this is an upper bound.
B1
Neglecting energy carried away by neutrinos, the number of fusion reactions per second is \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] Each reaction produces two neutrinos: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] For spherical emission, these neutrinos cross an area \(4\pi d_\odot^2\) at Earth’s distance. Thus \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] Substitution gives \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
After conversion, the fraction \(1-f\) remains as electron neutrinos and contributes \((1-f)N_1\) detections per year. The converted fraction contributes \(fN_1/6\), because its detection efficiency is six times smaller.
Consequently, \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] Solving, \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] The physical range \(0\leq f\leq1\) corresponds to \(N_1/6\leq N_2\leq N_1\).
B3
Cherenkov radiation ceases when the electron slows to the speed of light in water: \[ v_{\mathrm{th}}=\frac{c}{n}. \] The electron’s kinetic energy at that instant is relativistic: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
During the emission time \(\Delta t\), the electron loses energy \(\alpha\Delta t\). Its initial kinetic energy immediately after the neutrino interaction is therefore \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] Since the electron was initially at rest before the interaction, the energy imparted by the neutrino equals \(K_i\). Hence \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] The second term is essential: the electron still has kinetic energy when Cherenkov emission stops.
B4
Use the stated model in which the energy spread comes from thermal motion of the Be sources. Let \(v_\parallel\) be the source velocity component toward Earth.
For a massless neutrino, momentum in the source rest frame is \(p_\nu=E_\nu/c\). To first order in the nonrelativistic source velocity, the Lorentz transformation of energy gives \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] Thus the Doppler energy shift is \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
The thermal motion is isotropic, so \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] Taking the rms energy shift, \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] Therefore, \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] Numerically, \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] and hence \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] The ratio \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) validates the nonrelativistic approximation.
By equipartition for the three translational degrees of freedom, \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] Thus \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] Substituting the data, \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] which gives \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
The solution was checked against the official solution.
IPhO 2015 · Theorieaufgabe 1
Teilchen von der Sonne
Photonen von der Sonnenoberfläche und Neutrinos aus ihrem Kern können uns Aufschluss über die Temperaturen der Sonne geben und auch bestätigen, dass die Sonne aufgrund von Kernreaktionen leuchtet.
Nehmen Sie in dieser gesamten Aufgabe für die Masse der Sonne \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\), für ihren Radius \(R_\odot=7.00\times10^8\ \mathrm{m}\), für ihre Leuchtkraft (pro Zeiteinheit abgestrahlte Energie) \(L_\odot=3.85\times10^{26}\ \mathrm{W}\) und für den Abstand zwischen Erde und Sonne \(d_\odot=1.50\times10^{11}\ \mathrm{m}\) an.
Hinweis:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{Konstante} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{Konstante} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{Konstante} \]
A Strahlung von der Sonne :
| Teil | Frage | Punkte |
|---|---|---|
| A1 | Nehmen Sie an, dass die Sonne wie ein idealer schwarzer Körper strahlt. Berechnen Sie damit die Temperatur \(T_s\) der Sonnenoberfläche. | 0.3 |
Das Spektrum der Sonnenstrahlung lässt sich gut durch das Wiensche Strahlungsgesetz annähern. Demnach ist die pro Zeiteinheit und pro Frequenzintervall auf eine beliebige Fläche auf der Erde einfallende Sonnenenergie \(u(\nu)\) gegeben durch
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
wobei \(\nu\) die Frequenz und \(A\) der Flächeninhalt der zur Richtung der einfallenden Strahlung senkrechten Fläche ist.
Betrachten Sie nun eine Solarzelle, die aus einer dünnen Scheibe aus Halbleitermaterial mit dem Flächeninhalt \(A\) besteht und senkrecht zur Richtung der Sonnenstrahlen angeordnet ist.
| Teil | Frage | Punkte |
|---|---|---|
| A2 | Drücken Sie unter Verwendung der Wien-Näherung die gesamte auf die Oberfläche der Solarzelle einfallende Sonnenstrahlungsleistung \(P_{\mathrm{in}}\) durch \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) und die Naturkonstanten \(c\), \(h\), \(k_B\) aus. | 0.3 |
| A3 | Drücken Sie die Anzahl der pro Zeiteinheit und pro Frequenzintervall auf die Oberfläche der Solarzelle einfallenden Photonen \(n_\gamma(\nu)\) durch \(A\), \(R_\odot\), \(d_\odot\), \(T_s\), \(\nu\) und die Naturkonstanten \(c\), \(h\), \(k_B\) aus. | 0.2 |
Das Halbleitermaterial der Solarzelle hat eine „Bandlücke“ mit der Energie \(E_g\). Wir nehmen folgendes Modell an. Jedes Photon mit der Energie \(E\geq E_g\) regt ein Elektron über die Bandlücke hinweg an. Dieses Elektron liefert die Energie \(E_g\) als nutzbare Ausgangsenergie; jede darüber hinausgehende Energie wird als Wärme abgegeben (und nicht in nutzbare Energie umgewandelt).
| Teil | Frage | Punkte |
|---|---|---|
| A4 | Definieren Sie \(x_g=h\nu_g/k_BT_s\), wobei \(E_g=h\nu_g\). Drücken Sie die nutzbare Ausgangsleistung der Zelle \(P_{\mathrm{out}}\) durch \(x_g\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) und die Naturkonstanten \(c\), \(h\), \(k_B\) aus. | 1.0 |
| A5 | Drücken Sie den Wirkungsgrad \(\eta\) dieser Solarzelle durch \(x_g\) aus. | 0.2 |
| A6 | Skizzieren Sie qualitativ \(\eta\) in Abhängigkeit von \(x_g\). Die Werte bei \(x_g=0\) und für \(x_g\to\infty\) sollen deutlich erkennbar sein. Wie groß ist die Steigung von \(\eta(x_g)\) bei \(x_g=0\) und für \(x_g\to\infty\)? | 1.0 |
| A7 | Es sei \(x_0\) der Wert von \(x_g\), für den \(\eta\) maximal ist. Leiten Sie die kubische Gleichung her, die \(x_0\) bestimmt. Schätzen Sie den Wert von \(x_0\) mit einer Genauigkeit von \(\pm0.25\) ab. Berechnen Sie daraus \(\eta(x_0)\). | 1.0 |
| A8 | Die Bandlücke von reinem Silizium beträgt \(E_g=1.11\ \mathrm{eV}\). Berechnen Sie mit diesem Wert den Wirkungsgrad \(\eta_{\mathrm{Si}}\) einer Siliziumsolarzelle. | 0.2 |
\(^{1}\) Amol Dighe (TIFR), Anwesh Mazumdar (HBCSE-TIFR) und Vijay A. Singh (ehemaliger nationaler Koordinator der Wissenschaftsolympiaden) waren die Hauptautoren dieser Aufgabe. Für die Beiträge des Akademischen Komitees, der Akademischen Entwicklungsgruppe und des Internationalen Gremiums wird herzlich gedankt.
Im späten neunzehnten Jahrhundert stellten Kelvin und Helmholtz (KH) eine Hypothese auf, um zu erklären, wie die Sonne leuchtet. Sie nahmen an, dass die Sonne ursprünglich eine sehr große Materiewolke mit der Masse \(M_\odot\) und vernachlässigbarer Dichte war und seither kontinuierlich schrumpft. Das Leuchten der Sonne wäre dann auf die Freisetzung potenzieller Gravitationsenergie durch diese langsame Kontraktion zurückzuführen.
| Teil | Frage | Punkte |
|---|---|---|
| A9 | Nehmen wir an, dass die Materiedichte im Inneren der Sonne gleichmäßig ist. Bestimmen Sie die gesamte gegenwärtige potenzielle Gravitationsenergie \(\Omega\) der Sonne in Abhängigkeit von \(G\), \(M_\odot\) und \(R_\odot\). | 0.3 |
| A10 | Schätzen Sie die maximal mögliche Zeit \(\tau_{\mathrm{KH}}\) (in Jahren) ab, während der die Sonne gemäß der KH-Hypothese geleuchtet haben könnte. Nehmen Sie an, dass die Leuchtkraft der Sonne während dieses gesamten Zeitraums konstant war. | 0.5 |
Die oben berechnete Zeit \(\tau_{\mathrm{KH}}\) stimmt nicht mit dem aus Meteoritenuntersuchungen abgeschätzten Alter des Sonnensystems überein. Dies zeigt, dass die Energiequelle der Sonne nicht rein gravitativer Natur sein kann.
B Neutrinos von der Sonne :
Im Jahr 1938 schlug Hans Bethe vor, dass die Kernfusion von Wasserstoff zu Helium im Kern der Sonne ihre Energiequelle ist. Die Nettokernreaktion lautet:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
Die bei dieser Reaktion erzeugten „Elektron-Neutrinos“ \(\nu_e\) können als masselos betrachtet werden. Sie entweichen aus der Sonne, und ihr Nachweis auf der Erde bestätigt, dass im Inneren der Sonne Kernreaktionen stattfinden. Die von den Neutrinos abgeführte Energie kann in dieser Aufgabe vernachlässigt werden.
| Teil | Frage | Punkte |
|---|---|---|
| B1 | Berechnen Sie die Teilchenflussdichte \(\Phi_\nu\) der auf der Erde eintreffenden Neutrinos in Einheiten von \(\mathrm{m^{-2}\,s^{-1}}\). Die bei der obigen Reaktion freigesetzte Energie beträgt \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\). Nehmen Sie an, dass die von der Sonne abgestrahlte Energie vollständig aus dieser Reaktion stammt. | 0.6 |
Auf dem Weg vom Sonnenkern zur Erde wandelt sich ein Teil der Elektron-Neutrinos \(\nu_e\) in andere Neutrinoarten \(\nu_x\) um. Die Nachweiseffizienz des Detektors für \(\nu_x\) beträgt 1/6 seiner Nachweiseffizienz für \(\nu_e\). Ohne Neutrinoumwandlung erwarten wir, im Mittel \(N_1\) Neutrinos pro Jahr nachzuweisen. Aufgrund der Umwandlung werden jedoch tatsächlich im Mittel \(N_2\) Neutrinos (\(\nu_e\) und \(\nu_x\) zusammen) pro Jahr nachgewiesen.
| Teil | Frage | Punkte |
|---|---|---|
| B2 | Berechnen Sie in Abhängigkeit von \(N_1\) und \(N_2\), welcher Anteil \(f\) der \(\nu_e\) in \(\nu_x\) umgewandelt wird. | 0.4 |
Zum Nachweis von Neutrinos werden große, mit Wasser gefüllte Detektoren gebaut. Obwohl Wechselwirkungen von Neutrinos mit Materie sehr selten sind, schlagen sie gelegentlich Elektronen aus Wassermolekülen im Detektor heraus. Diese energiereichen Elektronen bewegen sich mit hohen Geschwindigkeiten durch das Wasser und senden dabei elektromagnetische Strahlung aus. Solange die Geschwindigkeit eines solchen Elektrons größer als die Lichtgeschwindigkeit im Wasser (Brechungsindex \(n\)) ist, wird diese als Tscherenkow-Strahlung bezeichnete Strahlung kegelförmig ausgesandt.
| Teil | Frage | Punkte |
|---|---|---|
| B3 | Nehmen Sie an, dass ein von einem Neutrino herausgeschlagenes Elektron beim Durchgang durch Wasser mit einer konstanten Rate von \(\alpha\) pro Zeiteinheit Energie verliert. Wenn dieses Elektron während einer Zeit \(\Delta t\) Tscherenkow-Strahlung aussendet, bestimmen Sie die vom Neutrino auf dieses Elektron übertragene Energie (\(E_{\mathrm{imparted}}\)) in Abhängigkeit von \(\alpha\), \(\Delta t\), \(n\), \(m_e\) und \(c\). (Nehmen Sie an, dass das Elektron vor seiner Wechselwirkung mit dem Neutrino ruht.) | 2.0 |
Die Fusion von H zu He im Inneren der Sonne erfolgt in mehreren Schritten. In einem dieser Zwischenschritte entsteht ein \({}^{7}\mathrm{Be}\)-Kern (Ruhemasse \(m_{\mathrm{Be}}\)). Anschließend kann er ein Elektron einfangen, wobei ein \({}^{7}\mathrm{Li}\)-Kern (Ruhemasse \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)) entsteht und ein \(\nu_e\) ausgesandt wird. Die entsprechende Kernreaktion lautet:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
Wenn ein Be-Kern (\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\)) ruht und ein ebenfalls ruhendes Elektron einfängt, hat das ausgesandte Neutrino die Energie \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\). Aufgrund der Temperatur \(T_c\) im Sonnenkern führen die Be-Kerne jedoch eine ungeordnete thermische Bewegung aus und wirken als bewegte Neutrinoquellen. Daher schwankt die Energie der ausgesandten Neutrinos mit einem Effektivwert (rms) der Abweichung von \(\Delta E_{\mathrm{rms}}\).
| Teil | Frage | Punkte |
|---|---|---|
| B4 | Berechnen Sie für \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\) die mittlere quadratische Geschwindigkeit \(V_{\mathrm{Be}}\) der Be-Kerne und schätzen Sie daraus \(T_c\) ab. (Hinweis: \(\Delta E_{\mathrm{rms}}\) hängt vom Effektivwert der Geschwindigkeitskomponente entlang der Sichtlinie ab.) | 2.0 |
Aufgabentext aus dem offiziellen englischen Text übersetzt.
Antworten
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 Die Kurve beginnt im Ursprung mit der Steigung \(1/3\), erreicht ein einziges Maximum nahe \((2.27,0.457)\) und nähert sich der horizontalen Achse von oben: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{von unten}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
Lösung
Für die numerische Auswertung verwenden wir die Standardkonstanten \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] mit \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) und \(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\).
A1
Nach dem Stefan–Boltzmann-Gesetz strahlt ein schwarzer Körper die Leistung \(\sigma T^4\) pro Flächeneinheit ab. Die abstrahlende Fläche der Sonne beträgt \(4\pi R_\odot^2\), also gilt \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] Daraus folgt \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
Die gesamte einfallende Leistung ergibt sich durch Integration über alle Frequenzen: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] Wir setzen \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] Dann gilt \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] Die angegebene Integrationsformel liefert \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] Somit ist \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] Für alle folgenden Wirkungsgrade wird dieselbe einfallende Leistung in der Wien-Näherung verwendet.
A3
Jedes Photon der Frequenz \(\nu\) trägt die Energie \(h\nu\). Division der spektralen Leistung durch diese Energie ergibt die spektrale Rate der eintreffenden Photonen: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] Damit folgt \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] Die Einheit ist Photonen pro Sekunde und pro Hertz.
A4
Nur Photonen mit \(\nu\geq\nu_g\) tragen bei, und jedes liefert die nutzbare Energie \(E_g\). Folglich gilt \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] Mit \(E_g=x_gk_BT_s\) und der Substitution aus A2 ergibt sich \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] Nun ist \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] Daher folgt \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
Der Wirkungsgrad ist die nutzbare Ausgangsleistung geteilt durch die einfallende Leistung: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] Die gemeinsamen dimensionsbehafteten und geometrischen Faktoren kürzen sich heraus, sodass sich ergibt \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \]
A6
Bei verschwindender Bandlücke gilt \[ \boxed{\eta(0)=0}. \] Obwohl Photonen dann Elektronen anregen können, ist die nutzbare Energie pro Elektron null.
Für große \(x_g\) gilt \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] also \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] Die exponentiell kleine Anzahl der Photonen mit Energien oberhalb der Bandlücke lässt die nutzbare Ausgangsleistung verschwinden.
Durch Ableiten erhalten wir \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] Damit gilt \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{von unten}}. \]
Die Skizze beginnt im Ursprung mit der Tangente \(\eta=x_g/3\), steigt bis zu einem einzigen Maximum nahe \(x_g=2.27\) an und fällt anschließend gegen die horizontale Asymptote \(\eta=0\) ab, wobei sie positiv bleibt. Weitere Informationen über die Kurvenform liefert \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] wonach die Kurve für \(0<x_g<4\) nach unten und für \(x_g>4\) nach oben gekrümmt ist.
Repräsentative Punkte sind:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
Am Maximum gilt \(\eta'(x_0)=0\). Da der Exponentialfaktor nie verschwindet, gilt \[ -x_0^3+x_0^2+2x_0+2=0, \] oder \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
Für \(x_0>0\) ist dies äquivalent zu \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] Die linke Seite wächst streng monoton mit \(x_0\), während die rechte Seite streng monoton abnimmt. Daher gibt es nur eine positive stationäre Stelle. Zusammen mit dem Randverhalten aus A6 zeigt dies, dass dort das globale Maximum liegt.
Sei \(F(x)=x^3-x^2-2x-2\). Dann gilt \[ F(2)=-2,\qquad F(2.5)=2.375. \] Damit ist \(2<x_0<2.5\), was bereits \(x_0=2.25\pm0.25\) ergibt. Eine genauere Bestimmung liefert \[ \boxed{x_0\simeq2.27}. \] Einsetzen in den Wirkungsgrad ergibt \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] Somit folgt \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
Bei der in A1 bestimmten Oberflächentemperatur gilt \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] Für Silizium ist \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] Daher ergibt sich \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] Somit ist \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] Dies liegt sehr nahe am Maximum, weil die Bandlücke von Silizium nahe an der optimalen Bandlücke \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\) liegt.
A9
Wir bauen die homogene Kugel auf, indem wir nacheinander dünne Kugelschalen aus dem Unendlichen heranbringen. Beim Radius \(r\) beträgt die bereits zusammengefügte Masse \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] und die Masse der Schale ist \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] Nach dem Schalentheorem beträgt die hinzugefügte potenzielle Energie \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] Daher gilt \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] Damit folgt \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] Das negative Vorzeichen entspricht der Festlegung, dass die potenzielle Energie bei unendlichem Abstand null ist.
A10
Die anfängliche, sehr ausgedehnte Wolke hat näherungsweise die potenzielle Gravitationsenergie null. Ihre Kontraktion setzt daher höchstens folgende Energie frei: \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] Die maximal mögliche Leuchtdauer ergibt sich unter der Annahme, dass diese gesamte freigesetzte Energie abgestrahlt wurde. Bei konstanter Leuchtkraft gilt \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] also \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] Numerisch ergibt sich \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] Jeder Anteil der freigesetzten Energie, der als innere Energie erhalten bleibt, würde die durch Gravitationskontraktion ermöglichte Leuchtdauer verkürzen; dies ist also eine obere Grenze.
B1
Wenn die von Neutrinos abgeführte Energie vernachlässigt wird, beträgt die Anzahl der Fusionsreaktionen pro Sekunde \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] Jede Reaktion erzeugt zwei Neutrinos: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] Bei kugelsymmetrischer Emission durchqueren diese Neutrinos im Abstand der Erde eine Fläche von \(4\pi d_\odot^2\). Daher gilt \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] Einsetzen ergibt \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
Nach der Umwandlung verbleibt der Anteil \(1-f\) als Elektron-Neutrinos und trägt \((1-f)N_1\) Nachweise pro Jahr bei. Der umgewandelte Anteil trägt \(fN_1/6\) bei, da seine Nachweiseffizienz um den Faktor sechs kleiner ist.
Folglich gilt \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] Auflösen ergibt \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] Der physikalische Bereich \(0\leq f\leq1\) entspricht \(N_1/6\leq N_2\leq N_1\).
B3
Die Tscherenkow-Strahlung endet, wenn das Elektron auf die Lichtgeschwindigkeit im Wasser abgebremst ist: \[ v_{\mathrm{th}}=\frac{c}{n}. \] Die kinetische Energie des Elektrons zu diesem Zeitpunkt ist relativistisch: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
Während der Emissionsdauer \(\Delta t\) verliert das Elektron die Energie \(\alpha\Delta t\). Seine anfängliche kinetische Energie unmittelbar nach der Wechselwirkung mit dem Neutrino beträgt daher \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] Da das Elektron vor der Wechselwirkung zunächst ruhte, ist die vom Neutrino übertragene Energie gleich \(K_i\). Damit folgt \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] Der zweite Term ist wesentlich: Das Elektron besitzt noch kinetische Energie, wenn die Tscherenkow-Emission endet.
B4
Wir verwenden das angegebene Modell, in dem die Energiestreuung durch die thermische Bewegung der Be-Quellen entsteht. Sei \(v_\parallel\) die Geschwindigkeitskomponente der Quelle in Richtung Erde.
Für ein masseloses Neutrino beträgt der Impuls im Ruhesystem der Quelle \(p_\nu=E_\nu/c\). In erster Ordnung in der nichtrelativistischen Quellengeschwindigkeit liefert die Lorentz-Transformation der Energie \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] Somit beträgt die Doppler-Energieverschiebung \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
Die thermische Bewegung ist isotrop, also gilt \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] Für den Effektivwert der Energieverschiebung ergibt sich \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] Daher gilt \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] Numerisch ist \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] und damit \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] Das Verhältnis \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) bestätigt die Gültigkeit der nichtrelativistischen Näherung.
Nach dem Gleichverteilungssatz für die drei Translationsfreiheitsgrade gilt \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] Somit folgt \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] Einsetzen der Daten ergibt \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] woraus folgt \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
Die Lösung wurde mit der offiziellen Lösung abgeglichen.
IPhO 2015 · Teori sorusu 1
Güneş'ten Gelen Parçacıklar
Güneş'in yüzeyinden gelen fotonlar ve çekirdeğinden gelen nötrinolar, Güneş'in sıcaklıkları hakkında bilgi verebilir ve Güneş'in nükleer tepkimeler sayesinde parladığını da doğrulayabilir.
Bu soru boyunca Güneş'in kütlesini \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\), yarıçapını \(R_\odot=7.00\times10^8\ \mathrm{m}\), ışınım gücünü (birim zamanda yayılan ışınım enerjisini) \(L_\odot=3.85\times10^{26}\ \mathrm{W}\) ve Dünya–Güneş uzaklığını \(d_\odot=1.50\times10^{11}\ \mathrm{m}\) alınız.
Not:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{sabit} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{sabit} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{sabit} \]
A Güneş'ten gelen ışınım :
| Bölüm | Soru | Puan |
|---|---|---|
| A1 | Güneş'in ideal bir kara cisim gibi ışınım yaptığını varsayınız. Bu bilgiyi kullanarak Güneş yüzeyinin sıcaklığını, \(T_s\), hesaplayınız. | 0.3 |
Güneş ışınımının spektrumu, Wien dağılım yasasıyla iyi bir yaklaşıklıkla ifade edilebilir. Buna göre, Dünya üzerindeki herhangi bir yüzeye birim zamanda ve birim frekans aralığında gelen Güneş enerjisi, \(u(\nu)\),
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
ile verilir; burada \(\nu\) frekans, \(A\) ise gelen ışınımın doğrultusuna dik yüzeyin alanıdır.
Şimdi, Güneş ışınlarının doğrultusuna dik yerleştirilmiş, alanı \(A\) olan ince bir yarı iletken malzeme diskinden oluşan bir güneş pili ele alınız.
| Bölüm | Soru | Puan |
|---|---|---|
| A2 | Wien yaklaşımını kullanarak güneş pilinin yüzeyine gelen toplam Güneş ışınım gücünü, \(P_{\mathrm{in}}\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) ve temel sabitler \(c\), \(h\), \(k_B\) cinsinden ifade ediniz. | 0.3 |
| A3 | Güneş pilinin yüzeyine birim zamanda ve birim frekans aralığında gelen foton sayısını, \(n_\gamma(\nu)\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\), \(\nu\) ve temel sabitler \(c\), \(h\), \(k_B\) cinsinden ifade ediniz. | 0.2 |
Güneş pilinin yarı iletken malzemesinin “bant aralığı” enerjisi \(E_g\)'dir. Aşağıdaki modeli varsayıyoruz. Enerjisi \(E\geq E_g\) olan her foton, bir elektronu bant aralığını aşacak şekilde uyarır. Bu elektron, yararlı çıkış enerjisine \(E_g\) kadar enerji katkısında bulunur ve fazla enerji ısı olarak dağılır (yararlı enerjiye dönüşmez).
| Bölüm | Soru | Puan |
|---|---|---|
| A4 | \(E_g=h\nu_g\) olmak üzere \(x_g=h\nu_g/k_BT_s\) tanımını yapınız. Pilin yararlı çıkış gücünü, \(P_{\mathrm{out}}\), \(x_g\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) ve temel sabitler \(c\), \(h\), \(k_B\) cinsinden ifade ediniz. | 1.0 |
| A5 | Bu güneş pilinin verimini, \(\eta\), \(x_g\) cinsinden ifade ediniz. | 0.2 |
| A6 | \(\eta\)'nın \(x_g\)'ye bağlı değişimini gösteren nitel bir grafik çiziniz. \(x_g=0\) ve \(x_g\to\infty\) için değerler açıkça gösterilmelidir. \(x_g=0\) ve \(x_g\to\infty\) için \(\eta(x_g)\)'nin eğimi nedir? | 1.0 |
| A7 | \(\eta\)'nın maksimum olduğu \(x_g\) değeri \(x_0\) olsun. \(x_0\)'ı veren üçüncü dereceden denklemi elde ediniz. \(x_0\)'ın değerini \(\pm0.25\) doğrulukla tahmin ediniz. Buradan \(\eta(x_0)\)'ı hesaplayınız. | 1.0 |
| A8 | Saf silisyumun bant aralığı \(E_g=1.11\ \mathrm{eV}\)'tur. Bu değeri kullanarak silisyumdan yapılmış bir güneş pilinin verimini, \(\eta_{\mathrm{Si}}\), hesaplayınız. | 0.2 |
\(^{1}\) Bu sorunun başlıca yazarları Amol Dighe (TIFR), Anwesh Mazumdar (HBCSE-TIFR) ve Vijay A. Singh'dir (Bilim Olimpiyatları eski Ulusal Koordinatörü). Akademik Komite, Akademik Geliştirme Grubu ve Uluslararası Kurulun katkıları için teşekkür ederiz.
On dokuzuncu yüzyılın sonlarında Kelvin ve Helmholtz (KH), Güneş'in nasıl parladığını açıklamak için bir hipotez öne sürdüler. Güneş'in, kütlesi \(M_\odot\) olan ve yoğunluğu ihmal edilebilen çok büyük bir madde bulutu olarak başlayıp sürekli büzüldüğünü varsaydılar. Bu durumda Güneş'in parlaması, bu yavaş büzülme sırasında kütleçekimsel potansiyel enerjinin açığa çıkmasından kaynaklanacaktı.
| Bölüm | Soru | Puan |
|---|---|---|
| A9 | Güneş'in içindeki madde yoğunluğunun düzgün olduğunu varsayalım. Güneş'in bugünkü toplam kütleçekimsel potansiyel enerjisini, \(\Omega\), \(G\), \(M_\odot\) ve \(R_\odot\) cinsinden bulunuz. | 0.3 |
| A10 | KH hipotezine göre Güneş'in parlıyor olabileceği mümkün olan en uzun süreyi, \(\tau_{\mathrm{KH}}\), (yıl cinsinden) tahmin ediniz. Güneş'in ışınım gücünün bu süre boyunca sabit kaldığını varsayınız. | 0.5 |
Yukarıda hesaplanan \(\tau_{\mathrm{KH}}\), göktaşları üzerindeki çalışmalardan tahmin edilen Güneş Sistemi'nin yaşıyla uyuşmaz. Bu, Güneş'in enerji kaynağının yalnızca kütleçekimi olamayacağını gösterir.
B Güneş'ten gelen nötrinolar :
1938'de Hans Bethe, Güneş'in enerji kaynağının, çekirdeğinde hidrojenin nükleer füzyonla helyuma dönüşmesi olduğunu öne sürdü. Net nükleer tepkime şöyledir:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
Bu tepkimede üretilen “elektron nötrinoları”, \(\nu_e\), kütlesiz kabul edilebilir. Güneş'ten kaçarlar ve Dünya'da algılanmaları, Güneş'in içinde nükleer tepkimeler gerçekleştiğini doğrular. Bu soruda nötrinoların taşıyıp götürdüğü enerji ihmal edilebilir.
| Bölüm | Soru | Puan |
|---|---|---|
| B1 | Dünya'ya ulaşan nötrinoların sayı akısı yoğunluğunu, \(\Phi_\nu\), \(\mathrm{m^{-2}\,s^{-1}}\) biriminde hesaplayınız. Yukarıdaki tepkimede açığa çıkan enerji \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\)'dür. Güneş'in yaydığı enerjinin tamamen bu tepkimeden kaynaklandığını varsayınız. | 0.6 |
Güneş'in çekirdeğinden Dünya'ya yol alırken elektron nötrinolarının, \(\nu_e\), bir kısmı başka tür nötrinolara, \(\nu_x\), dönüşür. Dedektörün \(\nu_x\)'i algılama verimi, \(\nu_e\)'yi algılama veriminin 1/6'sıdır. Nötrino dönüşümü olmasaydı, bir yılda ortalama \(N_1\) nötrino algılamayı beklerdik. Ancak dönüşüm nedeniyle gerçekte yılda ortalama \(N_2\) nötrino (\(\nu_e\) ve \(\nu_x\) birlikte) algılanır.
| Bölüm | Soru | Puan |
|---|---|---|
| B2 | \(\nu_e\)'lerin ne kadarının \(\nu_x\)'e dönüştüğünü gösteren \(f\) oranını \(N_1\) ve \(N_2\) cinsinden hesaplayınız. | 0.4 |
Nötrinoları algılamak için suyla dolu büyük dedektörler inşa edilir. Nötrinoların maddeyle etkileşimleri çok seyrek olsa da zaman zaman dedektördeki su moleküllerinden elektron koparırlar. Bu yüksek enerjili elektronlar su içinde yüksek süratlerle hareket eder ve bu sırada elektromanyetik ışınım yayarlar. Böyle bir elektronun sürati, ışığın sudaki (kırılma indisi \(n\)) süratinden büyük olduğu sürece, Çerenkov ışınımı adı verilen bu ışınım bir koni biçiminde yayılır.
| Bölüm | Soru | Puan |
|---|---|---|
| B3 | Bir nötrinonun kopardığı bir elektronun su içinde ilerlerken birim zamanda sabit \(\alpha\) miktarında enerji kaybettiğini varsayınız. Bu elektron \(\Delta t\) süresince Çerenkov ışınımı yayıyorsa, nötrinonun bu elektrona aktardığı enerjiyi (\(E_{\mathrm{imparted}}\)) \(\alpha\), \(\Delta t\), \(n\), \(m_e\) ve \(c\) cinsinden belirleyiniz. (Elektronun nötrinoyla etkileşmeden önce durgun olduğunu varsayınız.) | 2.0 |
Güneş'in içinde H'nin füzyonla He'ye dönüşmesi birkaç adımda gerçekleşir. Bu ara adımlardan birinde bir \({}^{7}\mathrm{Be}\) çekirdeği (durgun kütlesi \(m_{\mathrm{Be}}\)) üretilir. Ardından bu çekirdek bir elektron soğurarak bir \({}^{7}\mathrm{Li}\) çekirdeği (durgun kütlesi \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)) üretebilir ve bir \(\nu_e\) yayabilir. İlgili nükleer tepkime şöyledir:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
Bir Be çekirdeği (\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\)) durgunken yine durgun olan bir elektronu soğurduğunda, yayılan nötrinonun enerjisi \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\) olur. Ancak Be çekirdekleri, Güneş'in çekirdeğindeki \(T_c\) sıcaklığı nedeniyle rastgele ısıl hareket yapar ve hareketli nötrino kaynakları gibi davranır. Bunun sonucunda yayılan nötrinoların enerjisi, karekök ortalama kare (rms) değeri \(\Delta E_{\mathrm{rms}}\) olan dalgalanmalar gösterir.
| Bölüm | Soru | Puan |
|---|---|---|
| B4 | \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\) ise Be çekirdeklerinin rms süratini, \(V_{\mathrm{Be}}\), hesaplayınız ve buradan \(T_c\)'yi tahmin ediniz. (İpucu: \(\Delta E_{\mathrm{rms}}\), hızın görüş doğrultusundaki bileşeninin rms değerine bağlıdır.) | 2.0 |
Soru metni resmi İngilizce metinden çevrilmiştir.
Cevaplar
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 Eğri, başlangıç noktasından \(1/3\) eğimle başlar, \((2.27,0.457)\) yakınında tek bir maksimuma ulaşır ve yatay eksene yukarıdan yaklaşır: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{alttan yaklaşarak}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
Çözüm
Sayısal hesaplamalarda standart sabitleri \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] ve \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) ile \(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\) dönüşümlerini kullanınız.
A1
Stefan–Boltzmann yasasına göre bir kara cisim, birim yüzey alanı başına \(\sigma T^4\) gücünde ışınım yayar. Güneş'in ışınım yayan yüzeyinin alanı \(4\pi R_\odot^2\) olduğundan \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] Dolayısıyla, \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
Gelen toplam güç, tüm frekanslar üzerinden alınan integraldir: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] Şu değişken dönüşümünü yapalım: \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] Bu durumda \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] Verilen integral formülü şu sonucu verir: \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] Böylece \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] Aşağıdaki tüm verim hesaplarında, Wien yaklaşımıyla bulunan bu aynı giriş gücü kullanılır.
A3
Frekansı \(\nu\) olan her foton \(h\nu\) enerjisini taşır. Spektral gücü bu enerjiye bölmek, spektral foton geliş hızını verir: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] Buradan \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] Birimi, saniye ve hertz başına foton sayısıdır.
A4
Yalnızca \(\nu\geq\nu_g\) koşulunu sağlayan fotonlar katkıda bulunur ve her biri \(E_g\) kadar yararlı enerji sağlar. Dolayısıyla, \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] \(E_g=x_gk_BT_s\) bağıntısını ve A2'deki değişken dönüşümünü kullanarak \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] Şimdi \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] Dolayısıyla, \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
Verim, yararlı çıkış gücünün gelen güce oranıdır: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] Ortak boyutlu ve geometrik çarpanlar sadeleşir ve \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \] elde edilir.
A6
Bant aralığı sıfır olduğunda, \[ \boxed{\eta(0)=0}. \] Bu durumda fotonlar elektronları uyarabilse de elektron başına yararlı enerji sıfırdır.
Büyük \(x_g\) değerleri için, \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] dolayısıyla \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] Bant aralığının üzerindeki foton sayısının üstel olarak küçülmesi, yararlı çıkışın sıfıra gitmesine neden olur.
Türev alındığında, \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] Böylece \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{alttan yaklaşarak}}. \]
Grafik, \(\eta=x_g/3\) teğetiyle başlangıç noktasından başlar, \(x_g=2.27\) yakınında tek bir maksimuma yükselir ve ardından pozitif kalarak \(\eta=0\) yatay asimptotuna doğru azalır. Eğrinin biçimi hakkında ek bilgi için, \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] olduğundan eğri \(0<x_g<4\) için aşağı doğru, \(x_g>4\) için yukarı doğru içbükeydir.
Temsilî noktalar şunlardır:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
Maksimumda \(\eta'(x_0)=0\) olur. Üstel çarpan hiçbir zaman sıfır olmadığından, \[ -x_0^3+x_0^2+2x_0+2=0, \] ya da \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
\(x_0>0\) için bu denklem \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] ile eşdeğerdir. Sol taraf \(x_0\) ile kesin olarak artarken sağ taraf kesin olarak azalır. Dolayısıyla yalnızca bir pozitif durağan nokta vardır. A6'daki uç nokta davranışıyla birlikte bu, söz konusu noktanın mutlak maksimum olduğunu gösterir.
\(F(x)=x^3-x^2-2x-2\) olsun. Bu durumda \[ F(2)=-2,\qquad F(2.5)=2.375. \] Böylece \(2<x_0<2.5\) olur; bu da zaten \(x_0=2.25\pm0.25\) sonucunu verir. Daha hassas hesaplandığında \[ \boxed{x_0\simeq2.27}. \] Verim ifadesinde yerine koyarak, \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] Buradan \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
A1'de bulunan yüzey sıcaklığında, \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] Silisyum için, \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] Dolayısıyla, \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] Böylece \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] Silisyumun bant aralığı, en uygun bant aralığı olan \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\) değerine yakın olduğu için bu sonuç maksimuma çok yakındır.
A9
Sonsuzdan art arda ince küresel kabuklar getirerek düzgün yoğunluklu küreyi oluşturalım. \(r\) yarıçapında, o ana kadar bir araya getirilmiş kütle \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] ve kabuğun kütlesi \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] Kabuk teoremine göre eklenen potansiyel enerji \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] Dolayısıyla, \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] Buradan \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] Eksi işareti, sonsuz uzaklıkta potansiyel enerjinin sıfır alınmasına karşılık gelir.
A10
Başlangıçtaki çok geniş bulutun kütleçekimsel potansiyel enerjisi yaklaşık sıfırdır. Dolayısıyla büzülmesi en fazla \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] kadar enerji açığa çıkarır. Mümkün olan en uzun parlama süresi, açığa çıkan bu enerjinin tümünün ışınımla yayıldığı kabul edilerek bulunur. Sabit ışınım gücünde, \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] dolayısıyla \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] Sayısal olarak, \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] Açığa çıkan enerjinin bir kısmının iç enerji olarak tutulması, kütleçekimsel büzülmenin parlamayı sürdürebileceği süreyi kısaltır; dolayısıyla bu bir üst sınırdır.
B1
Nötrinoların taşıyıp götürdüğü enerji ihmal edildiğinde saniyede gerçekleşen füzyon tepkimesi sayısı \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] Her tepkime iki nötrino üretir: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] Küresel yayılımda bu nötrinolar, Dünya'nın uzaklığında \(4\pi d_\odot^2\) büyüklüğünde bir alandan geçer. Böylece \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] Değerler yerine konulduğunda \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
Dönüşümden sonra \(1-f\) oranındaki kısım elektron nötrinosu olarak kalır ve yılda \((1-f)N_1\) algılamaya katkıda bulunur. Dönüşen kısmın algılanma verimi altıda bir olduğundan bu kısmın katkısı \(fN_1/6\) olur.
Dolayısıyla, \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] Çözüldüğünde, \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] Fiziksel olarak mümkün olan \(0\leq f\leq1\) aralığı, \(N_1/6\leq N_2\leq N_1\) aralığına karşılık gelir.
B3
Elektron yavaşlayarak ışığın sudaki süratine ulaştığında Çerenkov ışınımı sona erer: \[ v_{\mathrm{th}}=\frac{c}{n}. \] Elektronun o andaki kinetik enerjisi görelilik bağıntısıyla verilir: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
Elektron, \(\Delta t\) ışınım yayma süresi boyunca \(\alpha\Delta t\) kadar enerji kaybeder. Dolayısıyla nötrino etkileşiminin hemen ardından sahip olduğu başlangıç kinetik enerjisi \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] Elektron etkileşimden önce başlangıçta durgun olduğundan, nötrinonun aktardığı enerji \(K_i\)'ye eşittir. Buradan \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] İkinci terim zorunludur: Çerenkov ışınımı durduğunda elektronun hâlâ kinetik enerjisi vardır.
B4
Enerji dağılımındaki genişlemenin Be kaynaklarının ısıl hareketinden kaynaklandığı, soruda belirtilen modeli kullanalım. Kaynağın hızının Dünya'ya doğru bileşeni \(v_\parallel\) olsun.
Kütlesiz bir nötrino için kaynağın durgunluk çerçevesindeki momentum \(p_\nu=E_\nu/c\)'dir. Göreli olmayan kaynak hızının birinci mertebesine kadar, enerjinin Lorentz dönüşümü \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] sonucunu verir. Böylece Doppler enerji kayması \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
Isıl hareket izotropiktir, dolayısıyla \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] Enerji kaymasının rms değeri alındığında, \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] Dolayısıyla, \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] Sayısal olarak, \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] ve buradan \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) oranı, göreli olmayan yaklaşımın geçerliliğini doğrular.
Üç öteleme serbestlik derecesi için enerjinin eşit paylaşımı ilkesine göre, \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] Böylece \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] Veriler yerine konulduğunda, \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] buradan \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
Çözüm, resmi çözümle karşılaştırılarak kontrol edildi.
IPhO 2015 · Problema teórico 1
Partículas procedentes del Sol
Los fotones procedentes de la superficie del Sol y los neutrinos procedentes de su núcleo pueden proporcionarnos información sobre las temperaturas solares y también confirmar que el Sol brilla gracias a las reacciones nucleares.
En todo este problema, tome la masa del Sol como \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\), su radio como \(R_\odot=7.00\times10^8\ \mathrm{m}\), su luminosidad (energía emitida en forma de radiación por unidad de tiempo) como \(L_\odot=3.85\times10^{26}\ \mathrm{W}\) y la distancia Tierra-Sol como \(d_\odot=1.50\times10^{11}\ \mathrm{m}\).
Nota:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{constante} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{constante} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{constante} \]
A Radiación procedente del Sol :
| Parte | Pregunta | Puntos |
|---|---|---|
| A1 | Suponga que el Sol emite radiación como un cuerpo negro perfecto. Utilice este hecho para calcular la temperatura, \(T_s\), de la superficie solar. | 0.3 |
El espectro de la radiación solar puede aproximarse bien mediante la ley de distribución de Wien. Según esta ley, la energía solar que incide sobre cualquier superficie de la Tierra por unidad de tiempo y por unidad de intervalo de frecuencia, \(u(\nu)\), viene dada por
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
donde \(\nu\) es la frecuencia y \(A\) es el área de la superficie perpendicular a la dirección de la radiación incidente.
Considere ahora una célula solar que consiste en un disco delgado de material semiconductor de área \(A\), colocado perpendicularmente a la dirección de los rayos solares.
| Parte | Pregunta | Puntos |
|---|---|---|
| A2 | Utilizando la aproximación de Wien, exprese la potencia total de la radiación solar, \(P_{\mathrm{in}}\), que incide sobre la superficie de la célula solar, en función de \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) y las constantes fundamentales \(c\), \(h\), \(k_B\). | 0.3 |
| A3 | Exprese el número de fotones, \(n_\gamma(\nu)\), que inciden sobre la superficie de la célula solar por unidad de tiempo y por unidad de intervalo de frecuencia, en función de \(A\), \(R_\odot\), \(d_\odot\), \(T_s\), \(\nu\) y las constantes fundamentales \(c\), \(h\), \(k_B\). | 0.2 |
El material semiconductor de la célula solar tiene una «banda prohibida» de energía \(E_g\). Suponemos el siguiente modelo. Cada fotón de energía \(E\geq E_g\) excita un electrón, haciéndolo atravesar la banda prohibida. Este electrón aporta una energía \(E_g\) como energía útil de salida, y cualquier energía adicional se disipa en forma de calor (no se convierte en energía útil).
| Parte | Pregunta | Puntos |
|---|---|---|
| A4 | Defina \(x_g=h\nu_g/k_BT_s\), donde \(E_g=h\nu_g\). Exprese la potencia útil de salida de la célula, \(P_{\mathrm{out}}\), en función de \(x_g\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) y las constantes fundamentales \(c\), \(h\), \(k_B\). | 1.0 |
| A5 | Exprese la eficiencia, \(\eta\), de esta célula solar en función de \(x_g\). | 0.2 |
| A6 | Dibuje una gráfica cualitativa de \(\eta\) en función de \(x_g\). Deben mostrarse claramente los valores en \(x_g=0\) y cuando \(x_g\to\infty\). ¿Cuál es la pendiente de \(\eta(x_g)\) en \(x_g=0\) y cuando \(x_g\to\infty\)? | 1.0 |
| A7 | Sea \(x_0\) el valor de \(x_g\) para el cual \(\eta\) es máxima. Obtenga la ecuación cúbica que determina \(x_0\). Estime el valor de \(x_0\) con una precisión de \(\pm0.25\). A partir de él, calcule \(\eta(x_0)\). | 1.0 |
| A8 | La energía de la banda prohibida del silicio puro es \(E_g=1.11\ \mathrm{eV}\). Calcule la eficiencia, \(\eta_{\mathrm{Si}}\), de una célula solar de silicio utilizando este valor. | 0.2 |
\(^{1}\) Amol Dighe (TIFR), Anwesh Mazumdar (HBCSE-TIFR) y Vijay A. Singh (excoordinador nacional de las Olimpiadas de Ciencias) fueron los autores principales de este problema. Se agradecen las contribuciones del Comité Académico, el Grupo de Desarrollo Académico y la Junta Internacional.
A finales del siglo XIX, Kelvin y Helmholtz (KH) propusieron una hipótesis para explicar por qué brilla el Sol. Postularon que el Sol, que inicialmente era una nube de materia muy extensa de masa \(M_\odot\) y densidad despreciable, se ha ido contrayendo continuamente. El brillo del Sol se debería entonces a la liberación de energía potencial gravitatoria mediante esta lenta contracción.
| Parte | Pregunta | Puntos |
|---|---|---|
| A9 | Supongamos que la densidad de materia es uniforme en el interior del Sol. Halle la energía potencial gravitatoria total, \(\Omega\), del Sol en la actualidad, en función de \(G\), \(M_\odot\) y \(R_\odot\). | 0.3 |
| A10 | Estime el tiempo máximo posible, \(\tau_{\mathrm{KH}}\) (en años), durante el cual el Sol podría haber estado brillando, según la hipótesis KH. Suponga que la luminosidad del Sol ha sido constante durante todo este periodo. | 0.5 |
El valor de \(\tau_{\mathrm{KH}}\) calculado anteriormente no coincide con la edad del sistema solar estimada a partir de estudios de meteoritos. Esto demuestra que la fuente de energía del Sol no puede ser puramente gravitatoria.
B Neutrinos procedentes del Sol :
En 1938, Hans Bethe propuso que la fusión nuclear del hidrógeno en helio en el núcleo del Sol es la fuente de su energía. La reacción nuclear neta es:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
Los «neutrinos electrónicos», \(\nu_e\), producidos en esta reacción pueden considerarse de masa nula. Escapan del Sol y su detección en la Tierra confirma que se producen reacciones nucleares en el interior del Sol. En este problema puede despreciarse la energía que se llevan los neutrinos.
| Parte | Pregunta | Puntos |
|---|---|---|
| B1 | Calcule la densidad de flujo del número de neutrinos que llegan a la Tierra, \(\Phi_\nu\), en unidades de \(\mathrm{m^{-2}\,s^{-1}}\). La energía liberada en la reacción anterior es \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\). Suponga que la energía radiada por el Sol se debe íntegramente a esta reacción. | 0.6 |
Durante su viaje desde el núcleo del Sol hasta la Tierra, algunos de los neutrinos electrónicos, \(\nu_e\), se convierten en otros tipos de neutrinos, \(\nu_x\). La eficiencia del detector para detectar \(\nu_x\) es 1/6 de su eficiencia para detectar \(\nu_e\). Si no hubiera conversión de neutrinos, esperaríamos detectar una media de \(N_1\) neutrinos en un año. Sin embargo, debido a la conversión, se detecta realmente una media de \(N_2\) neutrinos (\(\nu_e\) y \(\nu_x\) conjuntamente) al año.
| Parte | Pregunta | Puntos |
|---|---|---|
| B2 | En función de \(N_1\) y \(N_2\), calcule qué fracción, \(f\), de los \(\nu_e\) se convierte en \(\nu_x\). | 0.4 |
Para detectar neutrinos, se construyen grandes detectores llenos de agua. Aunque las interacciones de los neutrinos con la materia son muy poco frecuentes, en ocasiones arrancan electrones de las moléculas de agua del detector. Estos electrones energéticos se desplazan por el agua a gran velocidad y emiten radiación electromagnética durante el proceso. Mientras la rapidez de uno de estos electrones sea mayor que la velocidad de la luz en el agua (índice de refracción \(n\)), esta radiación, llamada radiación de Cherenkov, se emite en forma de cono.
| Parte | Pregunta | Puntos |
|---|---|---|
| B3 | Suponga que un electrón arrancado por un neutrino pierde energía a un ritmo constante de \(\alpha\) por unidad de tiempo mientras se desplaza por el agua. Si este electrón emite radiación de Cherenkov durante un tiempo \(\Delta t\), determine la energía transferida a este electrón (\(E_{\mathrm{imparted}}\)) por el neutrino, en función de \(\alpha\), \(\Delta t\), \(n\), \(m_e\) y \(c\). (Suponga que el electrón estaba en reposo antes de su interacción con el neutrino.) | 2.0 |
La fusión de H en He en el interior del Sol tiene lugar en varias etapas. En una de estas etapas intermedias se produce un núcleo de \({}^{7}\mathrm{Be}\) (masa en reposo \(m_{\mathrm{Be}}\)). Posteriormente, este puede absorber un electrón, produciendo un núcleo de \({}^{7}\mathrm{Li}\) (masa en reposo \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)) y emitiendo un \(\nu_e\). La reacción nuclear correspondiente es:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
Cuando un núcleo de Be (\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\)) está en reposo y absorbe un electrón que también está en reposo, el neutrino emitido tiene una energía \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\). Sin embargo, los núcleos de Be presentan un movimiento térmico aleatorio debido a la temperatura \(T_c\) en el núcleo del Sol, y actúan como fuentes de neutrinos en movimiento. Como resultado, la energía de los neutrinos emitidos fluctúa con un valor cuadrático medio (rms) \(\Delta E_{\mathrm{rms}}\).
| Parte | Pregunta | Puntos |
|---|---|---|
| B4 | Si \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\), calcule la rapidez cuadrática media de los núcleos de Be, \(V_{\mathrm{Be}}\), y, a partir de ella, estime \(T_c\). (Sugerencia: \(\Delta E_{\mathrm{rms}}\) depende del valor cuadrático medio de la componente de la velocidad a lo largo de la línea de visión). | 2.0 |
Enunciado traducido del texto oficial en inglés.
Respuestas
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 La curva parte del origen con pendiente \(1/3\), alcanza un único máximo cerca de \((2.27,0.457)\) y se aproxima al eje horizontal por arriba: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{por debajo}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
Solución
Para los cálculos numéricos, utilice las constantes estándar \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] con \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) y \(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\).
A1
Según la ley de Stefan–Boltzmann, un cuerpo negro emite una potencia \(\sigma T^4\) por unidad de superficie. El área emisora del Sol es \(4\pi R_\odot^2\), por lo que \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] Por tanto, \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
La potencia incidente total es la integral sobre todas las frecuencias: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] Hacemos el cambio \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] Entonces \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] La fórmula de integración proporcionada da \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] Así, \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] Todas las eficiencias calculadas a continuación utilizan esta misma potencia de entrada obtenida mediante la aproximación de Wien.
A3
Cada fotón de frecuencia \(\nu\) transporta una energía \(h\nu\). Al dividir la potencia espectral entre esta energía se obtiene la tasa espectral de llegada de fotones: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] Por tanto, \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] Sus unidades son fotones por segundo y por hercio.
A4
Solo contribuyen los fotones con \(\nu\geq\nu_g\), y cada uno aporta una energía útil \(E_g\). En consecuencia, \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] Utilizando \(E_g=x_gk_BT_s\) y el cambio de variable de A2, \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] Ahora, \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] Por tanto, \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
La eficiencia es la potencia útil de salida dividida entre la potencia incidente: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] Los factores dimensionales y geométricos comunes se cancelan, lo que da \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \]
A6
Cuando la energía de la banda prohibida es nula, \[ \boxed{\eta(0)=0}. \] Aunque en ese caso los fotones pueden excitar electrones, la energía útil por electrón es cero.
Para valores grandes de \(x_g\), \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] por lo que \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] El número exponencialmente pequeño de fotones con energía superior a la de la banda prohibida hace que la potencia útil de salida tienda a cero.
Derivando, \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] Así, \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{por debajo}}. \]
La gráfica parte del origen con tangente \(\eta=x_g/3\), crece hasta un único máximo cerca de \(x_g=2.27\) y después decrece hacia la asíntota horizontal \(\eta=0\), manteniéndose positiva. Para obtener información adicional sobre su forma, \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] por lo que la curva es cóncava hacia abajo para \(0<x_g<4\) y cóncava hacia arriba para \(x_g>4\).
Algunos puntos representativos son:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
En el máximo, \(\eta'(x_0)=0\). Como el factor exponencial nunca se anula, \[ -x_0^3+x_0^2+2x_0+2=0, \] o bien \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
Para \(x_0>0\), esto equivale a \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] El miembro izquierdo crece estrictamente con \(x_0\), mientras que el miembro derecho decrece estrictamente. Por tanto, solo existe un punto estacionario positivo. Junto con el comportamiento en los extremos descrito en A6, esto demuestra que se trata del máximo global.
Sea \(F(x)=x^3-x^2-2x-2\). Entonces \[ F(2)=-2,\qquad F(2.5)=2.375. \] Así, \(2<x_0<2.5\), lo que ya proporciona \(x_0=2.25\pm0.25\). Al refinar el cálculo se obtiene \[ \boxed{x_0\simeq2.27}. \] Sustituyendo en la expresión de la eficiencia, \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] Por tanto, \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
A la temperatura superficial hallada en A1, \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] Para el silicio, \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] Por tanto, \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] Así, \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] Este valor está muy cerca del máximo porque la energía de la banda prohibida del silicio está próxima a la energía óptima \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\).
A9
Construimos la esfera uniforme trayendo sucesivas capas esféricas delgadas desde el infinito. A un radio \(r\), la masa ya reunida es \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] y la masa de la capa es \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] Según el teorema de las capas esféricas, la energía potencial añadida es \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] Por tanto, \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] De ahí, \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] El signo negativo corresponde a tomar la energía potencial igual a cero para una separación infinita.
A10
La nube inicial, muy extensa, tiene una energía potencial gravitatoria aproximadamente nula. Por tanto, su contracción libera como máximo \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] El tiempo máximo posible durante el cual el Sol puede brillar se obtiene suponiendo que toda esta energía liberada se ha radiado. A luminosidad constante, \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] por lo que \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] Numéricamente, \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] Cualquier energía liberada que se retuviera como energía interna reduciría la duración que puede sostener la contracción gravitatoria, de modo que este resultado es una cota superior.
B1
Despreciando la energía que se llevan los neutrinos, el número de reacciones de fusión por segundo es \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] Cada reacción produce dos neutrinos: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] Para una emisión esférica, estos neutrinos atraviesan un área \(4\pi d_\odot^2\) a la distancia de la Tierra. Así, \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] Al sustituir se obtiene \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
Tras la conversión, la fracción \(1-f\) permanece como neutrinos electrónicos y aporta \((1-f)N_1\) detecciones al año. La fracción convertida aporta \(fN_1/6\), porque su eficiencia de detección es seis veces menor.
En consecuencia, \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] Despejando, \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] El intervalo físicamente admisible \(0\leq f\leq1\) corresponde a \(N_1/6\leq N_2\leq N_1\).
B3
La radiación de Cherenkov cesa cuando el electrón se frena hasta alcanzar la velocidad de la luz en el agua: \[ v_{\mathrm{th}}=\frac{c}{n}. \] La energía cinética del electrón en ese instante es relativista: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
Durante el tiempo de emisión \(\Delta t\), el electrón pierde una energía \(\alpha\Delta t\). Por tanto, su energía cinética inicial, inmediatamente después de la interacción con el neutrino, es \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] Como el electrón estaba inicialmente en reposo antes de la interacción, la energía transferida por el neutrino es igual a \(K_i\). Por tanto, \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] El segundo término es esencial: el electrón todavía tiene energía cinética cuando cesa la emisión de Cherenkov.
B4
Utilizamos el modelo indicado, en el que la dispersión de energía procede del movimiento térmico de las fuentes de Be. Sea \(v_\parallel\) la componente de la velocidad de la fuente dirigida hacia la Tierra.
Para un neutrino de masa nula, el momento en el sistema de referencia en reposo de la fuente es \(p_\nu=E_\nu/c\). A primer orden en la velocidad no relativista de la fuente, la transformación de Lorentz de la energía da \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] Por tanto, el desplazamiento Doppler de la energía es \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
El movimiento térmico es isótropo, por lo que \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] Tomando el valor cuadrático medio del desplazamiento de energía, \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] Por tanto, \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] Numéricamente, \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] y, por tanto, \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] El cociente \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) justifica la aproximación no relativista.
Por el principio de equipartición aplicado a los tres grados de libertad de traslación, \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] Así, \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] Sustituyendo los datos, \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] lo que da \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
La solución se comprobó con la solución oficial.
IPhO 2015 · 理论题 1
来自太阳的粒子
来自太阳表面的光子和来自太阳核心的中微子能够告诉我们太阳的温度,还能证实太阳发光是核反应的结果。
在本题中,取太阳质量为 \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\),半径为 \(R_\odot=7.00\times10^8\ \mathrm{m}\),光度(单位时间内辐射的能量)为 \(L_\odot=3.85\times10^{26}\ \mathrm{W}\),日地距离为 \(d_\odot=1.50\times10^{11}\ \mathrm{m}\)。
提示:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{常数} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{常数} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{常数} \]
A 太阳辐射:
| 部分 | 问题 | 分值 |
|---|---|---|
| A1 | 假设太阳像理想黑体一样辐射。利用这一点计算太阳表面的温度 \(T_s\)。 | 0.3 |
太阳辐射的光谱可以用维恩分布定律很好地近似。因此,单位时间、单位频率间隔内入射到地球上任意表面的太阳能 \(u(\nu)\) 为
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
其中,\(\nu\) 为频率,\(A\) 为垂直于入射辐射方向的表面面积。
现在考虑一个太阳能电池,它由面积为 \(A\) 的半导体薄圆片构成,放置方向垂直于太阳光线。
| 部分 | 问题 | 分值 |
|---|---|---|
| A2 | 利用维恩近似,用 \(A\)、\(R_\odot\)、\(d_\odot\)、\(T_s\) 以及基本常数 \(c\)、\(h\)、\(k_B\) 表示入射到太阳能电池表面的太阳辐射总功率 \(P_{\mathrm{in}}\)。 | 0.3 |
| A3 | 用 \(A\)、\(R_\odot\)、\(d_\odot\)、\(T_s\)、\(\nu\) 以及基本常数 \(c\)、\(h\)、\(k_B\) 表示单位时间、单位频率间隔内入射到太阳能电池表面的光子数 \(n_\gamma(\nu)\)。 | 0.2 |
太阳能电池的半导体材料具有能量宽度为 \(E_g\) 的“禁带”。我们采用以下模型:每个能量为 \(E\geq E_g\) 的光子都会激发一个电子,使其跨越禁带。这个电子贡献的有用输出能量为 \(E_g\),多余的能量则以热的形式耗散(不转化为有用能量)。
| 部分 | 问题 | 分值 |
|---|---|---|
| A4 | 定义 \(x_g=h\nu_g/k_BT_s\),其中 \(E_g=h\nu_g\)。用 \(x_g\)、\(A\)、\(R_\odot\)、\(d_\odot\)、\(T_s\) 以及基本常数 \(c\)、\(h\)、\(k_B\) 表示电池的有用输出功率 \(P_{\mathrm{out}}\)。 | 1.0 |
| A5 | 用 \(x_g\) 表示该太阳能电池的效率 \(\eta\)。 | 0.2 |
| A6 | 定性画出 \(\eta\) 随 \(x_g\) 变化的示意图。应清楚标出 \(x_g=0\) 和 \(x_g\to\infty\) 时的值。\(\eta(x_g)\) 在 \(x_g=0\) 和 \(x_g\to\infty\) 时的斜率是多少? | 1.0 |
| A7 | 设 \(x_0\) 为使 \(\eta\) 达到最大值的 \(x_g\) 值。求出确定 \(x_0\) 的三次方程。估算 \(x_0\) 的值,精度要求为 \(\pm0.25\)。由此计算 \(\eta(x_0)\)。 | 1.0 |
| A8 | 纯硅的禁带宽度为 \(E_g=1.11\ \mathrm{eV}\)。利用此值计算硅太阳能电池的效率 \(\eta_{\mathrm{Si}}\)。 | 0.2 |
\(^{1}\) 本题的主要命题人为 Amol Dighe(TIFR)、Anwesh Mazumdar(HBCSE-TIFR)和 Vijay A. Singh(前科学奥林匹克竞赛全国协调员)。谨向学术委员会、学术发展小组和国际委员会的贡献致谢。
十九世纪末,开尔文和亥姆霍兹(KH)提出了一个解释太阳如何发光的假说。他们假定,太阳最初是一个质量为 \(M_\odot\)、密度可忽略不计的巨大物质云,此后一直不断收缩。于是,太阳发光便可归因于这种缓慢收缩所释放的引力势能。
| 部分 | 问题 | 分值 |
|---|---|---|
| A9 | 假设太阳内部的物质密度均匀。用 \(G\)、\(M_\odot\) 和 \(R_\odot\) 求出太阳目前的总引力势能 \(\Omega\)。 | 0.3 |
| A10 | 根据 KH 假说,估算太阳可能已经发光的最长时间 \(\tau_{\mathrm{KH}}\)(以年为单位)。假设在此期间太阳的光度始终不变。 | 0.5 |
上面算出的 \(\tau_{\mathrm{KH}}\) 与通过陨石研究估算的太阳系年龄不符。这表明太阳的能量来源不可能仅仅是引力。
B 来自太阳的中微子:
1938 年,汉斯·贝特提出,太阳核心中氢聚变为氦的核反应是太阳的能量来源。净核反应为:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
此反应产生的“电子中微子”\(\nu_e\) 可视为无质量粒子。它们逃离太阳,在地球上探测到它们便证实了太阳内部发生核反应。本题中可忽略中微子带走的能量。
| 部分 | 问题 | 分值 |
|---|---|---|
| B1 | 计算到达地球的中微子数通量密度 \(\Phi_\nu\),单位为 \(\mathrm{m^{-2}\,s^{-1}}\)。上述反应释放的能量为 \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\)。假设太阳辐射的能量完全来自这一反应。 | 0.6 |
从太阳核心传播到地球的过程中,一部分电子中微子 \(\nu_e\) 转化为其他类型的中微子 \(\nu_x\)。探测器对 \(\nu_x\) 的探测效率是对 \(\nu_e\) 的探测效率的 1/6。如果没有中微子转化,我们预期每年平均探测到 \(N_1\) 个中微子。然而,由于发生转化,实际每年平均探测到的中微子数为 \(N_2\)(包括 \(\nu_e\) 和 \(\nu_x\))。
| 部分 | 问题 | 分值 |
|---|---|---|
| B2 | 用 \(N_1\) 和 \(N_2\) 计算 \(\nu_e\) 转化为 \(\nu_x\) 的比例 \(f\)。 | 0.4 |
为了探测中微子,人们建造了装满水的大型探测器。虽然中微子与物质的相互作用极为罕见,但它们偶尔会从探测器内的水分子中击出电子。这些高能电子在水中高速运动,并在此过程中发射电磁辐射。只要这种电子的速度大于水中光速(水的折射率为 \(n\)),就会以锥形发射这种称为切伦科夫辐射的辐射。
| 部分 | 问题 | 分值 |
|---|---|---|
| B3 | 假设被中微子击出的电子在水中运动时,以每单位时间 \(\alpha\) 的恒定速率损失能量。如果该电子发射切伦科夫辐射的时间为 \(\Delta t\),用 \(\alpha\)、\(\Delta t\)、\(n\)、\(m_e\) 和 \(c\) 确定中微子传递给该电子的能量(\(E_{\mathrm{imparted}}\))。(假设电子在与中微子相互作用之前处于静止状态。) | 2.0 |
太阳内部 H 聚变为 He 的过程分为若干步骤。其中一个中间步骤会产生 \({}^{7}\mathrm{Be}\) 原子核(静质量为 \(m_{\mathrm{Be}}\))。随后,它可以吸收一个电子,生成一个 \({}^{7}\mathrm{Li}\) 原子核(静质量为 \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)),并发射一个 \(\nu_e\)。相应的核反应为:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
当一个 Be 原子核(\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\))处于静止状态,并吸收一个同样处于静止状态的电子时,所发射中微子的能量为 \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\)。然而,由于太阳核心的温度为 \(T_c\),Be 原子核处于无规则热运动中,因而成为运动的中微子源。由此,发射出的中微子能量会发生涨落,其涨落的均方根(rms)值为 \(\Delta E_{\mathrm{rms}}\)。
| 部分 | 问题 | 分值 |
|---|---|---|
| B4 | 若 \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\),计算 Be 原子核的方均根速率 \(V_{\mathrm{Be}}\),并由此估算 \(T_c\)。(提示:\(\Delta E_{\mathrm{rms}}\) 取决于沿视线方向的速度分量的均方根值。) | 2.0 |
题目译自官方英文原文。
答案
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 曲线从原点出发,初始斜率为 \(1/3\),在 \((2.27,0.457)\) 附近达到唯一的最大值,随后从上方趋近横轴: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{从负值一侧趋近}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
解答
数值计算采用以下标准常数 \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] 以及 \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) 和 \(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\)。
A1
根据斯特藩–玻尔兹曼定律,黑体单位表面积的辐射功率为 \(\sigma T^4\)。太阳的辐射面积为 \(4\pi R_\odot^2\),因此 \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] 所以, \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
总入射功率由对所有频率积分得到: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] 令 \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] 则 \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] 由题目给出的积分公式可得 \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] 因此 \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] 下文的所有效率均使用这一维恩近似下的输入功率。
A3
每个频率为 \(\nu\) 的光子携带的能量为 \(h\nu\)。将光谱功率除以此能量,即得单位频率间隔内的光子到达率: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] 所以 \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] 其单位为每秒每赫兹的光子数。
A4
只有 \(\nu\geq\nu_g\) 的光子才有贡献,每个光子贡献的有用能量为 \(E_g\)。因此, \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] 利用 \(E_g=x_gk_BT_s\) 及 A2 中的变量代换, \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] 又有 \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] 所以, \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
效率等于有用输出功率与入射功率之比: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] 约去共同的有量纲因子和几何因子,得到 \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \]
A6
当禁带宽度为零时, \[ \boxed{\eta(0)=0}. \] 此时光子虽然能够激发电子,但每个电子提供的有用能量为零。
当 \(x_g\) 很大时, \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] 因此 \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] 能量超过禁带宽度的光子数按指数规律变得极小,使有用输出趋于零。
求导得 \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] 因此 \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{从负值一侧趋近}}. \]
示意曲线从原点出发,初始切线为 \(\eta=x_g/3\),上升至 \(x_g=2.27\) 附近的唯一最大值,随后下降,趋向水平渐近线 \(\eta=0\),且始终为正。为了进一步确定曲线形状,求得 \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] 因此,曲线在 \(0<x_g<4\) 时向下凹,在 \(x_g>4\) 时向上凹。
一些代表性点如下:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
在最大值处,\(\eta'(x_0)=0\)。由于指数因子始终不为零, \[ -x_0^3+x_0^2+2x_0+2=0, \] 即 \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
对于 \(x_0>0\),此式等价于 \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] 左端随 \(x_0\) 严格递增,右端则严格递减。因此,正半轴上只有一个驻点。结合 A6 中的端点行为,可知该点对应全局最大值。
令 \(F(x)=x^3-x^2-2x-2\)。则 \[ F(2)=-2,\qquad F(2.5)=2.375. \] 因此 \(2<x_0<2.5\),这已给出 \(x_0=2.25\pm0.25\)。进一步细化可得 \[ \boxed{x_0\simeq2.27}. \] 代入效率表达式, \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] 所以 \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
在 A1 求得的表面温度下, \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] 对于硅, \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] 所以, \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] 因此 \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] 这一数值非常接近最大效率,因为硅的禁带宽度接近最优禁带宽度 \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\)。
A9
将薄球壳从无穷远处逐层移来,构成均匀球体。在半径为 \(r\) 处,已经聚集的质量为 \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] 球壳的质量为 \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] 根据球壳定理,增加的势能为 \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] 所以, \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] 于是 \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] 负号对应于将各部分相距无穷远时的势能取为零。
A10
最初极为弥散的物质云,其引力势能近似为零。因此,它收缩时释放的能量至多为 \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] 若假设释放的能量全部以辐射形式发出,即可得到最长的可能发光时间。在光度恒定时, \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] 因此 \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] 数值上, \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] 释放的能量中,任何以内能形式保留下来的部分都会缩短引力收缩能够维持发光的时间,因此这是一个上限。
B1
忽略中微子带走的能量,每秒发生的聚变反应数为 \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] 每次反应产生两个中微子: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] 对于球对称发射,这些中微子在地球所在距离处穿过面积为 \(4\pi d_\odot^2\) 的球面。因此 \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] 代入可得 \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
转化之后,比例为 \(1-f\) 的部分仍为电子中微子,每年贡献 \((1-f)N_1\) 次探测计数。发生转化的部分贡献 \(fN_1/6\) 次计数,因为其探测效率仅为原来的六分之一。
因此, \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] 解得 \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] 物理上允许的范围 \(0\leq f\leq1\) 对应于 \(N_1/6\leq N_2\leq N_1\)。
B3
当电子减速至水中光速时,切伦科夫辐射停止: \[ v_{\mathrm{th}}=\frac{c}{n}. \] 此时电子的动能需用相对论表达式计算: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
在发射辐射的时间 \(\Delta t\) 内,电子损失的能量为 \(\alpha\Delta t\)。因此,电子刚与中微子相互作用后的初始动能为 \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] 由于电子在相互作用之前处于静止状态,中微子传递给它的能量等于 \(K_i\)。所以 \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] 第二项不可缺少:切伦科夫辐射停止时,电子仍具有动能。
B4
采用题目给定的模型,即能量展宽源于作为中微子源的 Be 原子核的热运动。设 \(v_\parallel\) 为源朝向地球的速度分量。
对于无质量中微子,在源的静止参考系中,其动量为 \(p_\nu=E_\nu/c\)。对非相对论的源速度保留至一阶,能量的洛伦兹变换给出 \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] 因此,多普勒能量偏移为 \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
热运动是各向同性的,因此 \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] 取能量偏移的均方根,得到 \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] 所以, \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] 数值上, \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] 由此可得 \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] 比值 \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) 证实了非相对论近似的有效性。
根据三个平动自由度的能量均分定理, \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] 因此 \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] 代入数据, \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] 得到 \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
解答已与官方解答核对。
IPhO 2015 · Теоретическая задача 1
Частицы от Солнца
Фотоны с поверхности Солнца и нейтрино из его ядра могут рассказать нам о температурах Солнца, а также подтвердить, что Солнце светит благодаря ядерным реакциям.
Во всей задаче считайте массу Солнца равной \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\), его радиус — \(R_\odot=7.00\times10^8\ \mathrm{m}\), его светимость (энергию излучения, испускаемую за единицу времени) — \(L_\odot=3.85\times10^{26}\ \mathrm{W}\), а расстояние от Земли до Солнца — \(d_\odot=1.50\times10^{11}\ \mathrm{m}\).
Примечание:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{постоянная} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{постоянная} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{постоянная} \]
A Излучение Солнца :
| Часть | Вопрос | Баллы |
|---|---|---|
| A1 | Предположите, что Солнце излучает как абсолютно чёрное тело. Используя это предположение, вычислите температуру \(T_s\) поверхности Солнца. | 0.3 |
Спектр солнечного излучения хорошо описывается законом распределения Вина. Согласно этому закону, энергия солнечного излучения, падающая на некоторую поверхность на Земле за единицу времени в единичном интервале частот, \(u(\nu)\), задаётся выражением
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
где \(\nu\) — частота, а \(A\) — площадь поверхности, перпендикулярной направлению падающего излучения.
Теперь рассмотрим солнечный элемент, представляющий собой тонкий диск из полупроводникового материала площадью \(A\), расположенный перпендикулярно направлению солнечных лучей.
| Часть | Вопрос | Баллы |
|---|---|---|
| A2 | Используя приближение Вина, выразите полную мощность солнечного излучения \(P_{\mathrm{in}}\), падающего на поверхность солнечного элемента, через \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) и фундаментальные постоянные \(c\), \(h\), \(k_B\). | 0.3 |
| A3 | Выразите число фотонов \(n_\gamma(\nu)\), падающих на поверхность солнечного элемента за единицу времени в единичном интервале частот, через \(A\), \(R_\odot\), \(d_\odot\), \(T_s\), \(\nu\) и фундаментальные постоянные \(c\), \(h\), \(k_B\). | 0.2 |
Полупроводниковый материал солнечного элемента имеет запрещённую зону шириной \(E_g\). Примем следующую модель. Каждый фотон с энергией \(E\geq E_g\) возбуждает электрон, переводя его через запрещённую зону. Этот электрон даёт полезную выходную энергию \(E_g\), а вся избыточная энергия рассеивается в виде тепла (не преобразуется в полезную энергию).
| Часть | Вопрос | Баллы |
|---|---|---|
| A4 | Введём \(x_g=h\nu_g/k_BT_s\), где \(E_g=h\nu_g\). Выразите полезную выходную мощность элемента \(P_{\mathrm{out}}\) через \(x_g\), \(A\), \(R_\odot\), \(d_\odot\), \(T_s\) и фундаментальные постоянные \(c\), \(h\), \(k_B\). | 1.0 |
| A5 | Выразите КПД \(\eta\) этого солнечного элемента через \(x_g\). | 0.2 |
| A6 | Изобразите качественный график зависимости \(\eta\) от \(x_g\). Значения при \(x_g=0\) и \(x_g\to\infty\) должны быть ясно указаны. Каков наклон графика \(\eta(x_g)\) при \(x_g=0\) и \(x_g\to\infty\)? | 1.0 |
| A7 | Пусть \(x_0\) — значение \(x_g\), при котором \(\eta\) максимален. Получите кубическое уравнение, определяющее \(x_0\). Оцените значение \(x_0\) с точностью \(\pm0.25\). Затем вычислите \(\eta(x_0)\). | 1.0 |
| A8 | Ширина запрещённой зоны чистого кремния равна \(E_g=1.11\ \mathrm{eV}\). Используя это значение, вычислите КПД \(\eta_{\mathrm{Si}}\) кремниевого солнечного элемента. | 0.2 |
\(^{1}\) Основными авторами этой задачи были Амол Диге (TIFR), Анвеш Мазумдар (HBCSE-TIFR) и Виджай А. Сингх (бывший национальный координатор научных олимпиад). Авторы с благодарностью отмечают вклад академического комитета, группы академического развития и международного совета.
В конце XIX века Кельвин и Гельмгольц (KH) предложили гипотезу, объясняющую, почему Солнце светит. Они предположили, что Солнце, первоначально представлявшее собой очень большое облако вещества массой \(M_\odot\) и пренебрежимо малой плотности, непрерывно сжималось. Тогда свечение Солнца было бы обусловлено высвобождением гравитационной потенциальной энергии при этом медленном сжатии.
| Часть | Вопрос | Баллы |
|---|---|---|
| A9 | Предположим, что плотность вещества внутри Солнца однородна. Найдите полную гравитационную потенциальную энергию \(\Omega\) Солнца в настоящее время через \(G\), \(M_\odot\) и \(R_\odot\). | 0.3 |
| A10 | Оцените максимально возможное время \(\tau_{\mathrm{KH}}\) (в годах), в течение которого Солнце могло светить согласно гипотезе KH. Считайте, что светимость Солнца оставалась постоянной на протяжении всего этого времени. | 0.5 |
Вычисленное выше значение \(\tau_{\mathrm{KH}}\) не соответствует возрасту Солнечной системы, оценённому по исследованиям метеоритов. Это показывает, что источник энергии Солнца не может быть исключительно гравитационным.
B Нейтрино от Солнца :
В 1938 году Ганс Бете предположил, что источником энергии Солнца является ядерный синтез гелия из водорода в его ядре. Суммарная ядерная реакция имеет вид:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
«Электронные нейтрино» \(\nu_e\), образующиеся в этой реакции, можно считать безмассовыми. Они покидают Солнце, и их регистрация на Земле подтверждает протекание ядерных реакций внутри Солнца. В этой задаче энергией, уносимой нейтрино, можно пренебречь.
| Часть | Вопрос | Баллы |
|---|---|---|
| B1 | Вычислите плотность потока числа нейтрино \(\Phi_\nu\), достигающих Земли, в единицах \(\mathrm{m^{-2}\,s^{-1}}\). Энергия, выделяющаяся в приведённой выше реакции, равна \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\). Считайте, что вся энергия, излучаемая Солнцем, обусловлена этой реакцией. | 0.6 |
На пути от ядра Солнца к Земле часть электронных нейтрино \(\nu_e\) превращается в нейтрино других типов \(\nu_x\). Эффективность регистрации \(\nu_x\) детектором составляет 1/6 от его эффективности регистрации \(\nu_e\). При отсутствии превращений нейтрино мы ожидали бы регистрировать в среднем \(N_1\) нейтрино за год. Однако из-за превращений фактически за год регистрируется в среднем \(N_2\) нейтрино (суммарно \(\nu_e\) и \(\nu_x\)).
| Часть | Вопрос | Баллы |
|---|---|---|
| B2 | Вычислите через \(N_1\) и \(N_2\), какая доля \(f\) нейтрино \(\nu_e\) превращается в \(\nu_x\). | 0.4 |
Для регистрации нейтрино сооружают большие детекторы, заполненные водой. Хотя взаимодействия нейтрино с веществом происходят крайне редко, иногда они выбивают электроны из молекул воды в детекторе. Эти электроны с высокой энергией движутся в воде с большими скоростями, испуская при этом электромагнитное излучение. Пока скорость такого электрона превышает скорость света в воде (показатель преломления \(n\)), это излучение, называемое излучением Черенкова, испускается в виде конуса.
| Часть | Вопрос | Баллы |
|---|---|---|
| B3 | Предположите, что электрон, выбитый нейтрино, при движении в воде теряет энергию с постоянной скоростью \(\alpha\) за единицу времени. Если этот электрон испускает излучение Черенкова в течение времени \(\Delta t\), определите энергию (\(E_{\mathrm{imparted}}\)), переданную ему нейтрино, через \(\alpha\), \(\Delta t\), \(n\), \(m_e\) и \(c\). (Считайте, что до взаимодействия с нейтрино электрон покоился.) | 2.0 |
Синтез He из H внутри Солнца происходит в несколько этапов. На одном из этих промежуточных этапов образуется ядро \({}^{7}\mathrm{Be}\) (масса покоя \(m_{\mathrm{Be}}\)). Впоследствии оно может захватить электрон, образуя ядро \({}^{7}\mathrm{Li}\) (масса покоя \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)) и испуская \(\nu_e\). Соответствующая ядерная реакция имеет вид:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
Когда покоящееся ядро Be (\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\)) захватывает также покоящийся электрон, испущенное нейтрино имеет энергию \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\). Однако из-за температуры \(T_c\) в ядре Солнца ядра Be совершают хаотическое тепловое движение и служат движущимися источниками нейтрино. В результате энергия испускаемых нейтрино флуктуирует со среднеквадратичным отклонением \(\Delta E_{\mathrm{rms}}\).
| Часть | Вопрос | Баллы |
|---|---|---|
| B4 | Если \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\), вычислите среднеквадратичную скорость ядер Be \(V_{\mathrm{Be}}\) и затем оцените \(T_c\). (Подсказка: \(\Delta E_{\mathrm{rms}}\) зависит от среднеквадратичного значения компоненты скорости вдоль луча зрения.) | 2.0 |
Условие переведено с официального английского текста.
Ответы
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 Кривая начинается в начале координат с наклоном \(1/3\), достигает единственного максимума вблизи \((2.27,0.457)\) и приближается к горизонтальной оси сверху: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{снизу}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
Решение
Для численных расчётов используем стандартные постоянные \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] а также \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) и \(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\).
A1
По закону Стефана — Больцмана абсолютно чёрное тело излучает мощность \(\sigma T^4\) с единицы площади поверхности. Площадь излучающей поверхности Солнца равна \(4\pi R_\odot^2\), поэтому \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] Следовательно, \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
Полная падающая мощность равна интегралу по всем частотам: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] Положим \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] Тогда \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] Приведённая в условии формула интегрирования даёт \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] Таким образом, \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] Во всех приведённых ниже расчётах КПД используется эта же входная мощность, полученная в приближении Вина.
A3
Каждый фотон частоты \(\nu\) несёт энергию \(h\nu\). Разделив спектральную мощность на эту энергию, получим число приходящих фотонов за единицу времени в единичном интервале частот: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] Отсюда \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] Эта величина измеряется в фотонах в секунду на герц.
A4
Вклад дают только фотоны с \(\nu\geq\nu_g\), причём каждый из них даёт полезную энергию \(E_g\). Следовательно, \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] Используя \(E_g=x_gk_BT_s\) и замену переменной из A2, получаем \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] При этом \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] Следовательно, \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
КПД равен отношению полезной выходной мощности к падающей мощности: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] Общие размерные и геометрические множители сокращаются, и получаем \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \]
A6
При нулевой ширине запрещённой зоны \[ \boxed{\eta(0)=0}. \] Хотя фотоны в этом случае могут возбуждать электроны, полезная энергия на один электрон равна нулю.
При больших \(x_g\) \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] поэтому \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] Экспоненциально малое число фотонов с энергией выше ширины запрещённой зоны приводит к исчезновению полезной выходной мощности.
Дифференцируя, получаем \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] Таким образом, \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{снизу}}. \]
График начинается в начале координат с касательной \(\eta=x_g/3\), возрастает до единственного максимума вблизи \(x_g=2.27\), а затем убывает к горизонтальной асимптоте \(\eta=0\), оставаясь положительным. Дополнительную информацию о форме кривой даёт \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] поэтому кривая выпукла вверх при \(0<x_g<4\) и выпукла вниз при \(x_g>4\).
Характерные точки:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
В точке максимума \(\eta'(x_0)=0\). Поскольку экспоненциальный множитель нигде не обращается в нуль, \[ -x_0^3+x_0^2+2x_0+2=0, \] или \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
При \(x_0>0\) это равносильно уравнению \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] Левая часть строго возрастает с \(x_0\), а правая — строго убывает. Следовательно, существует только одна положительная стационарная точка. Вместе с поведением на границах области, установленным в A6, это доказывает, что она является точкой глобального максимума.
Пусть \(F(x)=x^3-x^2-2x-2\). Тогда \[ F(2)=-2,\qquad F(2.5)=2.375. \] Таким образом, \(2<x_0<2.5\), что уже даёт \(x_0=2.25\pm0.25\). Уточнение даёт \[ \boxed{x_0\simeq2.27}. \] Подставляя в выражение для КПД, получаем \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] Следовательно, \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
При температуре поверхности, найденной в A1, \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] Для кремния \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] Следовательно, \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] Таким образом, \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] Это значение очень близко к максимальному, поскольку ширина запрещённой зоны кремния близка к оптимальной ширине \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\).
A9
Соберём однородный шар, последовательно принося из бесконечности тонкие сферические оболочки. При радиусе \(r\) уже собранная масса равна \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] а масса оболочки равна \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] По теореме о сферической оболочке добавляемая потенциальная энергия равна \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] Следовательно, \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] Отсюда \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] Отрицательный знак соответствует выбору нулевой потенциальной энергии при бесконечном удалении частей системы друг от друга.
A10
Изначальное очень протяжённое облако имеет приблизительно нулевую гравитационную потенциальную энергию. Поэтому при его сжатии высвобождается не более \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] Максимально возможное время свечения получается в предположении, что вся эта высвободившаяся энергия была излучена. При постоянной светимости \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] поэтому \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] Численно \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] Любая часть высвободившейся энергии, сохраняющаяся в виде внутренней энергии, уменьшила бы продолжительность свечения, обеспечиваемого гравитационным сжатием, поэтому это верхняя оценка.
B1
Пренебрегая энергией, уносимой нейтрино, получаем число реакций синтеза в секунду: \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] В каждой реакции образуются два нейтрино: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] При сферически симметричном испускании эти нейтрино пересекают поверхность площадью \(4\pi d_\odot^2\) на расстоянии Земли от Солнца. Таким образом, \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] Подстановка даёт \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
После превращений доля \(1-f\) остаётся электронными нейтрино и даёт \((1-f)N_1\) регистраций за год. Превратившаяся доля даёт \(fN_1/6\), поскольку эффективность её регистрации в шесть раз меньше.
Следовательно, \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] Решая уравнение, получаем \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] Физически допустимому диапазону \(0\leq f\leq1\) соответствует \(N_1/6\leq N_2\leq N_1\).
B3
Излучение Черенкова прекращается, когда электрон замедляется до скорости света в воде: \[ v_{\mathrm{th}}=\frac{c}{n}. \] Кинетическая энергия электрона в этот момент определяется релятивистским выражением: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
За время испускания \(\Delta t\) электрон теряет энергию \(\alpha\Delta t\). Поэтому его начальная кинетическая энергия сразу после взаимодействия с нейтрино равна \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] Поскольку до взаимодействия электрон первоначально покоился, энергия, переданная нейтрино, равна \(K_i\). Следовательно, \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] Второе слагаемое существенно: при прекращении черенковского излучения электрон всё ещё обладает кинетической энергией.
B4
Используем указанную модель, в которой разброс энергии обусловлен тепловым движением источников Be. Пусть \(v_\parallel\) — компонента скорости источника в направлении Земли.
Для безмассового нейтрино импульс в системе покоя источника равен \(p_\nu=E_\nu/c\). В первом порядке по нерелятивистской скорости источника преобразование Лоренца для энергии даёт \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] Таким образом, доплеровский сдвиг энергии равен \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
Тепловое движение изотропно, поэтому \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] Вычисляя среднеквадратичное значение сдвига энергии, получаем \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] Следовательно, \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] Численно \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] и, следовательно, \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] Отношение \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) подтверждает применимость нерелятивистского приближения.
По теореме о равнораспределении энергии для трёх поступательных степеней свободы \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] Таким образом, \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] Подставляя данные, получаем \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] откуда \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
Решение проверено по официальному решению.
IPhO 2015 · مسألة نظرية 1
جسيمات من الشمس
يمكن للفوتونات القادمة من سطح الشمس والنيوترينوات القادمة من لبّها أن تخبرنا عن درجات الحرارة الشمسية، وأن تؤكد أيضًا أن الشمس تشعّ بفعل التفاعلات النووية.
في جميع أجزاء هذه المسألة، اعتبر كتلة الشمس \(M_\odot=2.00\times10^{30}\ \mathrm{kg}\)، ونصف قطرها \(R_\odot=7.00\times10^8\ \mathrm{m}\)، وقدرتها الإشعاعية (الطاقة الإشعاعية المنبعثة في وحدة الزمن) \(L_\odot=3.85\times10^{26}\ \mathrm{W}\)، والمسافة بين الأرض والشمس \(d_\odot=1.50\times10^{11}\ \mathrm{m}\).
ملاحظة:
\[ \text{(i)}\quad \int xe^{ax}\,dx =\left(\frac{x}{a}-\frac{1}{a^2}\right)e^{ax} +\text{ثابت} \]
\[ \text{(ii)}\quad \int x^2e^{ax}\,dx =\left(\frac{x^2}{a}-\frac{2x}{a^2}+\frac{2}{a^3}\right)e^{ax} +\text{ثابت} \]
\[ \text{(iii)}\quad \int x^3e^{ax}\,dx =\left(\frac{x^3}{a}-\frac{3x^2}{a^2} +\frac{6x}{a^3}-\frac{6}{a^4}\right)e^{ax} +\text{ثابت} \]
A الإشعاع القادم من الشمس :
| الجزء | السؤال | الدرجات |
|---|---|---|
| A1 | افترض أن الشمس تشعّ كجسم أسود مثالي. استخدم هذه الحقيقة لحساب درجة حرارة سطح الشمس \(T_s\). | 0.3 |
يمكن تقريب طيف الإشعاع الشمسي تقريبًا جيدًا بقانون توزيع فين. ووفقًا لذلك، تُعطى الطاقة الشمسية الساقطة على أي سطح على الأرض في وحدة الزمن ولكل وحدة من مجال التردد، \(u(\nu)\)، بالعلاقة
\[ u(\nu)=A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \nu^3\exp(-h\nu/k_BT_s), \]
حيث \(\nu\) هو التردد، و\(A\) هي مساحة السطح العمودي على اتجاه الإشعاع الساقط.
لننظر الآن في خلية شمسية تتكون من قرص رقيق من مادة شبه موصلة، مساحته \(A\)، موضوع عموديًا على اتجاه أشعة الشمس.
| الجزء | السؤال | الدرجات |
|---|---|---|
| A2 | باستخدام تقريب فين، عبّر عن القدرة الإشعاعية الشمسية الكلية \(P_{\mathrm{in}}\) الساقطة على سطح الخلية الشمسية بدلالة \(A\)، و\(R_\odot\)، و\(d_\odot\)، و\(T_s\)، والثوابت الأساسية \(c\)، و\(h\)، و\(k_B\). | 0.3 |
| A3 | عبّر عن عدد الفوتونات \(n_\gamma(\nu)\) الساقطة على سطح الخلية الشمسية في وحدة الزمن ولكل وحدة من مجال التردد، بدلالة \(A\)، و\(R_\odot\)، و\(d_\odot\)، و\(T_s\)، و\(\nu\)، والثوابت الأساسية \(c\)، و\(h\)، و\(k_B\). | 0.2 |
للمادة شبه الموصلة في الخلية الشمسية «فجوة نطاق» طاقتها \(E_g\). نفترض النموذج الآتي: كل فوتون طاقته \(E\geq E_g\) يثير إلكترونًا عبر فجوة النطاق. ويسهم هذا الإلكترون بطاقة \(E_g\) بوصفها طاقة خرج مفيدة، وتتبدد أي طاقة زائدة على هيئة حرارة (ولا تتحول إلى طاقة مفيدة).
| الجزء | السؤال | الدرجات |
|---|---|---|
| A4 | عرّف \(x_g=h\nu_g/k_BT_s\)، حيث \(E_g=h\nu_g\). عبّر عن قدرة الخرج المفيدة للخلية \(P_{\mathrm{out}}\) بدلالة \(x_g\)، و\(A\)، و\(R_\odot\)، و\(d_\odot\)، و\(T_s\)، والثوابت الأساسية \(c\)، و\(h\)، و\(k_B\). | 1.0 |
| A5 | عبّر عن كفاءة هذه الخلية الشمسية \(\eta\) بدلالة \(x_g\). | 0.2 |
| A6 | ارسم رسمًا نوعيًا لـ\(\eta\) بدلالة \(x_g\). ينبغي إظهار القيمتين عند \(x_g=0\) وعندما \(x_g\to\infty\) بوضوح. ما ميل \(\eta(x_g)\) عند \(x_g=0\) وعندما \(x_g\to\infty\)؟ | 1.0 |
| A7 | لتكن \(x_0\) قيمة \(x_g\) التي تكون عندها \(\eta\) عظمى. أوجد المعادلة التكعيبية التي تحدد \(x_0\). قدّر قيمة \(x_0\) بدقة في حدود \(\pm0.25\). ومن ثم احسب \(\eta(x_0)\). | 1.0 |
| A8 | فجوة النطاق للسيليكون النقي هي \(E_g=1.11\ \mathrm{eV}\). احسب كفاءة خلية شمسية من السيليكون \(\eta_{\mathrm{Si}}\) باستخدام هذه القيمة. | 0.2 |
\(^{1}\) كان أمول ديغي (TIFR)، وأنويش مازومدار (HBCSE-TIFR)، وفيجاي أ. سينغ (المنسق الوطني السابق للأولمبيادات العلمية) المؤلفين الرئيسيين لهذه المسألة. ونتقدم بالشكر والتقدير على مساهمات اللجنة الأكاديمية، ومجموعة التطوير الأكاديمي، والمجلس الدولي.
في أواخر القرن التاسع عشر، اقترح كلفن وهلمهولتز (KH) فرضية لتفسير كيفية إشعاع الشمس. فقد افترضا أن الشمس بدأت كسحابة هائلة من المادة، كتلتها \(M_\odot\) وكثافتها مهملة، ثم أخذت تنكمش باستمرار. وعندئذ يكون إشعاع الشمس ناتجًا عن تحرر طاقة الوضع الثقالية خلال هذا الانكماش البطيء.
| الجزء | السؤال | الدرجات |
|---|---|---|
| A9 | لنفترض أن كثافة المادة منتظمة داخل الشمس. أوجد طاقة الوضع الثقالية الكلية للشمس في الوقت الحاضر \(\Omega\)، بدلالة \(G\)، و\(M_\odot\)، و\(R_\odot\). | 0.3 |
| A10 | قدّر أقصى مدة ممكنة \(\tau_{\mathrm{KH}}\) (بالسنوات) يمكن أن تكون الشمس قد ظلت تشعّ خلالها وفقًا لفرضية KH. افترض أن القدرة الإشعاعية للشمس كانت ثابتة طوال هذه المدة. | 0.5 |
لا تتوافق \(\tau_{\mathrm{KH}}\) المحسوبة أعلاه مع عمر النظام الشمسي المقدّر من دراسات النيازك. وهذا يبيّن أن مصدر طاقة الشمس لا يمكن أن يكون ثقاليًا بحتًا.
B النيوترينوات القادمة من الشمس :
في عام 1938، اقترح هانس بيته أن الاندماج النووي للهيدروجين إلى الهيليوم في لبّ الشمس هو مصدر طاقتها. والتفاعل النووي الصافي هو:
\[ 4\,{}^{1}\mathrm{H}\longrightarrow{}^{4}\mathrm{He} +2e^++2\nu_e \]
يمكن اعتبار «النيوترينوات الإلكترونية» \(\nu_e\) الناتجة في هذا التفاعل عديمة الكتلة. وهي تفلت من الشمس، ويؤكد الكشف عنها على الأرض حدوث تفاعلات نووية داخل الشمس. يمكن إهمال الطاقة التي تحملها النيوترينوات بعيدًا في هذه المسألة.
| الجزء | السؤال | الدرجات |
|---|---|---|
| B1 | احسب كثافة فيض عدد النيوترينوات الواصلة إلى الأرض \(\Phi_\nu\)، بوحدات \(\mathrm{m^{-2}\,s^{-1}}\). الطاقة المتحررة في التفاعل أعلاه هي \(\Delta E=4.0\times10^{-12}\ \mathrm{J}\). افترض أن الطاقة التي تشعّها الشمس ناتجة بالكامل عن هذا التفاعل. | 0.6 |
أثناء انتقال بعض النيوترينوات الإلكترونية \(\nu_e\) من لبّ الشمس إلى الأرض، تتحول إلى أنواع أخرى من النيوترينوات \(\nu_x\). كفاءة الكاشف في الكشف عن \(\nu_x\) تساوي 1/6 من كفاءته في الكشف عن \(\nu_e\). إذا لم يحدث تحول للنيوترينوات، نتوقع الكشف عن متوسط قدره \(N_1\) نيوترينو في السنة. ولكن بسبب التحول، يُكشف فعليًا عن متوسط قدره \(N_2\) نيوترينو (من \(\nu_e\) و\(\nu_x\) معًا) في السنة.
| الجزء | السؤال | الدرجات |
|---|---|---|
| B2 | بدلالة \(N_1\) و\(N_2\)، احسب النسبة \(f\) من \(\nu_e\) التي تتحول إلى \(\nu_x\). | 0.4 |
للكشف عن النيوترينوات، تُنشأ كواشف كبيرة مملوءة بالماء. وعلى الرغم من ندرة تفاعلات النيوترينوات مع المادة، فإنها تقتلع أحيانًا إلكترونات من جزيئات الماء في الكاشف. تتحرك هذه الإلكترونات عالية الطاقة عبر الماء بسرعات عالية، وتصدر أثناء ذلك إشعاعًا كهرومغناطيسيًا. ما دامت سرعة أحد هذه الإلكترونات أكبر من سرعة الضوء في الماء (معامل انكساره \(n\))، ينبعث هذا الإشعاع، المسمى إشعاع تشيرينكوف، على شكل مخروط.
| الجزء | السؤال | الدرجات |
|---|---|---|
| B3 | افترض أن إلكترونًا اقتلعه نيوترينو يفقد الطاقة بمعدل ثابت مقداره \(\alpha\) في وحدة الزمن أثناء انتقاله عبر الماء. إذا أصدر هذا الإلكترون إشعاع تشيرينكوف لمدة \(\Delta t\)، فحدد الطاقة التي منحها النيوترينو لهذا الإلكترون (\(E_{\mathrm{imparted}}\)) بدلالة \(\alpha\)، و\(\Delta t\)، و\(n\)، و\(m_e\)، و\(c\). (افترض أن الإلكترون كان ساكنًا قبل تفاعله مع النيوترينو.) | 2.0 |
يحدث اندماج H إلى He داخل الشمس عبر عدة خطوات. وتنتج نواة \({}^{7}\mathrm{Be}\) (كتلتها السكونية \(m_{\mathrm{Be}}\)) في إحدى هذه الخطوات الوسيطة. ويمكنها بعد ذلك أن تمتص إلكترونًا، فتنتج نواة \({}^{7}\mathrm{Li}\) (كتلتها السكونية \(m_{\mathrm{Li}}<m_{\mathrm{Be}}\)) وتصدر \(\nu_e\). والتفاعل النووي الموافق هو:
\[ {}^{7}\mathrm{Be}+e^-\longrightarrow{}^{7}\mathrm{Li}+\nu_e \]
عندما تكون نواة Be (\(m_{\mathrm{Be}}=11.65\times10^{-27}\ \mathrm{kg}\)) ساكنة وتمتص إلكترونًا ساكنًا أيضًا، تكون طاقة النيوترينو المنبعث \(E_\nu=1.44\times10^{-13}\ \mathrm{J}\). إلا أن نوى Be تتحرك حركة حرارية عشوائية بفعل درجة الحرارة \(T_c\) في لبّ الشمس، وتعمل كمصادر نيوترينوات متحركة. ونتيجة لذلك، تتذبذب طاقة النيوترينوات المنبعثة، ويكون الجذر التربيعي لمتوسط مربع هذه التذبذبات (rms) هو \(\Delta E_{\mathrm{rms}}\).
| الجزء | السؤال | الدرجات |
|---|---|---|
| B4 | إذا كان \(\Delta E_{\mathrm{rms}}=5.54\times10^{-17}\ \mathrm{J}\)، فاحسب السرعة الجذرية المتوسطة لنوى Be، أي \(V_{\mathrm{Be}}\)، ومن ثم قدّر \(T_c\). (تلميح: تعتمد \(\Delta E_{\mathrm{rms}}\) على الجذر التربيعي لمتوسط مربع مركّبة السرعة على امتداد خط النظر). | 2.0 |
نص المسألة مترجم عن النص الإنجليزي الرسمي.
الإجابات
-
A1 \[ T_s=\left(\frac{L_\odot}{4\pi R_\odot^2\sigma}\right)^{1/4} \simeq 5.76\times10^3\ \mathrm{K}. \]
-
A2 \[ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3}. \]
-
A3 \[ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)}. \]
-
A4 \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g}. \]
-
A5 \[ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6}. \]
-
A6 يبدأ المنحنى من نقطة الأصل بميل \(1/3\)، ويبلغ قيمة عظمى واحدة قرب \((2.27,0.457)\)، ويقترب من المحور الأفقي من الأعلى: \[ \eta(0)=0,\qquad \lim_{x_g\to\infty}\eta=0,\qquad \eta'(0)=\frac13,\qquad \lim_{x_g\to\infty}\eta'=0 \quad\text{من الأسفل}. \]
-
A7 \[ x_0^3-x_0^2-2x_0-2=0,\qquad x_0\simeq2.27,\qquad \eta(x_0)\simeq0.457=45.7\%. \]
-
A8 \[ x_{g,\mathrm{Si}}\simeq2.24,\qquad \eta_{\mathrm{Si}}\simeq0.457=45.7\%. \]
-
A9 \[ \Omega=-\frac{3GM_\odot^2}{5R_\odot} \simeq-2.29\times10^{41}\ \mathrm{J}. \]
-
A10 \[ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} \simeq1.88\times10^7\ \mathrm{years}. \]
-
B1 \[ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} \simeq6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}. \]
-
B2 \[ f=\frac65\left(1-\frac{N_2}{N_1}\right). \]
-
B3 \[ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right). \]
-
B4 \[ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} \simeq2.00\times10^5\ \mathrm{m\,s^{-1}}, \] \[ T_c=\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2 \simeq1.12\times10^7\ \mathrm{K}. \]
الحل
لإجراء الحسابات العددية، نستخدم الثوابت القياسية \[ \begin{aligned} \sigma&=5.6704\times10^{-8}\ \mathrm{W\,m^{-2}\,K^{-4}},\\ G&=6.6743\times10^{-11}\ \mathrm{m^3\,kg^{-1}\,s^{-2}},\\ k_B&=1.380649\times10^{-23}\ \mathrm{J\,K^{-1}},\\ c&=2.99792458\times10^8\ \mathrm{m\,s^{-1}}, \end{aligned} \] مع \(1\ \mathrm{eV}=1.602176634\times10^{-19}\ \mathrm{J}\) و\(1\ \mathrm{year}=3.15576\times10^7\ \mathrm{s}\).
A1
وفقًا لقانون ستيفان–بولتزمان، يصدر الجسم الأسود قدرة \(\sigma T^4\) لكل وحدة مساحة من سطحه. مساحة السطح المشعّ للشمس هي \(4\pi R_\odot^2\)، ولذلك \[ L_\odot=4\pi R_\odot^2\sigma T_s^4. \] ومن ثم، \[ T_s= \left[ \frac{3.85\times10^{26}} {4\pi(7.00\times10^8)^2(5.6704\times10^{-8})} \right]^{1/4} =\boxed{5.76\times10^3\ \mathrm{K}}. \]
A2
القدرة الكلية الساقطة هي التكامل على جميع الترددات: \[ P_{\mathrm{in}} =A\frac{R_\odot^2}{d_\odot^2}\frac{2\pi h}{c^2} \int_0^\infty\nu^3e^{-h\nu/(k_BT_s)}\,d\nu. \] نضع \[ x=\frac{h\nu}{k_BT_s}, \qquad d\nu=\frac{k_BT_s}{h}\,dx. \] وعندئذ \[ P_{\mathrm{in}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} \int_0^\infty x^3e^{-x}\,dx. \] تعطي صيغة التكامل المعطاة \[ \int_0^\infty x^3e^{-x}\,dx =\left[-e^{-x}(x^3+3x^2+6x+6)\right]_0^\infty =6. \] إذًا \[ \boxed{ P_{\mathrm{in}} =\frac{12\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} }. \] تستخدم جميع حسابات الكفاءة فيما يأتي قدرة الدخل نفسها المحسوبة بتقريب فين.
A3
يحمل كل فوتون تردده \(\nu\) طاقة \(h\nu\). وبقسمة القدرة الطيفية على هذه الطاقة نحصل على المعدل الطيفي لوصول الفوتونات: \[ n_\gamma(\nu)=\frac{u(\nu)}{h\nu}. \] ومن ثم \[ \boxed{ n_\gamma(\nu) =\frac{2\pi A R_\odot^2}{d_\odot^2c^2} \nu^2e^{-h\nu/(k_BT_s)} }. \] وحداته هي فوتونات في الثانية لكل هرتز.
A4
لا تسهم إلا الفوتونات التي تحقق \(\nu\geq\nu_g\)، ويسهم كل منها بطاقة مفيدة \(E_g\). وبناءً على ذلك، \[ P_{\mathrm{out}} =E_g\int_{\nu_g}^\infty n_\gamma(\nu)\,d\nu. \] باستخدام \(E_g=x_gk_BT_s\) والتعويض الوارد في A2، \[ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g\int_{x_g}^\infty x^2e^{-x}\,dx. \] والآن \[ \int_{x_g}^\infty x^2e^{-x}\,dx =\left[-(x^2+2x+2)e^{-x}\right]_{x_g}^\infty =(x_g^2+2x_g+2)e^{-x_g}. \] ومن ثم، \[ \boxed{ P_{\mathrm{out}} =\frac{2\pi A R_\odot^2(k_BT_s)^4} {d_\odot^2c^2h^3} x_g(x_g^2+2x_g+2)e^{-x_g} }. \]
A5
الكفاءة هي قدرة الخرج المفيدة مقسومة على القدرة الساقطة: \[ \eta=\frac{P_{\mathrm{out}}}{P_{\mathrm{in}}}. \] تُختزل العوامل المشتركة ذات الأبعاد الفيزيائية والعوامل الهندسية، فنحصل على \[ \boxed{ \eta(x_g)=\frac{x_g(x_g^2+2x_g+2)e^{-x_g}}{6} },\qquad x_g\geq0. \]
A6
عندما تكون فجوة النطاق صفرًا، \[ \boxed{\eta(0)=0}. \] على الرغم من أن الفوتونات تستطيع حينئذ إثارة الإلكترونات، فإن الطاقة المفيدة لكل إلكترون تساوي صفرًا.
للقيم الكبيرة من \(x_g\)، \[ \eta(x_g)\sim\frac{x_g^3e^{-x_g}}6, \] ولذلك \[ \boxed{\lim_{x_g\to\infty}\eta(x_g)=0}. \] إن عدد الفوتونات التي تتجاوز طاقتها الفجوة يتناقص أُسّيًا، مما يجعل الخرج المفيد يؤول إلى الصفر.
بالتفاضل، \[ \eta'(x_g) =\frac{e^{-x_g}}6 \left(-x_g^3+x_g^2+2x_g+2\right). \] إذًا \[ \boxed{\eta'(0)=\frac13}, \qquad \boxed{\lim_{x_g\to\infty}\eta'(x_g)=0\ \text{من الأسفل}}. \]
يبدأ الرسم من نقطة الأصل بمماس معادلته \(\eta=x_g/3\)، ويرتفع إلى قيمة عظمى وحيدة قرب \(x_g=2.27\)، ثم يتناقص مقتربًا من المقارب الأفقي \(\eta=0\)، مع بقائه موجبًا. ولمزيد من المعلومات عن شكل المنحنى، \[ \eta''(x_g)=\frac{e^{-x_g}}6x_g^2(x_g-4), \] ولذلك يكون المنحنى مقعّرًا نحو الأسفل عندما \(0<x_g<4\)، ومقعّرًا نحو الأعلى عندما \(x_g>4\).
فيما يأتي نقاط تمثيلية:
| \(x_g\) | \(\eta(x_g)\) |
|---|---|
| 0 | 0 |
| 1 | 0.307 |
| 2 | 0.451 |
| 2.27 | 0.457 |
| 4 | 0.317 |
| 6 | 0.124 |
| 10 | 0.00923 |
| \(\to\infty\) | \(\to0\) |
A7
عند القيمة العظمى، \(\eta'(x_0)=0\). وبما أن العامل الأُسّي لا ينعدم أبدًا، \[ -x_0^3+x_0^2+2x_0+2=0, \] أو \[ \boxed{x_0^3-x_0^2-2x_0-2=0}. \]
عندما \(x_0>0\)، تكافئ هذه المعادلة \[ x_0=1+\frac{2}{x_0}+\frac{2}{x_0^2}. \] يتزايد الطرف الأيسر تزايدًا صارمًا مع \(x_0\)، بينما يتناقص الطرف الأيمن تناقصًا صارمًا. لذلك لا توجد إلا نقطة سكون واحدة عند قيمة موجبة. ومع السلوك عند طرفي المجال المبين في A6، يثبت هذا أنها نقطة القيمة العظمى المطلقة.
لتكن \(F(x)=x^3-x^2-2x-2\). عندئذ \[ F(2)=-2,\qquad F(2.5)=2.375. \] إذًا \(2<x_0<2.5\)، وهو ما يعطينا بالفعل \(x_0=2.25\pm0.25\). ويعطي تحسين التقريب \[ \boxed{x_0\simeq2.27}. \] بالتعويض في الكفاءة، \[ \eta(x_0) =\frac{x_0(x_0^2+2x_0+2)e^{-x_0}}6 \simeq0.45703. \] ومن ثم \[ \boxed{\eta_{\max}\simeq45.7\%}. \]
A8
عند درجة حرارة السطح التي أوجدناها في A1، \[ k_BT_s\simeq0.4966\ \mathrm{eV}. \] بالنسبة إلى السيليكون، \[ x_{g,\mathrm{Si}} =\frac{E_g}{k_BT_s} =\frac{1.11}{0.4966} \simeq2.235. \] ومن ثم، \[ \eta_{\mathrm{Si}} =\frac{2.235\left(2.235^2+2(2.235)+2\right)e^{-2.235}}6 \simeq0.45694. \] إذًا \[ \boxed{\eta_{\mathrm{Si}}\simeq45.7\%}. \] هذه القيمة قريبة جدًا من القيمة العظمى، لأن فجوة السيليكون قريبة من الفجوة المثلى \(x_0k_BT_s\simeq1.13\ \mathrm{eV}\).
A9
نبني الكرة المتجانسة بإحضار قشور كروية رقيقة متتالية من اللانهاية. عند نصف القطر \(r\)، تكون الكتلة المجمّعة بالفعل \[ M(r)=M_\odot\frac{r^3}{R_\odot^3}, \] وتكون كتلة القشرة \[ dM=\frac{3M_\odot r^2}{R_\odot^3}\,dr. \] وفقًا لمبرهنة القشرة، تكون طاقة الوضع المضافة \[ d\Omega=-\frac{GM(r)\,dM}{r}. \] ومن ثم، \[ \Omega =-\frac{3GM_\odot^2}{R_\odot^6} \int_0^{R_\odot}r^4\,dr =-\frac{3GM_\odot^2}{R_\odot^6} \frac{R_\odot^5}{5}. \] وعليه \[ \boxed{\Omega=-\frac{3GM_\odot^2}{5R_\odot}} \simeq-2.29\times10^{41}\ \mathrm{J}. \] تقابل الإشارة السالبة اختيار طاقة الوضع مساوية للصفر عند تباعد لانهائي.
A10
تكون طاقة الوضع الثقالية للسحابة الابتدائية الشديدة الامتداد مساوية للصفر تقريبًا. ولذلك يحرر انكماشها، في أقصى تقدير، طاقة مقدارها \[ E_{\mathrm{available}}=-\Omega =\frac{3GM_\odot^2}{5R_\odot}. \] نحصل على أقصى مدة ممكنة للإشعاع بافتراض أن كل هذه الطاقة المتحررة قد أُشعّت. وعند ثبات القدرة الإشعاعية، \[ L_\odot\tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot}, \] ولذلك \[ \boxed{ \tau_{\mathrm{KH}} =\frac{3GM_\odot^2}{5R_\odot L_\odot} }. \] عدديًا، \[ \tau_{\mathrm{KH}} =\frac{2.288\times10^{41}}{3.85\times10^{26}} \simeq5.94\times10^{14}\ \mathrm{s} =\boxed{1.88\times10^7\ \mathrm{years}}. \] أي جزء من الطاقة المتحررة يُحتفظ به على هيئة طاقة داخلية سيقلل المدة التي يمكن للانكماش الثقالي أن يدعم فيها الإشعاع؛ لذا فإن هذه المدة حد أعلى.
B1
بإهمال الطاقة التي تحملها النيوترينوات بعيدًا، يكون عدد تفاعلات الاندماج في الثانية \[ \dot N_{\mathrm{reaction}}=\frac{L_\odot}{\Delta E}. \] ينتج كل تفاعل نيوترينوين: \[ \dot N_\nu=\frac{2L_\odot}{\Delta E}. \] في حالة الانبعاث الكروي، تعبر هذه النيوترينوات مساحة \(4\pi d_\odot^2\) عند بُعد الأرض. إذًا \[ \boxed{ \Phi_\nu=\frac{2L_\odot}{4\pi d_\odot^2\Delta E} }. \] بالتعويض نحصل على \[ \Phi_\nu =\frac{2(3.85\times10^{26})} {4\pi(1.50\times10^{11})^2(4.0\times10^{-12})} =\boxed{6.8\times10^{14}\ \mathrm{m^{-2}\,s^{-1}}}. \]
B2
بعد التحول، تبقى النسبة \(1-f\) على هيئة نيوترينوات إلكترونية، وتسهم في \((1-f)N_1\) عملية كشف في السنة. وتسهم النسبة المتحولة في \(fN_1/6\) عملية كشف، لأن كفاءة الكشف عنها تساوي السدس.
وبناءً على ذلك، \[ N_2=(1-f)N_1+\frac{fN_1}{6} =N_1\left(1-\frac{5f}{6}\right). \] بحل المعادلة، \[ \boxed{ f=\frac65\left(1-\frac{N_2}{N_1}\right) }. \] يقابل المجال الفيزيائي \(0\leq f\leq1\) المجال \(N_1/6\leq N_2\leq N_1\).
B3
يتوقف إشعاع تشيرينكوف عندما يتباطأ الإلكترون حتى تبلغ سرعته سرعة الضوء في الماء: \[ v_{\mathrm{th}}=\frac{c}{n}. \] تُعطى طاقة حركة الإلكترون عند تلك اللحظة بالعلاقة النسبية: \[ K_{\mathrm{th}} =(\gamma_{\mathrm{th}}-1)m_ec^2, \qquad \gamma_{\mathrm{th}} =\frac{1}{\sqrt{1-v_{\mathrm{th}}^2/c^2}} =\frac{1}{\sqrt{1-1/n^2}}. \]
خلال زمن الانبعاث \(\Delta t\)، يفقد الإلكترون طاقة \(\alpha\Delta t\). ولذلك فإن طاقة حركته الابتدائية مباشرة بعد التفاعل مع النيوترينو هي \[ K_i=K_{\mathrm{th}}+\alpha\Delta t. \] بما أن الإلكترون كان ساكنًا في البداية قبل التفاعل، فإن الطاقة التي يمنحها النيوترينو تساوي \(K_i\). ومن ثم \[ \boxed{ E_{\mathrm{imparted}} =\alpha\Delta t+ m_ec^2\left(\frac{1}{\sqrt{1-1/n^2}}-1\right) }. \] الحد الثاني أساسي: إذ لا يزال للإلكترون طاقة حركية عندما يتوقف انبعاث تشيرينكوف.
B4
نستخدم النموذج المذكور الذي ينشأ فيه تشتت الطاقة من الحركة الحرارية لمصادر Be. لتكن \(v_\parallel\) مركّبة سرعة المصدر باتجاه الأرض.
بالنسبة إلى نيوترينو عديم الكتلة، تكون كمية حركته في إطار سكون المصدر \(p_\nu=E_\nu/c\). إلى الرتبة الأولى في سرعة المصدر غير النسبية، يعطي تحويل لورنتز للطاقة \[ E_{\mathrm{obs}}\simeq E_\nu+v_\parallel\frac{E_\nu}{c}. \] إذًا، إزاحة الطاقة بفعل دوبلر هي \[ \delta E\simeq E_\nu\frac{v_\parallel}{c}. \]
الحركة الحرارية متساوية الخواص، ولذلك \[ \langle v_\parallel\rangle=0, \qquad \langle v_\parallel^2\rangle =\frac{\langle v^2\rangle}{3} =\frac{V_{\mathrm{Be}}^2}{3}. \] بأخذ الجذر التربيعي لمتوسط مربع إزاحة الطاقة، \[ \Delta E_{\mathrm{rms}} =\frac{E_\nu}{c}\sqrt{\langle v_\parallel^2\rangle} =\frac{E_\nu V_{\mathrm{Be}}}{\sqrt3\,c}. \] ومن ثم، \[ \boxed{ V_{\mathrm{Be}} =\sqrt3\,c\frac{\Delta E_{\mathrm{rms}}}{E_\nu} }. \] عدديًا، \[ \frac{\Delta E_{\mathrm{rms}}}{E_\nu} =\frac{5.54\times10^{-17}}{1.44\times10^{-13}} =3.847\times10^{-4}, \] ومن ثم \[ \boxed{ V_{\mathrm{Be}}\simeq2.00\times10^5\ \mathrm{m\,s^{-1}} }. \] تؤكد النسبة \(V_{\mathrm{Be}}/c\simeq6.66\times10^{-4}\) صلاحية التقريب غير النسبي.
وفقًا لمبدأ التوزيع المتساوي للطاقة على درجات الحرية الانتقالية الثلاث، \[ \frac12m_{\mathrm{Be}}V_{\mathrm{Be}}^2 =\frac32k_BT_c. \] إذًا \[ T_c=\frac{m_{\mathrm{Be}}V_{\mathrm{Be}}^2}{3k_B} =\frac{m_{\mathrm{Be}}c^2}{k_B} \left(\frac{\Delta E_{\mathrm{rms}}}{E_\nu}\right)^2. \] بالتعويض بالمعطيات، \[ T_c =\frac{(11.65\times10^{-27})(2.9979\times10^8)^2} {1.380649\times10^{-23}} \left(\frac{5.54\times10^{-17}}{1.44\times10^{-13}}\right)^2, \] فنحصل على \[ \boxed{T_c\simeq1.12\times10^7\ \mathrm{K}}. \]
تمت مراجعة الحل بمقارنته بالحل الرسمي.