Are these distributions the same?












1












$begingroup$


Consider the following distribution, where $delta$ is the Dirac delta:



$$f(x,y)=delta(x)+delta(y).tag1$$



This can be viewed as a limit of the following sequence of smooth functions:



$$newcommand{sech}{operatorname{sech}}
%
f_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2right),tag2$$



with $atoinfty$.



Looking at the vicinity of $x=y=0$, we can see that this function always has some "extra" peak there due to the sum of the wavefronts, regardless of $a$. Suppose we remove this central peak, getting the following sequence:



$$g_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2-sech(ax)^2sech(ay)^2right).tag3$$



Will this still result in the same distribution $(1)$, or will it be different in some regards? In particular, if $(3)$ is used as (minus) potential in a Schrödinger equation, will the eigenfunctions be the same as those with $(1)$ used instead?










share|cite|improve this question









$endgroup$

















    1












    $begingroup$


    Consider the following distribution, where $delta$ is the Dirac delta:



    $$f(x,y)=delta(x)+delta(y).tag1$$



    This can be viewed as a limit of the following sequence of smooth functions:



    $$newcommand{sech}{operatorname{sech}}
    %
    f_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2right),tag2$$



    with $atoinfty$.



    Looking at the vicinity of $x=y=0$, we can see that this function always has some "extra" peak there due to the sum of the wavefronts, regardless of $a$. Suppose we remove this central peak, getting the following sequence:



    $$g_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2-sech(ax)^2sech(ay)^2right).tag3$$



    Will this still result in the same distribution $(1)$, or will it be different in some regards? In particular, if $(3)$ is used as (minus) potential in a Schrödinger equation, will the eigenfunctions be the same as those with $(1)$ used instead?










    share|cite|improve this question









    $endgroup$















      1












      1








      1





      $begingroup$


      Consider the following distribution, where $delta$ is the Dirac delta:



      $$f(x,y)=delta(x)+delta(y).tag1$$



      This can be viewed as a limit of the following sequence of smooth functions:



      $$newcommand{sech}{operatorname{sech}}
      %
      f_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2right),tag2$$



      with $atoinfty$.



      Looking at the vicinity of $x=y=0$, we can see that this function always has some "extra" peak there due to the sum of the wavefronts, regardless of $a$. Suppose we remove this central peak, getting the following sequence:



      $$g_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2-sech(ax)^2sech(ay)^2right).tag3$$



      Will this still result in the same distribution $(1)$, or will it be different in some regards? In particular, if $(3)$ is used as (minus) potential in a Schrödinger equation, will the eigenfunctions be the same as those with $(1)$ used instead?










      share|cite|improve this question









      $endgroup$




      Consider the following distribution, where $delta$ is the Dirac delta:



      $$f(x,y)=delta(x)+delta(y).tag1$$



      This can be viewed as a limit of the following sequence of smooth functions:



      $$newcommand{sech}{operatorname{sech}}
      %
      f_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2right),tag2$$



      with $atoinfty$.



      Looking at the vicinity of $x=y=0$, we can see that this function always has some "extra" peak there due to the sum of the wavefronts, regardless of $a$. Suppose we remove this central peak, getting the following sequence:



      $$g_a(x,y)=frac a2left(sech(ax)^2+sech(ay)^2-sech(ax)^2sech(ay)^2right).tag3$$



      Will this still result in the same distribution $(1)$, or will it be different in some regards? In particular, if $(3)$ is used as (minus) potential in a Schrödinger equation, will the eigenfunctions be the same as those with $(1)$ used instead?







      functional-analysis distribution-theory dirac-delta






      share|cite|improve this question













      share|cite|improve this question











      share|cite|improve this question




      share|cite|improve this question










      asked Jan 18 at 11:38









      RuslanRuslan

      3,72721633




      3,72721633






















          1 Answer
          1






          active

          oldest

          votes


















          1












          $begingroup$

          That subtraction indeed doesn't matter, because the term $frac{a}{2}operatorname{sech}(ax^2)operatorname{sech}(ay^2)$ goes to zero (weakly). If we multiplied it by $frac{a^2}{4}$ instead of $frac a2$ to properly compensate for two dimensions of stretching, we would get the identity $delta(x)delta(y)$ in the limit. It doesn't matter whether we have a double-height peak at the origin, a peak to match the ridge lines, or even a valley at that scale; all that matters is that we don't have something infinitely larger there.



          The thing is, when we're talking about Dirac deltas? It just doesn't matter how we got to them. Work with the $delta$s directly, not some arbitrary choice of functionals that converge to them.






          share|cite|improve this answer









          $endgroup$













          • $begingroup$
            It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
            $endgroup$
            – Ruslan
            Jan 18 at 12:19










          • $begingroup$
            $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
            $endgroup$
            – jmerry
            Jan 18 at 12:24










          • $begingroup$
            You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
            $endgroup$
            – Ruslan
            Jan 18 at 12:30











          Your Answer





          StackExchange.ifUsing("editor", function () {
          return StackExchange.using("mathjaxEditing", function () {
          StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
          StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
          });
          });
          }, "mathjax-editing");

          StackExchange.ready(function() {
          var channelOptions = {
          tags: "".split(" "),
          id: "69"
          };
          initTagRenderer("".split(" "), "".split(" "), channelOptions);

          StackExchange.using("externalEditor", function() {
          // Have to fire editor after snippets, if snippets enabled
          if (StackExchange.settings.snippets.snippetsEnabled) {
          StackExchange.using("snippets", function() {
          createEditor();
          });
          }
          else {
          createEditor();
          }
          });

          function createEditor() {
          StackExchange.prepareEditor({
          heartbeatType: 'answer',
          autoActivateHeartbeat: false,
          convertImagesToLinks: true,
          noModals: true,
          showLowRepImageUploadWarning: true,
          reputationToPostImages: 10,
          bindNavPrevention: true,
          postfix: "",
          imageUploader: {
          brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
          contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
          allowUrls: true
          },
          noCode: true, onDemand: true,
          discardSelector: ".discard-answer"
          ,immediatelyShowMarkdownHelp:true
          });


          }
          });














          draft saved

          draft discarded


















          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3078134%2fare-these-distributions-the-same%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown

























          1 Answer
          1






          active

          oldest

          votes








          1 Answer
          1






          active

          oldest

          votes









          active

          oldest

          votes






          active

          oldest

          votes









          1












          $begingroup$

          That subtraction indeed doesn't matter, because the term $frac{a}{2}operatorname{sech}(ax^2)operatorname{sech}(ay^2)$ goes to zero (weakly). If we multiplied it by $frac{a^2}{4}$ instead of $frac a2$ to properly compensate for two dimensions of stretching, we would get the identity $delta(x)delta(y)$ in the limit. It doesn't matter whether we have a double-height peak at the origin, a peak to match the ridge lines, or even a valley at that scale; all that matters is that we don't have something infinitely larger there.



          The thing is, when we're talking about Dirac deltas? It just doesn't matter how we got to them. Work with the $delta$s directly, not some arbitrary choice of functionals that converge to them.






          share|cite|improve this answer









          $endgroup$













          • $begingroup$
            It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
            $endgroup$
            – Ruslan
            Jan 18 at 12:19










          • $begingroup$
            $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
            $endgroup$
            – jmerry
            Jan 18 at 12:24










          • $begingroup$
            You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
            $endgroup$
            – Ruslan
            Jan 18 at 12:30
















          1












          $begingroup$

          That subtraction indeed doesn't matter, because the term $frac{a}{2}operatorname{sech}(ax^2)operatorname{sech}(ay^2)$ goes to zero (weakly). If we multiplied it by $frac{a^2}{4}$ instead of $frac a2$ to properly compensate for two dimensions of stretching, we would get the identity $delta(x)delta(y)$ in the limit. It doesn't matter whether we have a double-height peak at the origin, a peak to match the ridge lines, or even a valley at that scale; all that matters is that we don't have something infinitely larger there.



          The thing is, when we're talking about Dirac deltas? It just doesn't matter how we got to them. Work with the $delta$s directly, not some arbitrary choice of functionals that converge to them.






          share|cite|improve this answer









          $endgroup$













          • $begingroup$
            It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
            $endgroup$
            – Ruslan
            Jan 18 at 12:19










          • $begingroup$
            $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
            $endgroup$
            – jmerry
            Jan 18 at 12:24










          • $begingroup$
            You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
            $endgroup$
            – Ruslan
            Jan 18 at 12:30














          1












          1








          1





          $begingroup$

          That subtraction indeed doesn't matter, because the term $frac{a}{2}operatorname{sech}(ax^2)operatorname{sech}(ay^2)$ goes to zero (weakly). If we multiplied it by $frac{a^2}{4}$ instead of $frac a2$ to properly compensate for two dimensions of stretching, we would get the identity $delta(x)delta(y)$ in the limit. It doesn't matter whether we have a double-height peak at the origin, a peak to match the ridge lines, or even a valley at that scale; all that matters is that we don't have something infinitely larger there.



          The thing is, when we're talking about Dirac deltas? It just doesn't matter how we got to them. Work with the $delta$s directly, not some arbitrary choice of functionals that converge to them.






          share|cite|improve this answer









          $endgroup$



          That subtraction indeed doesn't matter, because the term $frac{a}{2}operatorname{sech}(ax^2)operatorname{sech}(ay^2)$ goes to zero (weakly). If we multiplied it by $frac{a^2}{4}$ instead of $frac a2$ to properly compensate for two dimensions of stretching, we would get the identity $delta(x)delta(y)$ in the limit. It doesn't matter whether we have a double-height peak at the origin, a peak to match the ridge lines, or even a valley at that scale; all that matters is that we don't have something infinitely larger there.



          The thing is, when we're talking about Dirac deltas? It just doesn't matter how we got to them. Work with the $delta$s directly, not some arbitrary choice of functionals that converge to them.







          share|cite|improve this answer












          share|cite|improve this answer



          share|cite|improve this answer










          answered Jan 18 at 12:03









          jmerryjmerry

          10.8k1225




          10.8k1225












          • $begingroup$
            It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
            $endgroup$
            – Ruslan
            Jan 18 at 12:19










          • $begingroup$
            $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
            $endgroup$
            – jmerry
            Jan 18 at 12:24










          • $begingroup$
            You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
            $endgroup$
            – Ruslan
            Jan 18 at 12:30


















          • $begingroup$
            It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
            $endgroup$
            – Ruslan
            Jan 18 at 12:19










          • $begingroup$
            $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
            $endgroup$
            – jmerry
            Jan 18 at 12:24










          • $begingroup$
            You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
            $endgroup$
            – Ruslan
            Jan 18 at 12:30
















          $begingroup$
          It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
          $endgroup$
          – Ruslan
          Jan 18 at 12:19




          $begingroup$
          It's not possible to work with $delta$s directly when doing numerics, which is where my question originates.
          $endgroup$
          – Ruslan
          Jan 18 at 12:19












          $begingroup$
          $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
          $endgroup$
          – jmerry
          Jan 18 at 12:24




          $begingroup$
          $delta$ isn't a function - it's an operator. It only really becomes meaningful when we convolve it with something; here, $f*h(x,y)=int_{mathbb{R}}h(x,v),dv+int_{mathbb{R}}h(u,y),du$. That's something you can work with.
          $endgroup$
          – jmerry
          Jan 18 at 12:24












          $begingroup$
          You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
          $endgroup$
          – Ruslan
          Jan 18 at 12:30




          $begingroup$
          You do have to approximate it somehow e.g. when it's a potential in a Schrödinger equation you're solving numerically.
          $endgroup$
          – Ruslan
          Jan 18 at 12:30


















          draft saved

          draft discarded




















































          Thanks for contributing an answer to Mathematics Stack Exchange!


          • Please be sure to answer the question. Provide details and share your research!

          But avoid



          • Asking for help, clarification, or responding to other answers.

          • Making statements based on opinion; back them up with references or personal experience.


          Use MathJax to format equations. MathJax reference.


          To learn more, see our tips on writing great answers.




          draft saved


          draft discarded














          StackExchange.ready(
          function () {
          StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmath.stackexchange.com%2fquestions%2f3078134%2fare-these-distributions-the-same%23new-answer', 'question_page');
          }
          );

          Post as a guest















          Required, but never shown





















































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown

































          Required, but never shown














          Required, but never shown












          Required, but never shown







          Required, but never shown







          Popular posts from this blog

          MongoDB - Not Authorized To Execute Command

          Npm cannot find a required file even through it is in the searched directory

          in spring boot 2.1 many test slices are not allowed anymore due to multiple @BootstrapWith