Register to reply 
A question about Dirac delta function... 
Share this thread: 
#1
Feb313, 12:16 AM

P: 175

Hello,
Is this correct: [tex] \int [f_j(x)\delta (xx_i) f_k(x)\delta (xx_i)]dx = f_j(x_i)f_k(x_i)[/tex] If it is not, what must the left hand side look like in order to obtain the right handside, where the right hand side multiplies two constants? Thanks! 


#2
Feb313, 11:37 AM

Sci Advisor
PF Gold
P: 2,080

First the range of the integral must include the point x=x_i or all bets are off. Second, it's cleaner to use only one delta function.



#3
Feb313, 12:41 PM

Mentor
P: 18,244

Things like [itex]\delta(xx_i)\delta(xx_i)[/itex] or [itex]\delta^2(x)[/itex] don't exist (as far as I know). 


#4
Feb313, 08:14 PM

P: 175

A question about Dirac delta function...
Okay, thanks. I would be interested to know why exactly [itex]\delta^2(x)[/itex] doesn't exist. Assuming it doesn't, then how do I derive the right hand side of the equals sign? Can the two test functions simply be placed beside each other (as if the spike on the graph were transformed by the test function j and then the second test function k)?
That is, is this correct: [tex] \int^{\infty}_{\infty}[f_j(x) f_k(x)\delta (xx_i)]dx = f_j(x_i)f_k(x_i)[/tex] (Marcus: integral range is now explicit.) 


#5
Feb313, 08:32 PM

Sci Advisor
P: 829




#6
Feb313, 09:12 PM

P: 175

Thank you pwsnafu. I just have a quick follow up question: can one of the test functions be a function from x to a vector quantity? I understand that if ##f_j(x_i)## is a scalar function, say, 2, then the height of the spike on the graph doubles, from ##\infty## to ##2\infty##. But I'm not sure what happens to the spike if we multiply by a vector, is this well defined?
I'm not sure of a good example, but imagine for some reason you knew the mass (scalar) and acceleration (vector) of a Newtonian particle, and you wanted to determine its force (vector). Let's say you also wanted to show off your grasp of delta functions in the process; does this make sense: [tex] \int^{\infty}_{\infty}[m(x) \vec{a}(x)\delta (xx_i)]dx = m_i*\vec{a}_i = \vec{F}_i[/tex] ...where ##\vec{F}_i## is the force on the particle located at ##x_i##, and ##m_i*\vec{a}_i## is that particle's mass times its instantaneous acceleration? 


#7
Feb313, 09:56 PM

Sci Advisor
P: 829

We aren't doing calculus in the ordinals. ##\int_{\mathbb{R}^n} \mathbf{F}(\mathbb{x}) \, \delta(\mathbf{x}\mathbf{x}_i) \, d\mathbf{x} := \mathbf{F}(\mathbf{x}_i)## and develop the theory from scratch. ##\mathbf{R}^m## is not field for ##m>1##, so that is no longer a functional. Note that in the literature you'll sometimes see ##\delta^3(xx_i)##, they mean "the Dirac with variable in 3D space" and not "Dirac cubed". 


#8
Feb313, 11:10 PM

P: 175

[tex] \int^{\infty}_{\infty}[m(x) \vec{a}(x)\delta (xx_i)]dx = m_i*\vec{a}_i[/tex] And we also know that the right hand side is equivalent to the right result (that is ##\vec{F}_i##). I don't understand how this result is blocked if we haven't yet defined a DDF for the force. 


#9
Feb413, 12:44 AM

Sci Advisor
P: 829

A generalised function is a functional: by definition that means the codomain is the base field, i.e. ℝ^{1}. If you want other spaces you'll need to do it yourself. 


#10
Feb413, 01:18 AM

P: 260

Hello again James!
The equation I wrote for g(x) in the thread you mention is actually onedimensional, though it is easily generalized to three dimensions. It tells you the gravitational field at x due to a massdensity function [itex]\rho=\sum_i m_i~ \delta(\zeta x_i)[/itex]. I.e. it tells you the gravitational field at point x due to a finite number of point particles of mass [itex]m_i[/itex] located at the points [itex]x_i[/itex]. The analogous vector equation you're looking for is: [tex]\mathbf{g}(\mathbf{x})=G\int \sum_i \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] This tells you the gravitational field at the position vector [itex]\mathbf{x}[/itex] due to a number of point particles of mass [itex]m_i[/itex] located at the position vectors [itex]\mathbf{x}_i[/itex]. Note that this equation is well known and can be found here: http://en.wikipedia.org/wiki/Gauss's...Newton.27s_law 


#11
Feb413, 01:35 AM

Sci Advisor
P: 829

James: you need to be very careful because a "vector function" could mean ℝ^{3} → ℝ, or ℝ^{3} → ℝ^{3} depending on the authors.
elfmotat: you've bolded g, so I'm assuming it's codomain is ℝ^{3}. If that's true, I'm guessing you are applying Dirac componentwise, is that correct? 


#12
Feb413, 01:45 AM

P: 260

[tex]\int \delta (\mathbf{r}\mathbf{x}_i)\mathrm{d}^3 \mathbf{r}=\iiint \delta (xx_i) \delta (yy_i) \delta (zz_i)~\mathrm{d}x\mathrm{d}y\mathrm{d}z[/tex] Where [itex]\mathbf{r}=(x,y,z)[/itex] and [itex]\mathbf{x}_i=(x_i,y_i,z_i)[/itex]. (Though this assumes we're working in Cartesian coordinates.) 


#13
Feb413, 06:41 AM

P: 175

Thanks elfmotat! I was wondering if one just needed more spatial dimensions. I'm curious about your formulation, is it important that you have the sum exactly where you put it? For example, are these three equivalent:
[tex] \mathbf{g}(\mathbf{x})=G\int \sum_i \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] [tex] \mathbf{g}(\mathbf{x})=G \sum_i\int \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] [tex] \mathbf{g}(\mathbf{x})=G\int \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}\sum_i ~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] (Reason I ask: I've always used the second equation to represent the many particle law, and I use the third equation to show that by swapping in the sum, one can recover the form of the single particle law in a manner that is relevant to proving mass addditivity. So I'm wondering if they're both equivalent to your formulation (i.e. the first equation).) And now relating this to the inertial law, am I right in thinking that to account for acceleration in DDF form, I just need to make clear that I'm working in three dimensions: [tex] \int m(\mathbf{x}) \mathbf{a}(\mathbf{x})\delta (\mathbf{x}\mathbf{x}_i)]d^3\mathbf{x} = m_i*\mathbf{a}_i = \mathbf{F}_i[/tex] And for the total force ##F_T## on multiple particles: [tex] \sum_i \int m_i(\mathbf{x}) \mathbf{a}_i(\mathbf{x})\delta (\mathbf{x}\mathbf{x}_i)]d^3\mathbf{x} = \mathbf{F}_T[/tex] Thanks for the terminological hint pwsnafu, and sorry for the misunderstanding about F  I see what you mean now. 


#14
Feb413, 12:16 PM

P: 260




#15
Feb413, 09:23 PM

P: 175

I do appreciate you trying to help me, Elfmotat. The background to all this is somewhat involved. What follows is a trade off between brevity and coherence. I aim to meet you somewhere in the middle. Okay here goes...
The project concerns the principle of mass additivity, which states that the mass of a composite object is the sum of the masses of its elementary parts. In most textbook discussions of Newtonian mechanics, this principle is simply presupposed, despite playing important roles (e.g. in the proof that centres of mass obey Newton's laws). However, some textbooks are more rigorous and set out to prove mass additivity, to show that "mass additivity is a consequence of Newtonian laws". I have attached an example (taken from Kibble & Berkshire's 'Classical Mechanics' (p12)). The proof is only ever given in the context of inertial mass, and involves something like the following move. Where ##F_T## is the total force on a number of point masses indexed by i (due to an external source, which I do not specify, so as to simplify the equations), we begin with: [tex] \sum_i m_i\mathbf{a}_i = \mathbf{F}_T[/tex] It is then assumed that the particles indexed by i compose a composite body and that the force on the composite is equivalent to the total force on the parts (##F_T=F_c##). A physical situation is then assumed in which the accelerations of the particles indexed by i are identical, which allows the following transformation: [tex] \mathbf{a}_i\sum_i m_i = \mathbf{F}_c[/tex] It is then assumed that the acceleration of the composite ##\mathbf{a}_c## just is the acceleration of its parts so that: [tex] \mathbf{a}_c\sum_i m_i = \mathbf{F}_c[/tex] Which is taken to show that given the force and acceleration of the composite, the composite's mass must be additive, that is: [tex] \sum_i m_i = m_c[/tex] Hence mass is additive. There is an objection to this proof, which focuses on the simplification that the acceleration of the parts are all identical. One way to show that this simplification is problematic is to show that the analogous simplification required to prove gravitational mass additivity is unphysical. The analogous simplification requires, not identical accelerations, but identical positions. The proof starts with the fundamental gravitation law: [tex] \mathbf{g}(\mathbf{x})=G\sum_i \frac{\mathbf{x}\mathbf{x}_i}{\mathbf{x}\mathbf{x}_i^3}~m_i~[/tex] Now the analogue of the acceleration simplification: the positions of the particles indexed by i are identical, so that the following transformation is possible: [tex] \mathbf{g}(\mathbf{x})=G \frac{\mathbf{x}\mathbf{x}_1}{\mathbf{x}\mathbf{x}_1^3}\sum_i~m_i~[/tex] Assuming that the composite determines the same field as the parts (##\mathbf{g}(\mathbf{x})=\mathbf{g}_c(\mathbf{x})##) and that the composite is positioned where its parts are positioned, we can derive gravitational mass additivity in the same way. But, so the objection goes, the proof fails because it assumes the unphysical assumption of particle position overlap. This enables one to object to the inertial proof with the following argument: (1) If inertial mass additivity is a consequence of Newton's inertial laws then gravitational mass additivity is a consequence of Newton's gravitation law. (2) It is not the case that gravitational mass additivity is a consequence of Newton's gravitation law. (From above.) (3) Therefore, it is not the case that inertial mass additivity is a consequence of Newton's inertial laws. Okay, so we've seen (i) a standard textbook proof for mass additivity and (ii) an objection to that proof. There are surely many ways of responding to the objection. I am looking at one type of response. This response states that the simplifications are only a product of the simple formalism chosen in the proofs. The response sets out to show that if we can formulate the equations so that they can cope with composites whose parts have varying positions and accelerations then a more general proof is possible. Enter Dirac delta functions. The gravitational proof can be formulated as follows: Fundamental manyparticle law: [tex] \mathbf{g}(\mathbf{x})=G \sum_i\int \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] Recover the form of the single particle law: [tex] \mathbf{g}(\mathbf{x})=G\int \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}\sum_i ~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] Assume that the field determined by the composite is the field determined by its parts: [tex] \mathbf{g}_c(\mathbf{x})=G\int \frac{\mathbf{x}\mathbf{r}}{\mathbf{x}\mathbf{r}^3}\sum_i ~m_i~\delta(\mathbf{r}\mathbf{x}_i)~\mathrm{d}^3\mathbf{r}[/tex] Infer gravitational mass additivity. Then in the inertial case: Fundamental manyparticle law: [tex] \sum_i \int m_i(\mathbf{x}) \mathbf{a}_i(\mathbf{x})\delta (\mathbf{x}\mathbf{x}_i)d^3\mathbf{x} = \mathbf{F}_T[/tex] Recover the form of the single particle law (this involves deriving ##\mathbf{a}_c(\mathbf{x})##, the acceleration distribution of the composite, from the acceleration distributions of the parts, which I'll just assume here): [tex] \int \mathbf{a}_c(\mathbf{x})\sum_i m_i(\mathbf{x}) \delta (\mathbf{x}\mathbf{x}_i)d^3\mathbf{x} = \mathbf{F}_T[/tex] Assume that the force on the composite is the force on the parts: [tex] \int \mathbf{a}_c(\mathbf{x})\sum_i m_i(\mathbf{x}) \delta (\mathbf{x}\mathbf{x}_i)d^3\mathbf{x} = \mathbf{F}_c[/tex] Infer inertial mass additivity. Sorry about the length of the post. I suspect people are having trouble with my inertial laws because I'm not defining the external force that determines the acceleration of all the parts? I thought leaving the external force undefined simplifies the equations, but that may be the thing that's causing problems. Anyway, I'll stop now. I am very interested in your thoughts on this! 


#16
Feb413, 10:14 PM

P: 260

When talking about a "composite" body, it is assumed that there are internal forces which keep the individual particles at the same positions relative to each other. This is known as a "rigid body." So when an external force is applied to the (center of the) composite the force is "transferred" to all of the individual particles, resulting in each particle obtaining the same acceleration vector.
In the gravitational case each individual particle in a composite body certainly experiences a different gravitational force. This results in what is known as "tidal force." In most cases, it is assumed that the composite is either (1) too small to be significantly affected by tidal forces, or (2) the internal forces in the composite offset tidal forces so that the body isn't deformed. I suspect your troubles are coming from these simplifying assumptions. 


#17
Feb413, 10:44 PM

P: 175

I'm not assuming that composite bodies must be rigid bodies, so I am under no constraint, as far as I can tell, to keep them from deforming. Are you saying that this is precisely what the problem is? That is, are you saying that because I do not stipulate a mechanism that keeps the composite rigid, the total force ##F_T## on all the particles is not well defined (whether or not we use deltas)? This would locate the problem at the very first equation! (I take it there is no corresponding problem in the gravity case, given that we are only defining a potential at a point due to many particles; and not a distributed force.)



Register to reply 
Related Discussions  
Dirac Delta function question  Calculus & Beyond Homework  3  
Dirac Delta Function question(s)  General Physics  26  
Quick Question on the Dirac Delta Function  Calculus  12  
A question of Dirac Delta function  Advanced Physics Homework  11  
The Dirac delta function question  Calculus  5 