2025年大学《数学与应用数学》专业题库- 数学物理学与宇宙奥秘探寻

2025年大学《数学与应用数学》专业题库——数学物理学与宇宙奥秘探寻考试时间：______分钟总分：______分姓名：______一、设$\mathbf{r}=(x,y,z)$为三维空间中的位置矢量，$\nabla$为梯度算子，$\nabla\cdot$为散度算子，$\nabla\times$为旋度算子。证明在直角坐标系下，下列矢量恒等式成立：1.$\nabla\times(\nablaf)=\mathbf{0}$，其中$f$是标量函数。2.$\nabla\cdot(\nabla\times\mathbf{A})=0$，其中$\mathbf{A}$是矢量函数。3.$\nabla\times(\mathbf{A}\times\mathbf{B})=(\nabla\cdot\mathbf{B})\mathbf{A}-(\nabla\cdot\mathbf{A})\mathbf{B}+(\mathbf{B}\cdot\nabla)\mathbf{A}-(\mathbf{A}\cdot\nabla)\mathbf{B}$，其中$\mathbf{A}$和$\mathbf{B}$是矢量函数。二、考虑在区域$\Omega$内定义的二维标量场$u(x,y)$和矢量场$\mathbf{A}(x,y)=(A_1(x,y),A_2(x,y))$，其中$\Omega$由$x^2+y^2\leq1$定义，边界$\partial\Omega$为单位圆周。已知$u$和$\mathbf{A}$满足以下方程：1.$\Deltau=0$在$\Omega$内。2.$\nablau\cdot\mathbf{n}=A_2$在$\partial\Omega$上，其中$\mathbf{n}$是$\partial\Omega$上的外法向单位矢量。3.$\nabla\cdot\mathbf{A}=0$在$\Omega$内。4.$\mathbf{A}\cdot\mathbf{n}=0$在$\partial\Omega$上。请推导出在$\Omega$内$u$和$\mathbf{A}$满足的方程组，并说明其物理意义（例如，$u$和$\mathbf{A}$可能分别代表什么物理量）。三、薛定谔方程是量子力学描述微观粒子状态演化的基本方程。一维定态薛定谔方程为：$-\frac{\hbar^2}{2m}\frac{d^2\psi(x)}{dx^2}+V(x)\psi(x)=E\psi(x)$其中$\psi(x)$是波函数，$E$是能量，$V(x)$是势能，$m$是粒子质量，$\hbar$是约化普朗克常数。假设粒子被限制在一个无限深势阱中，势阱边界为$x=0$和$x=a$，势能$V(x)$在$0<x<a$区域内为0，在$x\leq0$或$x\geqa$区域内为无穷大。求该粒子的能级$E$和对应的本征态波函数$\psi_n(x)$。四、广义相对论描述了引力的几何理论。爱因斯坦场方程为：$R_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R+\Lambdag_{\mu\nu}=\frac{8\piG}{c^4}T_{\mu\nu}$其中$R_{\mu\nu}$是里奇曲率张量，$R$是标量曲率，$g_{\mu\nu}$是度规张量，$\Lambda$是宇宙学常数，$G$是万有引力常数，$c$是光速，$T_{\mu\nu}$是能量-动量张量。考虑一个仅受自身引力作用的、处于静态的、无限长的圆柱体。简述如何利用广义相对论来估算该圆柱体产生的引力场（引力势或时空曲率）的特性。不需要求解具体的场方程，只需说明关键的思想和步骤。五、宇宙膨胀是现代宇宙学的核心概念之一。弗里德曼方程描述了宇宙尺度的几何随时间演化的动态方程：$$\left(\frac{\dot{a}}{a}\right)^2=\frac{8\piG}{3}\rho-\frac{k}{a^2}+\frac{\Lambda}{3}$$其中$a(t)$是宇宙标度因子，$\dot{a}$是其时间导数，$\rho$是宇宙物质密度，$k$是宇宙曲率常数，$\Lambda$是宇宙学常数。方程左侧是宇宙的哈勃参数的平方，右侧包含物质密度、曲率和宇宙学常数项。讨论方程中各项的物理意义，并说明当$\rho,k,\Lambda$取不同值时，宇宙可能的命运（例如，开放宇宙、封闭宇宙、平坦宇宙或暗能量主导的加速膨胀宇宙）。六、数学在理解暗物质现象中扮演了重要角色。暗物质通过其引力效应被间接观测到。假设宇宙中存在一种简单的、不与普通物质发生电磁相互作用的暗物质，其密度$\rho_{\mathrm{DM}}(r)$仅是径向距离$r$的函数。普通物质密度$\rho_m(r)$也仅是$r$的函数。请推导出仅考虑引力作用时，暗物质晕（一个围绕星系或星系团分布的暗物质球壳或球体）的质量密度分布$\rho_{\mathrm{DM}}(r)$所需满足的方程（即引力平衡方程），并说明其物理意义。假设宇宙是静态的，可以忽略宇宙膨胀的影响。七、考虑一个由两个点电荷$q_1$和$q_2$组成的系统，它们位于空间中固定位置。根据经典电磁学，这个系统在真空中产生的总能量可以表示为：$$U=\frac{1}{4\pi\epsilon_0}\left(\frac{q_1q_2}{r_{12}}\right)$$其中$r_{12}$是两点电荷之间的距离，$\epsilon_0$是真空介电常数。现在假设该系统被放置在一个由相对介电常数为$\epsilon_r>1$的无限大均匀介质中。请计算此时系统的总静电能，并解释为什么计算结果与在真空中计算的结果不同。试卷答案一、1.证明：在直角坐标系下，$\nablaf=\left(\frac{\partialf}{\partialx},\frac{\partialf}{\partialy},\frac{\partialf}{\partialz}\right)$。则$\nabla\times(\nablaf)=\left(\frac{\partial^2f}{\partialy\partialz}-\frac{\partial^2f}{\partialz\partialy},\frac{\partial^2f}{\partialz\partialx}-\frac{\partial^2f}{\partialx\partialz},\frac{\partial^2f}{\partialx\partialy}-\frac{\partial^2f}{\partialy\partialx}\right)$。由于混合偏导数与求导次序无关，即$\frac{\partial^2f}{\partialy\partialz}=\frac{\partial^2f}{\partialz\partialy}$，$\frac{\partial^2f}{\partialz\partialx}=\frac{\partial^2f}{\partialx\partialz}$，$\frac{\partial^2f}{\partialx\partialy}=\frac{\partial^2f}{\partialy\partialx}$。因此，$\nabla\times(\nablaf)=(0,0,0)=\mathbf{0}$。2.证明：在直角坐标系下，$\nabla\times\mathbf{A}=\left(\frac{\partialA_3}{\partialy}-\frac{\partialA_2}{\partialz},\frac{\partialA_1}{\partialz}-\frac{\partialA_3}{\partialx},\frac{\partialA_2}{\partialx}-\frac{\partialA_1}{\partialy}\right)$。则$\nabla\cdot(\nabla\times\mathbf{A})=\frac{\partial}{\partialx}\left(\frac{\partialA_3}{\partialy}-\frac{\partialA_2}{\partialz}\right)+\frac{\partial}{\partialy}\left(\frac{\partialA_1}{\partialz}-\frac{\partialA_3}{\partialx}\right)+\frac{\partial}{\partialz}\left(\frac{\partialA_2}{\partialx}-\frac{\partialA_1}{\partialy}\right)$。展开后，所有各项的偏导数均为零，因为它们是二阶导数，且混合偏导数与求导次序无关。因此，$\nabla\cdot(\nabla\times\mathbf{A})=0$。3.证明：使用直角坐标系下矢量三重积的展开式：$\mathbf{C}\times(\mathbf{A}\times\mathbf{B})=(\mathbf{C}\cdot\mathbf{B})\mathbf{A}-(\mathbf{C}\cdot\mathbf{A})\mathbf{B}$。令$\mathbf{C}=\nabla$，则$\mathbf{C}\cdot\mathbf{B}=\nabla\cdot\mathbf{B}$，$\mathbf{C}\cdot\mathbf{A}=\nabla\cdot\mathbf{A}$。同时，$\mathbf{C}\times\mathbf{A}=\nabla\times\mathbf{A}$，$\mathbf{C}\times\mathbf{B}=\nabla\times\mathbf{B}$。代入展开式：$\nabla\times(\mathbf{A}\times\mathbf{B})=(\nabla\cdot\mathbf{B})\mathbf{A}-(\nabla\cdot\mathbf{A})\mathbf{B}+(\nabla\times\mathbf{B})\times\mathbf{A}-(\nabla\times\mathbf{A})\times\mathbf{B}$。现在处理$\nabla\times\mathbf{B}\times\mathbf{A}$和$\nabla\times\mathbf{A}\times\mathbf{B}$。利用矢量恒等式$\nabla\times(\mathbf{D}\times\mathbf{E})=(\mathbf{E}\cdot\nabla)\mathbf{D}-(\mathbf{D}\cdot\nabla)\mathbf{E}$：$(\nabla\times\mathbf{B})\times\mathbf{A}=(\mathbf{A}\cdot\nabla)\mathbf{B}-(\mathbf{B}\cdot\nabla)\mathbf{A}$。$(\nabla\times\mathbf{A})\times\mathbf{B}=(\mathbf{B}\cdot\nabla)\mathbf{A}-(\mathbf{A}\cdot\nabla)\mathbf{B}$。将这两个结果代入前面的式子：$\nabla\times(\mathbf{A}\times\mathbf{B})=(\nabla\cdot\mathbf{B})\mathbf{A}-(\nabla\cdot\mathbf{A})\mathbf{B}+(\mathbf{A}\cdot\nabla)\mathbf{B}-(\mathbf{B}\cdot\nabla)\mathbf{A}-(\mathbf{B}\cdot\nabla)\mathbf{A}+(\mathbf{A}\cdot\nabla)\mathbf{B}$。合并同类项，得到最终结果：$\nabla\times(\mathbf{A}\times\mathbf{B})=(\nabla\cdot\mathbf{B})\mathbf{A}-(\nabla\cdot\mathbf{A})\mathbf{B}+(\mathbf{B}\cdot\nabla)\mathbf{A}-(\mathbf{A}\cdot\nabla)\mathbf{B}$。二、推导：1.$\Deltau=0$在$\Omega$内意味着$u$是调和函数。2.$\nabla\cdot\mathbf{A}=0$在$\Omega$内意味着$\mathbf{A}$是无源场（或称为保守场的一部分）。3.在边界$\partial\Omega$上，$\nablau\cdot\mathbf{n}=A_2$是$u$的法向导数等于$\mathbf{A}$在法向的分量。4.在边界$\partial\Omega$上，$\mathbf{A}\cdot\mathbf{n}=0$是$\mathbf{A}$垂直于边界。利用格林第一恒等式$\int_{\Omega}u\Deltav\,d\Omega=\int_{\Omega}\nablau\cdot\nablav\,d\Omega-\int_{\partial\Omega}u\frac{\partialv}{\partialn}\,ds$，令$v=A_2$：$\int_{\Omega}u\DeltaA_2\,d\Omega=\int_{\Omega}\nablau\cdot\nablaA_2\,d\Omega-\int_{\partial\Omega}uA_2\,ds$。由于$\DeltaA_2=0$（因为$\nabla\cdot\mathbf{A}=0$，对时间求导或做其他运算通常不改变此性质，此处假设静态），第一项为0。$\int_{\Omega}\nablau\cdot\nablaA_2\,d\Omega=\int_{\partial\Omega}uA_2\,ds$。利用散度定理$\int_{\Omega}\nabla\cdot\mathbf{A}\,d\Omega=\int_{\partial\Omega}\mathbf{A}\cdot\mathbf{n}\,ds$，得到：$\int_{\Omega}\nabla\cdot(\mathbf{A}u)\,d\Omega-\int_{\Omega}u\nabla\cdot\mathbf{A}\,d\Omega=\int_{\partial\Omega}uA_2\,ds$。由于$\nabla\cdot\mathbf{A}=0$，第二项为0。又因为$\mathbf{A}\cdot\mathbf{n}=0$，第三项也为0。$\int_{\Omega}\nabla\cdot(\mathbf{A}u)\,d\Omega=\int_{\partial\Omega}uA_2\,ds$。根据高斯散度定理，左侧可化为边界积分：$\int_{\partial\Omega}(\mathbf{A}u)\cdot\mathbf{n}\,ds=\int_{\partial\Omega}uA_2\,ds$。由于$\mathbf{A}\cdot\mathbf{n}=0$，左侧$\int_{\partial\Omega}(\mathbf{A}u)\cdot\mathbf{n}\,ds=\int_{\partial\Omega}u(\mathbf{A}\cdot\mathbf{n})\,ds=0$。因此，$\int_{\partial\Omega}uA_2\,ds=0$。结合边界条件$\int_{\partial\Omega}uA_2\,ds=\int_{\partial\Omega}u\frac{\partialu}{\partialn}\,ds$，得到：$\int_{\partial\Omega}u\frac{\partialu}{\partialn}\,ds=0$。由于$u$和$\frac{\partialu}{\partialn}$在$\partial\Omega$上连续且非零（否则不构成物理问题），这只有在$u$在$\partial\Omega$上恒等于某个常数时才成立。设该常数为0（可以通过平移$u\tou+C$实现）。因此，$u=0$在$\partial\Omega$上。现在，考虑$u$和$\mathbf{A}$满足的方程组：从$\Deltau=0$和$\nabla\cdot\mathbf{A}=0$出发，可以构造方程。利用拉普拉斯算子$\Delta$是线性算子，可以写成$\Delta(uA_1)=(\Deltau)A_1+u\DeltaA_1=0$（因为$\Deltau=0$且$\DeltaA_1=0$）。同样，$\Delta(uA_2)=0$。利用格林第二恒等式$\int_{\Omega}(u\Deltav-v\Deltau)\,d\Omega=\int_{\partial\Omega}(u\frac{\partialv}{\partialn}-v\frac{\partialu}{\partialn})\,ds$，令$v=A_1$，$u$为任意试探函数：$\int_{\Omega}u\DeltaA_1\,d\Omega=\int_{\partial\Omega}u\frac{\partialA_1}{\partialn}\,ds$。由于$\DeltaA_1=0$，左侧为0。利用边界条件$u=0$在$\partial\Omega$，得到$\int_{\partial\Omega}u\frac{\partialA_1}{\partialn}\,ds=0$。同理，令$v=A_2$：$\int_{\Omega}u\DeltaA_2\,d\Omega=\int_{\partial\Omega}u\frac{\partialA_2}{\partialn}\,ds$。由于$\DeltaA_2=0$，左侧为0。利用边界条件$u=0$在$\partial\Omega$，得到$\int_{\partial\Omega}u\frac{\partialA_2}{\partialn}\,ds=0$。现在，考虑矢量场$\mathbf{A}$，利用矢量格林恒等式$\int_{\Omega}u(\nabla\cdot\mathbf{A})\,d\Omega=\int_{\partial\Omega}u(\mathbf{A}\cdot\mathbf{n})\,ds$。由于$\nabla\cdot\mathbf{A}=0$，左侧为0。利用边界条件$\mathbf{A}\cdot\mathbf{n}=0$在$\partial\Omega$，得到$\int_{\partial\Omega}u(\mathbf{A}\cdot\mathbf{n})\,ds=0$。这些推导似乎并未直接给出新的方程。更直接的思路是考虑将$\mathbf{A}$表示为势场的旋度：$\mathbf{A}=\nabla\times\mathbf{F}$，其中$\mathbf{F}$是某个标量势函数的旋度。则有$\nabla\cdot\mathbf{A}=\nabla\cdot(\nabla\times\mathbf{F})=0$，且$\nabla\times\mathbf{A}=\nabla\times(\nabla\times\mathbf{F})=\nabla(\nabla\cdot\mathbf{F})-\nabla^2\mathbf{F}$。要求$\mathbf{A}\cdot\mathbf{n}=0$在$\partial\Omega$，即$(\nabla\times\mathbf{F})\cdot\mathbf{n}=0$在$\partial\Omega$。这通常意味着$\nabla\cdot\mathbf{F}=0$在$\partial\Omega$（例如，如果$\mathbf{F}$在$\partial\Omega$上取常值）。现在，利用$\nabla\cdot(\nabla\times\mathbf{A})=0$：$\nabla\cdot(\nabla\times(\nabla\times\mathbf{F}))=\nabla\cdot(\nabla(\nabla\cdot\mathbf{F})-\nabla^2\mathbf{F})=\nabla^3(\nabla\cdot\mathbf{F})-\nabla\cdot(\nabla^2\mathbf{F})=0$。如果$\nabla\cdot\mathbf{F}=0$，则$\nabla^3(\nabla\cdot\mathbf{F})=0$。因此，要求$\nabla\cdot(\nabla^2\mathbf{F})=0$，即$\nabla^2(\nabla\cdot\mathbf{F})=0$。如果$\nabla\cdot\mathbf{F}=0$，则$\nabla^2(\nabla\cdot\mathbf{F})=0$自动满足。因此，$\mathbf{A}=\nabla\times\mathbf{F}$且$\nabla\cdot\mathbf{F}=0$是一个可能的解结构。此时，$\mathbf{A}$满足$\nabla\cdot\mathbf{A}=0$和$\mathbf{A}\cdot\mathbf{n}=0$在$\partial\Omega$。$u$满足$\Deltau=0$和$u=0$在$\partial\Omega$。物理意义：$u$可能代表某种标量势（如静电势），$\mathbf{A}$可能代表与某种矢量势相关的磁场分量（如磁矢势的旋度，在特定规范下）。边界条件$u=0$可能对应于势在边界上为零或无穷远处的条件，$\mathbf{A}\cdot\mathbf{n}=0$可能对应于通过边界的净磁通量为零。$\nabla\cdot\mathbf{A}=0$对应于无磁单极子。三、求解：1.边界条件：在$x=0$和$x=a$处，$\psi(0)=0$，$\psi(a)=0$。2.方程：$-\frac{\hbar^2}{2m}\frac{d^2\psi(x)}{dx^2}+0\cdot\psi(x)=E\psi(x)$，即$-\frac{\hbar^2}{2m}\frac{d^2\psi(x)}{dx^2}=E\psi(x)$。令$k^2=\frac{2mE}{\hbar^2}$，则方程变为$\frac{d^2\psi(x)}{dx^2}+k^2\psi(x)=0$。3.通解：通解为$\psi(x)=A\sin(kx)+B\cos(kx)$。4.应用边界条件：*$\psi(0)=A\sin(0)+B\cos(0)=B=0$。因此，$\psi(x)=A\sin(kx)$。*$\psi(a)=A\sin(ka)=0$。要使$\psi(x)\neq0$，需要$A\neq0$且$\sin(ka)=0$。这意味着$ka=n\pi$，其中$n$是非零整数（$n=1,2,3,\dots$）。如果$n=0$，则$k=0$，$\psi(x)=0$，无物理意义。因此，$k=\frac{n\pi}{a}$。5.能级：由$k=\frac{n\pi}{a}$和$k^2=\frac{2mE}{\hbar^2}$，得到能级$E_n$：$E_n=\frac{\hbar^2k^2}{2m}=\frac{\hbar^2}{2m}\left(\frac{n\pi}{a}\right)^2=\frac{n^2\pi^2\hbar^2}{2ma^2}$。能量是量子化的，只能取分立值$E_n$。6.波函数：对应的本征态波函数为$\psi_n(x)=A\sin\left(\frac{n\pix}{a}\right)$。由于$\psi_n(x)$需要归一化，通常取$A=\sqrt{\frac{2}{a}}$（通过计算$\int_0^a|\psi_n(x)|^2dx=1$验证）。归一化后的波函数为$\psi_n(x)=\sqrt{\frac{2}{a}}\sin\left(\frac{n\pix}{a}\right)$，$n=1,2,3,\dots$。四、推导与说明：1.引力场描述：圆柱体产生的引力场可以通过求解引力势$\Phi(\mathbf{r})$来描述。引力势满足泊松方程$\nabla^2\Phi=4\piG\rho_{\mathrm{cm}}(\mathbf{r})$，其中$\rho_{\mathrm{cm}}(\mathbf{r})$是圆柱体的质量密度（$\rho_{\mathrm{cm}}=\rho_{\mathrm{DM}}$，如果只考虑暗物质）。2.静态假设：静态意味着时空度规$g_{\mu\nu}$不随时间变化，时空度规张量$g^{\mu\nu}$的导数为零。3.度规形式：对于无限长圆柱体，其产生的引力场在柱坐标系$(t,r,\theta,z)$中，度规可以近似为$ds^2=-c^2dt^2+dr^2+r^2d\theta^2+dz^2$。更精确地，如果考虑引力透镜效应，需要使用静态度规，通常包含$g_{tt}$和$g_{rr}$项，形式可能为$ds^2=-c^2(1-2\Phi/c^2)dt^2+(1+2\Phi/c^2)dr^2+r^2d\theta^2+dz^2$，其中$\Phi$是引力势。4.场方程：在静态时空中，爱因斯坦场方程简化为$R_{\mu\nu}-\frac{1}{2}g_{\mu\nu}R=\frac{8\piG}{c^4}T_{\mu\nu}$。对于只受自身引力作用的孤立系统（如暗物质圆柱体），其能量-动量张量$T_{\mu\nu}$主要由其质量密度$\rho_{\mathrm{cm}}$和动量密度$\mathbf{p}_{\mathrm{cm}}$决定，且通常满足$\nabla\cdot\mathbf{p}_{\mathrm{cm}}=0$（对于静态物体），有时近似为$T_{\mu\nu}\approx(\rho_{\mathrm{cm}},0,0,0)$。5.引力势与时空曲率：*一种方法是直接求解简化后的爱因斯坦场方程。对于静态、轴对称（旋转对称）、可解度的圆柱体，其度规$g_{\mu\nu}$可以精确求解（如Bianchi度规或Tolman-Oppenheimer-Volkoff度规的简化形式）。一旦度规$g_{\mu\nu}$确定，可以通过计算$R_{\mu\nu}$和$R$来描述时空曲率。*另一种方法是利用弱场近似或只考虑引力势$\Phi$。在弱场近似下，度规可以写为$g_{\mu\nu}=\eta_{\mu\nu}+h_{\mu\nu}$，其中$\eta_{\mu\nu}$是Minkowski度规，$h_{\mu\nu}$是小的扰动量。爱因斯坦场方程可以近似为$\Boxh_{\mu\nu}=-\frac{8\piG}{c^4}T_{\mu\nu}$，其中$\Box=\partial^\mu\partial_\mu$。对于引力势$\Phi$，通常有$h_{tt}\approx-\Phi$，$h_{rr}\approx\Phi$。此时，场方程近似为$\nabla^2\Phi=4\piG\rho_{\mathrm{cm}}(\mathbf{r})$，即泊松方程。求解此方程即可得到引力势$\Phi$，进而可以得到引力加速度$-\nabla\Phi$。6.关键步骤：*明确物理模型：无限长圆柱形暗物质分布，密度$\rho_{\mathrm{DM}}(r)$。*选择合适的时空度规：可能是近似Minkowski度规（弱场）或精确的静态轴对称度规。*写出简化后的爱因斯坦场方程（静态、轴对称、物质为暗物质）。*如果使用弱场近似，推导出引力势$\Phi$满足的泊松方程。*求解泊松方程（或场方程）得到$\Phi(r)$。*由$\Phi(r)$计算引力场（加速度或时空曲率张量的部分分量）。五、讨论：1.各项物理意义：*$\frac{\dot{a}}{a}$：哈勃参数$H(t)$，描述宇宙膨胀速率，即单位时间内空间距离的增长率。*$\frac{8\piG}{3}\rho$：物质密度项，包括普通物质、暗物质和暗能量的总密度。$G$是引力常数，$\rho$是总物质能量密度。这个项描述了物质对时空曲率的影响（引力效应）。*$-\frac{k}{a^2}$：曲率项。$k$是宇宙曲率常数。*$k>0$：封闭宇宙，总能量密度大于临界密度，宇宙有限无界。*$k=0$：平坦宇宙，总能量密度等于临界密度，宇宙无限且平坦。*$k<0$：开放宇宙，总能量密度小于临界密度，宇宙无限且膨胀加速。*$\frac{\Lambda}{3}$：宇宙学常数项。$\Lambda$代表真空能量密度或标量场势能密度，通常被认为具有斥力，导致宇宙加速膨胀。2.宇宙命运讨论：*平坦宇宙($k=0$)：如果总能量密度$\rho=\rho_{\text{crit}}$(临界密度)，弗里德曼方程简化为$(H^2)=\frac{8\piG}{3}\rho_{\text{crit}}-\frac{\Lambda}{3}$。宇宙的膨胀会持续下去，但空间几何保持平坦。如果存在正的宇宙学常数$\Lambda>0$，则宇宙将加速膨胀，最终可能永远膨胀下去，趋向于热寂或大撕裂（取决于暗能量的具体性质）。如果$\Lambda=0$，宇宙会膨胀，但减速膨胀，最终可能停止膨胀并开始收缩（大坍缩），除非存在其他导致加速膨胀的机制。*开放宇宙($k<0$)：如果总能量密度$\rho<\rho_{\text{crit}}$且$\Lambda\geq0$。哈勃参数$H(t)$随时间$t$趋于零，但始终保持正值（$H>0$）。这意味着宇宙的膨胀速率会持续减小，但永远不会变为零。宇宙将永远膨胀下去，但膨胀会越来越慢。空间几何是双曲的。如果$\Lambda>0$，存在斥力，将导致加速膨胀。如果$\Lambda=0$，则减速膨胀。*封闭宇宙($k>0$)：如果总能量密度$\rho>\rho_{\text{crit}}$。哈勃参数$H(t)$随时间$t$减小，并在某个时间点（大坍缩时间）变为零并开始反向增大。这意味着宇宙膨胀最终会停止，并开始收缩。空间几何是球面的。如果存在正的宇宙学常数$\Lambda>0$，斥力会减缓收缩，但宇宙最终仍会坍缩。如果$\Lambda=0$，则宇宙会在达到最大尺度后开始坍缩。*暗能量主导的加速膨胀宇宙：当今观测表明，宇宙正在加速膨胀。这通常归因于暗能量的存在（$\Lambda>0$或某种具有斥力的暗能量模型）。在暗能量主导的宇宙中，弗里德曼方程中的$\frac{\Lambda}{3}$项起主导作用，使得$H^2$保持一个相对稳定的正值或随时间缓慢变化，导致宇宙加速膨胀。无论初始的$k$和$\rho$如何（只要$\rho$不是无限大），暗能量的存在都会在足够晚的时间尺度上主导宇宙的动力学，导致加速膨胀。六、推导：1.物理背景：暗物质不与电磁相互作用，其引力效应通过质量密度$\rho_{\mathrm{DM}}(r)$体现。普通物质密度$\rho_m(r)$也产生引力效应。总引力场由总质量密度$\rho_{\mathrm{tot}}(r)=\rho_{\mathrm{DM}}(r)+\rho_m(r)$决定。如果只考虑暗物质，则$\rho_{\mathrm{tot}}(r)=\rho_{\mathrm{DM}}(r)$。2.引力平衡条件：对于稳定的、不发生坍缩或无限膨胀的暗物质分布，其内部的每个质元都处于引力平衡状态。这意味着该质元受到的引力（来自内部和外部暗物质）与其自身的引力场（由其自身质量分布产生）以及可能的压力梯度力达到平衡。3.引力势与引力场：在引力作用下，引力势$\Phi(r)$满足泊松方程$\nabla^2\Phi=4\piG\rho_{\mathrm{tot}}(r)$。在仅考虑暗物质的情况下，$\nabla^2\Phi=4\piG\rho_{\mathrm{DM}}(r)$。4.引力加速度：暗物质质元受到的引力加速度$\mathbf{a}_{\mathrm{grav}}=-\nabla\Phi$。5.压力梯度力：对于非相对论性暗物质，其运动方程为$\mathbf{F}=m\mathbf{a}=m\mathbf{a}_{\mathrm{grav}}+\nablaP\mathbf{n}$，其中$P$是暗物质的压力，$\mathbf{n}$是单位法向矢量。在引力平衡状态下，假设暗物质处于静态或准静态状态，压力梯度力$\nablaP\mathbf{n}$和引力$\mathbf{F}=m\mathbf{a}_{\mathrm{grav}}$大小相等，方向相反。6.引力平衡方程推导：将引力加速度$\mathbf{a}_{\mathrm{grav}}=-\nabla\Phi$代入平衡条件$\mathbf{F}=m\mathbf{a}_{\mathrm{grav}}$，得到$m(-\nabla\Phi)=-m\nabla\Phi=\nablaP\mathbf{n}$。7.分量形式：写成分量形式，$-m\frac{\partial\Phi}{\partialr}=P\frac{\partialP}{\partialr}$（假设暗物质分布对径向对称，压力$P$仅是$r$的函数）。8.整理：$P\frac{\partialP}{\partialr}=m\frac{\partial\Phi}{\partialr}$。9.引力势$\Phi(r)$的求解：首先需要求解泊松方程$\nabla^2\Phi=4\piG\rho_{\mathrm{DM}}(r)$。这通常需要根据$\rho_{\mathrm{DM}}(r)$的具体形式（例如，指数分布、核球模型等）进行积分求解。例如，对于球对称的指数分布暗物质模型$\rho_{\mathrm{DM}}(r)=\rho_0\exp(-r/R)$，可以解出$\Phi(r)$。10.平衡方程应用：假设已解得$\Phi(r)$。将$\Phi(r)$代入平衡方程$P\frac{\partialP}{\partialr}=m\frac{\partial\Phi}{\partialr}$。11.求解压力$P(r)$：通过积分求解上述平衡方程，可以得到压力$P(r)$的表达式。例如，对于球对称模型，可以得到$P(r)$的解析解或数值解。12.物理意义：推导出的平衡方程$\nabla\cdot(P\mathbf{n})=m\nabla\Phi$表明，在引力平衡状态下，暗物质内部的径向压力梯度力与引力相等。这揭示了暗物质在引力场中的稳定分布需要满足的条件。它将暗物质密度$\rho_{\mathrm{DM}}(r)$、引力势$\Phi(r)$和压力$P(r)$联系起来，是描述暗物质稳定分布（如暗物质晕）的常用数学工具。13.简化假设：推导过程中可能做了以下简化假设：*暗物质分布球对称。*暗物质非相对论性（压力$P$小于引力项，允许推导简化）。*忽略自引力效应（对于稀疏暗物质分布）。*假设已知或能求解$\rho_{\mathrm{DM}}(r)$和$\Phi(r)$的表达式。*忽略磁场、自相互作用等复杂因素。七、计算与解释：1.计算总能量：系统的总静电能$U$由各部分电荷间的相互作用能构成。对于两个点电荷$q_1$和$q_ularsystem,thetotalelectrostaticenergyisgivenbythesumofthepotentialenergyofinteractionbetweenthechargesandtheexternalelectricfield,orsimplytheinteractionenergyintheabsenceofexternalfields.Fortwocharges$q_1$and$q_2$separatedbyadistance$r_{12}$,thepotentialenergyis:$U=\frac{1}{4\pi\epsilon_0}\frac{q_1q_2}{r_{12}}$.这是对真空中两个点电荷系统静电能的标准计算公式，来源于库仑定律和电场能量公式。它表示两个点电荷间的相互作用能。推导：真空中，两个点电荷$q_1$和$q_2$产生的电场分别为$\mathbf{E}_1$和$\mathbf{E}_2$。系统的总能量$U$可以通过计算每个电荷在另一电荷产生的电场中所具有的电势能之和得到：$U=q_1\phi_2+q_2\phi_1$，其中$\phi_1$和$\phi_2$分别是电荷$q_1$所在位置（