Converting the v3 to v4 conflict detection example (D3)

I am trying to convert the collision detection example https://bl.ocks.org/mbostock/3231298

up to D3 V4.

I converted the force, but I can't seem to get it to work. Can anyone point out what I am doing wrong?

thank

var width = window.innerWidth
    var height = window.innerHeight

    var nodes = d3.range(200).map(function() { return {radius: Math.random() * 12 + 4}; }),
        root = nodes[0];
    var color = d3.scaleOrdinal().range(d3.schemeCategory20)

    root.radius = 0;
    root.fixed = true;

    const forceX = d3.forceX(width / 2).strength(0.015)
    const forceY = d3.forceY(height / 2).strength(0.015)

    var force = d3.forceSimulation()
        .force('x', forceX)
        .force('y',  forceY)
        .force('charge', function(d, i) { return i ? 0 : -2000; })
        .on('tick', ticked)
        .nodes(nodes);

    // var force = d3.layout.force()
    //     .gravity(0.05)
    //     .charge(function(d, i) { return i ? 0 : -2000; })
    //     .nodes(nodes)
    //     .size([width, height]);


    // force.start();

    var svg = d3.select(ctr)
        .attr("width", width)
        .attr("height", height);

    svg.selectAll("circle")
        .data(nodes.slice(1))
        .enter().append("circle")
        .attr("r", function(d) { return d.radius; })
        .style("fill", function(d, i) { return color(i % 3); });

    function ticked(e) {
        var q = d3.quadtree(nodes),
            i = 0,
            n = nodes.length;

        while (++i < n) q.visit(collide(nodes[i]));

        svg.selectAll("circle")
            .attr("cx", function(d) { return d.x; })
            .attr("cy", function(d) { return d.y; });
    };

    svg.on("mousemove", function() {
        var p1 = d3.mouse(this);
        root.px = p1[0];
        root.py = p1[1];
        force.restart();
    });

    function collide(node) {
        var r = node.radius + 16,
            nx1 = node.x - r,
            nx2 = node.x + r,
            ny1 = node.y - r,
            ny2 = node.y + r;
        return function(quad, x1, y1, x2, y2) {
            if (quad.point && (quad.point !== node)) {
                var x = node.x - quad.point.x,
                    y = node.y - quad.point.y,
                    l = Math.sqrt(x * x + y * y),
                    r = node.radius + quad.point.radius;
                if (l < r) {
                    l = (l - r) / l * .5;
                    node.x -= x *= l;
                    node.y -= y *= l;
                    quad.point.x += x;
                    quad.point.y += y;
                }
            }
            return x1 > nx2 || x2 < nx1 || y1 > ny2 || y2 < ny1;
        };
    }